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Abstract 

Young massive clusters are dense aggregates of young stars that form the funda- 
mental building blocks of galaxies. Several examples exist in the Milky Way Galaxy 
and the Local Group, but they are particularly abundant in starburst and interacting 
galaxies. The few young massive clusters that are close enough to resolve are of prime 
interest for studying the stellar mass function and the ecological interplay between stel- 
lar evolution and stellar dynamics. The distant unresolved clusters may be effectively 
used to study the star-cluster mass function, and they provide excellent constraints on 
the formation mechanisms of young cluster populations. Young massive clusters are 
expected to be the nurseries for many unusual objects, including a wide range of exotic 
stars and binaries. So far only a few such objects have been found in young massive 
clusters, although their older cousins, the globular clusters, are unusually rich in stellar 
exotica. 

In this review we focus on star clusters younger than ~ 100 Myr, more than a few 
current crossing times old, and more massive than ~ 10 4 Mq, irrespective of cluster 
size or environment. We describe the global properties of the currently known young 
massive star clusters in the Local Group and beyond, and discuss the state of the art 
in observations and dynamical modeling of these systems. In order to make this review 
readable by observers, theorists, and computational astrophysicists, we also review the 
cross-disciplinary terminology. 
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1 Introduction 



Stars form in clustered environments (ILada fc Lada 20031 ). In the Milky Way Galaxy, evi- 



dence for this statement comes from the global clustering of spectral O -type stars (IParker fc Goodwin 20071 ) 

of which ~ 70% reside in young clusters or associations (IGies 1987|). and ~ 5 % of the re- 

maining field population are directly identified as runaways ( Ide Wit et al. 20051 ). Ide Wit et al. (20051 ) 
found that only ~ 4% of O-type stars can be considered as having formed outside a clustered 
environment ; further analysis has shown that a few of even those 4 % may actually be run- 



away stars (IGvaramadze &: Bomans 20081 . ISchilbach &: Roser 20081 ). further strengthening 



the case for clustered formation. Additional evidence that clusters are the primary mode of 
star formation comes from the observed formation rate of stars in embedded clusters (~3x 



10 3 M Myr 1 kpc 2 ; lLada fc Lada 2003) which is compa rable to the formation rate of field 



stars (~ 3-7 x 10 3 M ? Mvr" 1 kpc" 2 : iMiller fc Scalo 1979h. Finally, some 96% of the stars i n 
the nearby Orion B star-forming region are clustered (IClarke. Bonnell fc Hillenbrand 2000T ) . 

In nearby young starburst galaxies at least 20%, and possibly all, of the ultraviolet 
light appears to come from young star clusters (jMeurer et al. 19951 ) ; this also seems to be 
the case fo r the observed Ha and B-band l uminosities in inte racting galaxies , such as the 
Antennae (frail. Chandar fc Whitmore 20051 ) and NGC 3256 flZepf et al. 19991 ). The fossil 
record of an early episode of star formation is evidenced by the present-day population of 
old glob ular clusters, although the firs t (population III) stars seem not to have formed in 
clusters (lAbel. Bryan fc Norman 20021 ). 



1.1 Scope of this review 

In this review we focus on the young massive clusters (hereafter YMCs) found in an increasing 
number of galaxies. We adopt a deliberately broad definition of this term, concentrating on 
observations of star clusters younger than about 100 Myr and more massive than 1O 4 M . 
In addition, implicit throughout most of our discussion is the assumption that the clusters 
under study are actually bound. This may seem obvious but, as we will discuss in ^2j some 
very young objects satisfying our age and mass criteria may well be unbound, expanding 
freely into space following the loss of the intracluster gas out of which they formed. We find 
that imposing a third requirement, that the age of the cluster exceed its current dynamical 
time (the orbit time of a typical star — see §1.3.3p by a factor of a few, effectively distinguishes 
between bound clusters and unbound associations. 

Our adopted age limit is somewhat arbitrary, but represents roughly the epoch at which a 
cluster can be said to have survived its birth (phase 1, see §|3]) and the subsequent early phase 
of vigorous stellar evolution (phase 2, see §H]), and to be entering the long-term evolutionary 
phase during which its lifetime is determined principally by stellar dynamical processes and 
external influences (phase 3, see §1]). This latter, "stellar dynamical" phase generally starts 
after about 100 Myr and is not the principal focus of this review. 

The mass limit is such that lower-mass clusters are unlikely to survive for more than 
1 Gyr. Based on the lifetimes presented in §H] (in particular see Eq. [T9|) . we estimate that 
a cluster with an initial number of stars iV ~ 10 5 will survive for ~ 10 Gyr. Since young 
clusters more massive than 10 5 M Q are relatively rare, we relax our criterion to include star 
clusters with masses as low as 10 4 M Q . We place no limits on cluster size, metallicity, or 
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galactic location, for the practical reason that this would further reduce our already small 
sample of YMCs, even though it seems likely that clusters such as the Arches and Quintuplet 
systems near the Galactic center (see Tab. [2]) are likely to dissolve within a gigayear. 

Thus, any young cluster massive enough to survive for a significant fraction of a Hub- 
ble time — regardless of its current location — meets our criteria for inclusion in this review. 
Within the current context we cannot make a strong connection between cluster properties 
and environment, but from the discussion around Fig. [TD] below it is evident that environment 
has at least some influence on global cluster characteristics (see §2.4.2p . 

The masses and projected lifetimes of YMCs coincide with those of the old globular 
clusters (hereafter GCs) that populate the bulges and halos of many galaxies, including 
our own. Indeed, YMCs are sometimes referred to in the literature as "young globular 
clusters." This possible connection between YMCs and GCs offers the exciting prospect of 
studying in the local universe physical processes that may have occurred during the otherwise 
practically unobservable formative stages of the GC population^ However, although the 
characterization is highly suggestive, the extent to which today's YMCs will someday come 
to resemble GCs remains unclear. 

As we will see, the size distribution of YMCs does appear to be consistent with them 



evolving into GCs (iMafz-Apellaniz 20021 ). and we can reasonably expect that after (say) 10 
Gyr they will be roughly spherical in shape and will have surface brightness distributions 
similar in character to those of the GCs, but other properties are not so easy to assess. 
Except for a few nearby cases, such as Westerlund 1 and the Arches cluster (see Tab. [2] in 
$2]), observations of YMCs are limited to stellar masses ^ 1 M , whereas the most massive 
stars observed in GCs are less than 1M . Thus there is no guarantee that the stellar 
mass function in YMCs below ~ 1M resembles the initial mass functions of known GCs, 
although we note that the stellar mass funct ions in nearby open clusters are consistent wit h 



the distribution of low-mass stars in GCs ( IChabrier 20031 ; iBastian. Covey fc Meyer 20101 ). 



Similar uncertainties apply to cluster binary populations. The initial binary fractions inferred 
for GCs are generally small, whereas YMCs appear to be binary-rich (see §5. 1 1) . But the 
known binaries in massive clusters tend to be found among massive stars, while the binary 
properties of low-mass stars in YMCs are unknown. Again, any comparison is complicated 
by the absence of any significant overlap in the observed stellar mass spectra. 

Our objective in this review is to summarize the current state of knowledge of YMCs, to 
describe the key physical processes governing their evolution and survival, and to assess the 
extent to which we expect them to evolve into systems comparable to the old GCs observed 
today. 

In this review we will use the terms young, dense, and massive in relation to star clusters. 
Although not precise, these descriptions do have specific connotations. As already indicated, 
"young" means star clusters that are still in the early, violent mass-loss phase during the 
first 100 Myr (see §1]). "Dense" indicates that in some clusters the stars are packed together 
so closely that stellar collisions start to play an important role (see §3.4.2j) . In terms of 
the Safronov number^ young dense clusters have O ^ 10 2 . As a practical matter, we 
consider a cluster to be dense if its half-mass relaxation time (Eq. [l~6]) is less than ~ 10 s years 



a "Real" young globular clusters at z ~ 5 are expected to be ~2 magnitudes brighter than the detection 
limit of the James Webb Space Telescope. 

The Saf ronov (19691) number 9 is defined as the square of the ratio of the escape velocity from the stellar 
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(IPortegies Zwart. McMillan fc Makino 20071 ) . "Massive" indicates that we expect the cluster 
to survive for ~ 10 Gyr, into the "old globular cluster" regime. Tab. [T] summarizes the main 
parameters of the three different populations of star clusters. 
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Table 1: Comparison of fundamental parameters for star cluster families relevant to this 
review: open cluster (OC), globular cluster (GC), and young massive cluster (YMC). The 
numbers in the columns are intended to be indicative of the population and are rounded-off, 
and should be used with care, but they provide some flavor of the various cluster types. The 
second column gives cluster age, followed by the turn-off mass (in M Q ), the total cluster 
mass (in M ), the virial radius r vir (see §1.3.2p . the core density, and the metallicity. The 
last three columns give the location in the Galaxy where these clusters are found, and the 
dynamical and relaxation time scales, defined in §1.3.31 

Figs. [I] and [2] compare the distributions of massive open star clusters, YMCs, and GCs 
in the Milky Way Galaxy. The spatial distribution of the YMCs (Fig.[T]) clearly identify 
them as a disk population comparable to the open clusters, but in the mass-radius diagram 
(Fig. [2]), YMCs seem more closely related to GCs. 



1.2 Properties of Cluster Systems 

The Milky Way Galaxy contains some 150 GCs, with mass estimates ranging from ~ 10 3 M 
(for AM4, a member of the Sgr dwarf spheroidal galaxy) to 2.2 x 1O 6 M (for NGC5139, 
Omega Cent aurifl If we assume constant M/L = 2, the current total mass in GCs in the 
Harris (19961 ) catalog is ~ 3.5 x 1O 7 M , or ~ 0.07% of the baryonic mass of the Galaxy and 
0.005% of the total mass (including dark matter). The luminosity function of GCs in the 
Galaxy peaks at My rs —7.4 mag, correspondin g to a typical mass of ~ 2 x 10 5 M , and has 
a (Gaussian) width of ~ 1 .2 mag flHarris 200lh . The total initial mass for GCs is estimated 
to be ~ 4 - 8 x 10 s M flFall fc Zhang 200lf ). Apparentl y more than 90% of all glo bular 
star clusters have been disrupted during the last ~ 12 Gyr (jChernoff &: Weinberg 19901 ); the 
inferred total mass in disrupted clusters is comparable to the total mass of the Galactic 



surface to the rms ((v 2 ) 1 / 2 ) velocity in the cluster: 







,2x1/2 



(1) 



c T hese estima tes are made assuming a constant mass-to-light ratio M/L = 2, with data from 
the (lHarris 19 96) catalog of Milky Way GCs http:// www. physics .mcmaster . ca/Globular .html 
Another useful catalog for Milky Way globular cluster data is a vailable online 
http : //www. astro . caltech. edu/^george/glob/data.html (Diorgovski fc Mevlan 1993 ). 
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Figure 1: Left: Distribution of young (< 100 Myr, bullet s) and old (> 3 G yr, circles) open 
clusters in the Galactic plane, based on the catalog of iDias et al. (20021 ). The old open 
clusters are found preferentially towards the Galactic anti-center and above the plane. The 
young massive clusters (squares) seem to be concentrated in the same quadrant as the Sun, 
which is an observatio nal selection effect. Right: Distribution of old globular clusters, data 
from the lHarris (19961 ) catalog. 
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Figure 2: Radius-mass diagram of Milky Way open clusters, young massive clusters and old 
globular clusters . Ope n cluster half-mass radii r^ m (see §1.3.2)) and masses are taken from 
Dias et al. (20021 ) and lLamers et al. (20051 private communication), respectively. Data for 
the YMCs are discussed in more detail in §5J Globular cluster data are taken from the Harris 
catalog. Dashed and dotted lines represent constant half-mass density ph = 3M/87rr hm 3 and 
half- mass relaxation time t r h (Eq. U7|) . respectively. 



stellar halo (IHut fc Djorgovski 19921 . iBell et al. 20081 see also §U), although w e note that the 
spatial distributions of halo stars and GCs differ (jForte. Vega fc Faifer 20091). 



The open cluster databases of iKharchenko et al. (20051 ) and lPiskunov et al. (20081 ). with 
81 clusters, are probably complete to a distance of ~ 600 pc and h ave a mean cluster mass of 



500 M^. With an o pen cluster birth rate of 0.2-0.5 Myr" 1 kpc" 2 flBattinelli fc Capuzzo-Dolcetta 1991 



Piskunov et al. 20061 ). the implied total star formation rate in open clusters is ~ 2x 10 2 M Myr" 1 kpc" 2 
For an average cluster age of ~ 250 Myr these estimates imply a total of about 23,000-37,000 
open star clusters currently in the Galaxy. However, the formation rate of embedded clus- 
ters (still partly or completely enshrouded in the molecular clo ud from which they formed) 
is considerably higher — 2 — 4 Myr -1 kpc - ( Lada fc Lada 20031 ) — and with a similar cluster 
mean mass the total star formation rate in embedded clusters is ~ 3 x 10 3 M Myr -1 kpc -2 , 
comparable to the f ormation rate of field stars in the disk (3 — 7 x 10 3 M Q Myr -1 kpc -2 ; 
Miller fc Scalo 19791 ). Although the uncertainties are large, this suggests that the majority 
of stars form in embedded clusters, but only a relatively small fraction (~ 10%) of clusters 
survive the embedded phase. 

These estimates are sensitive to the underlying assumptions made about the star-formation 
history of the Galaxy and the duration of the embedded ph ase, as well as to observational se- 
lection effects. For example, the open cluster sample used bv lBattinelli Capuzzo-Dolcetta (19911 ) 
is based on a luminosity limited sample of 100 clusters from the lLynga (19821 ) catalog, which 
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Figure 3: Top view of the Milky Way Galaxy, with the known spiral pattern ( jVallee 20081 ) 
and young star clusters more massive than > 10 4 M identified. Basic cluster parameters 
are listed in Tab. El The location of the Sun is indicated by 0, and its orbit by the dotted 
circle. Dashed lines indicate circles of 1 kpc and 2 kpc around the sun. 



claims to be complete to a distance of 2 kpc, but the mass of open clusters in this catalog be- 
tween 600 pc and 2 kpc averages several thousand solar masses. The much higher mean mass 
of open clusters at large distance indicates that care has to be taken in using these catalogs, 
as there appear to be selection effects with respect to distance. Another problem arises from 
confining the analysis to a distance of 600 pc around the Sun, since th e cluster sample does 



not include any nearby spiral arms, where many young clusters form; lLada fc Lada (20031 ) 
considered a sample of clusters within 2 kpc of the sun, which therefore includes many objects 
in the Perseus and Sagittarius arms (Fig. [3j) . These differences complicate direct comparison 
of cluster samples. 

1.3 Terminology 

The study of star clusters has suffered from conflicting terminology used by theorists and 
observers. In this section we attempt to clarify some terms, with the goal of making this 
review more readable by all. 
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1.3.1 Cluster center 



Determining the center of a star cluster sounds like a trivial exercise, but in practice it is not 
easy. The cluster center is not well defined observationally, although theorists have reached 
some consensus about its definition. 



von Hoerner (19631 ) defined the center of a simulated (iV-body) cluster as a density- 
weighted average of stellar positions: 

Y, ■>:/>' J 



X d,j — ^ j (2) 



where pf^ 1 is the density estimator of order j around star i, and Xi is the (3- dimensional) 
position vector of star i. 

In direct iV-body simulations (see §A.3I) . alternatives to Eq.[2] are preferred due to the 
computational expense of determining the local density p\ . The center of mass, often 
used in simple estimates, is generally not a good measure of the cluster center, as distant 
stars tend to dominate. This has led to approximate, but m ore efficient, estimators, such 



as the "modified" density center (IPortegies Zwart et al. 20011 ) . which iteratively determines 



the weighted mean of the positions of a spe cified Lagrangian fraction (typicall y ~ 90%) of 



stars, relative to the modified density center ( iHeggie. Inagaki fc McMillan 1994T ). In general, 



it agrees well with the density center defined above. 

Observationally, the cluster center is considerably harder to define, both because of the 
lack of full 3-dimensional stellar positions, and also because of observational selection effects, 
including the influence of low-luminosity stars and remnants, crowding, and the broad range 
in luminosities of individual stars. Both number-averaged and luminosity-averaged estima- 
tors are found in the literature. In principle, the 2-dimensional equivalent of Eq.|2] could be 
used, but observers often prefer the point of maximal symmet ry of the observed projec ted 



stellar distribution. An example is the technique adopted by iMcLaughlin et al. (20061 ) to 
determine the center of GC 47 Tuc. 



1.3.2 Size scales 

Massive star clusters tend to be approximately spherically symmetric in space, or at least 
circular on the sky, so the radius of a cluster is a meaningful measure of its size. Theo- 
rists often talk in terms of Lagrangian radii — distances from the center containing specific 
fractions of the total cluster mass. For observers, a similar definition can be formulated in 
terms of isophotes containing given fractions of the total luminosity The half-mass radius 
(7h m ; the 50% Lagrangian radius) is the distance from the cluster center containing half of 
the total mass. Observationally, the projected half-light radius — the effective radius r eff — is 
often used, although the total cluster light, which is obviously required to define the La- 
grangian radii, can be hard to determine and is not the same as rh m in projection when the 
mass-to-light ratio varies with the distance to the cluster centre. 

Arguably more useful cluster scales are the virial radius r vir , the core radius r c (see Eq.[6]), 
and the tidal radius r t . We now define them in turn. 
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The virial radius is defined as 

rvir " 2\U\ ■ [ ' 

Here M is the total cluster mass, U is the total potential energy and G is the gravitational 
constant. This is clearly a theoretical definition, as neither the total mass nor the potential 
energy are actually observed. The potential energy may be obtained directly from the stellar 
masses and positions in an iV-body simulation (see §4.ip . or from a potential-density pair 
by U = 2tt J p(r)<f>(r) r 2 dr. 

From an observational point of view, the parameter 77 = 6r vir /r e g is generally introduced 
to determine the dynamical mass of star clusters. In virial equilibrium (U = — 2T, where T is 
the total kinetic energy of the cluster stars), T/M = \{v 2 ) = §o"id 2 f° r an isotropic system. 
The line-of-sight velocity dispersion ctid can be directly measured, yielding the cluster mass 

(4) 

In Fig. H] we present the dependence of rj on the parameters of some typical density profiles: 



the concentration parameter c = log(r t /r c ) of a lKing (19661 1 model or the parameter 7 in an 



Elson. Fall &: Freeman (19871 hereafter EFF87) surface brightness profile, 

2 \ -7/2 
2 



S(r) = E ( 1 + r - ) , (5) 



where a is a scale parameter, the 3-dimensional density profile has a logarithmic slope of 
— 73D = —(7 + 1) for r ^> a, and r c and r t are, respectively, the core radius and tidal radius 
of the cluster (to be d efined below). 



A iPlummer (19111 ) density profile has 7 = 4, r v ; r /r e g = 16/3-7T, and therefore 77 ~ 10. 
The value 77 = 9.75, corresponding to r vir = 1.625r e ff, is a reasonable and widely used choice 
for clusters with relatively shallow density profiles — 7 > 4 or c < 1.8. For 7 < 2 the 
EFF87 profile has infinite mass, and the ratio r w - ir /r e ^ drops sharply for 7 < 2.5. The choice 
for 7] = 9.75 should be made cautiously, since many young clusters tend to have relatively 
shallow density profiles with 2 < 7 < 3, for which rj £ 9 (see Fig.H]and also §5]). In addition, 
mass s egregation can ha ve a severe effect on rj, resulting in a variation of more than a factor 



of ~3 flFleck et al. 20051 ) . 

Observers generally define the cluster core radius, r c , as the distance from the cluster 
center at which the surface brightness drops by a factor of two from the central value. 
Unfortunately, theorists use at least two distinct definitions of r c , depending on context. 
When the central density p an d velocity dispersion (v 2 ) are easily and stably defined, as 
is often the case for analytic, Fokker-Planck, and Monte Carlo models (see §4.11 and the 
Appendix IA1 for more details), the definition 



( iKing 19661 ) is often adopted. For typical cluster models this corresponds roughly to the 
radius at which the three-dimensional stellar density drops by a factor of 3, and the surface 
density by ~ 2. 
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Figure 4: The ratio r v i r /r e g and the parameter rj used to convert an observed 1-D velocity 
disp ersion and h alf-light radius into a dynamical mass (Eq.HJ) for Eq.[5](left) and lKing (1962l ) 
and King (19661) models (ri ght). The dashed line in the left panel indicates the analytical 
result for a lPlummer (19111 ) model (7 = 4 in Eq. [5]). 



In TV-body simulations, however, both po and (v 2 )o are difficult to determine, as they are 
subject to substantial stochastic fluctuations. As a result, a density-weighted core radius is 
used instea d. Specifically, for each star a local density pi is defined using the star's k nearest 
neighbors (jCasertano fc Hut 19851 ). where k = 12 is a common choice. A density center is 
then determined, either simply the location of the star having the greatest neighbor density, 
or as a mean stellar position, as in Eq.[2j except that the density estimator is pf. [The sq uare 
is used rather than the first power, as originally suggested by ICasertano fc Hut (19850 . to 
stabilize the algorithm and make it less sensitive to outliers.] The core radius then is the 
p^-weighted rms stellar distance from the density center: 



£ 



% III 



£i p\ 



(7) 



Despite their rather different definitions, in practice the two "theoretical" core radii (Eqs. [6] 
and [7j) behave quite comparably in simulations. 

For simple models, the values of r c and r vir determine the density profile, which is gen- 
erally assumed to be s pherically sy mmetric. This is the case for the empirical King (19621 ) 
profiles and dynamical iKing (19661 ) models, both of which fit the observed surface brightness 
distribution of many Milky Way GCs. The dynamical King models are often parameterized 
by a quantity Wq representing the dimensionless depth of the cluster potential well. Centrally 
concentrated clusters have W £ 8 (c £ 1.8), whereas shallow models have W ^,4 (c^ 0.8). 

1.5). 



The empirical and dynamical King profiles are in good agreement for W$ £ 7 (c 



< 
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c = log 10 (r t /r c ) 

Figure 5: EFF87 model fits (with power-law index 7) to King surface brightness profiles 
(characterized by King concentration parameter c = log 10 r t /r c ). The fits of Eq.[5] to the 
King profiles/models are done within half the cluster tidal radius. 



Galactic GCs are well fit by King models, but the obse rved surface brightness pr ofiles of 
young clusters in, for example, the LMC are not (see e.g. iMackev fc Gilmore 2003h . They 
are much better represented by EFF87 profiles (Eq. [5]), which have cores (different from 
the King-model cores) and power-law halos. For a King model with concentration c > 1 
(W > 5), the surface brightness drops to approximately half of its central value at r = r c , 
as defined by Eq.0 so the observed core radius is a good measure of the core radius of 
the underlying three-dimensional stellar density distribution. The "King" core radius of the 
EFF87 surface brightness profile often adopted by observers is 



( 2 2/7 _ !) 



1/2 



Here a is the scale parameter in the EFF87 profile. Thus, when an EFF87 surface brightness 
profile is fit to an observed cluster, Eq.[S]can be used to determine with good confidence the 
3-dimensional core radius r c of Eq.[6j 

In Fig. [5] we compare the EFF87 profiles with the empirical King profiles and (pro- 
ject ed) King mo dels, by fitting Eq.|5]to each within the inner half-mass radius. For c — > 00 
the King (1962h surface brightness profile tends to a power-law with index —2, which has 
(logarithmically) infinite mass. The King (19661 ) model in that case becomes an isothermal 
sphere (p oc r~ 2 ), also with (linearly) infinite mass, corresponding to 7 = 1 in Eq.[5J 

The tidal radius, r t , is the distance from the center of a star cluster where the grav- 
i tational accelerati on due to the cluster equals the tidal aceleration of the parent Galaxy 
( jvon Hoerner 19571 ). In Fig.[6]we show the equipotential curves of a cluster in the tidal po- 
tential of its parent galaxy (the galactic center is to the left in the figure). The equipotential 
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Figure 6: Equipotential surfaces for a Wq = 3 star cluster in an external point-mass tidal 
field. The galactic center is to the left. The Jacobi surface is delineated by the two crossing 
equipotentials, passing thro ugh the Li Lagrangian po int to the left and the L 2 point to the 
right. (Data from Fig. 2 of iFukushige fc Heggie 20001 ) . 



surface through the two Lagrangian points (Li and L 2 at x ~ ±3.1 in the figure) is the Jacobi 
surface, defining the effective extent of the cluster in the ex ternal field. Stars within this 
surface may be regarded as cluster members (IRead et al. 20061 ) . Stars do not e scape from th e 
cluster in random directions, but instead do so through the L\ and L 2 points (IHeggie 20011 ) . 
As a practical matter, the Jacobi radius rj, the distance from the center to the L\ point, is 
a commonly used measure of the cluster's "size." For clusters on circular orbits rj is defined 
by the cluster mass, M, th e orbital angular frequency in the galaxy, u, and the galactic 
potential, 0, by ( iKing 19621 ) 



Here, u = Vq/Rg, where Rq is galactocentric distance and Vq is the circular orbital speed 
around the galaxy center, and we have assumed a flat rotation curve: 4>(Rg) = VqI^Rq)- 
Note that a factor 2/3 is sometimes included to correct for the elonga tion in the direction 
along the line connecting L\ and L 2 ( jlnnanen. Harris fc Webbink 19831 ). 

For "Roche lobe filling" clusters that exactly fill their Jac obi surfaces, the Jacobi radius is 
often identified with the truncation radius of a lKing (19661 ) model. However, we emphasize 
that there is no compelling reason why a star cluster should exactly fill its Jacobi surface, nor 
is there necessarily any connection in general between the King truncation radius and the 
properties (or even the existence) of an external tidal field. Observationally, the truncation 
radius is not measured directly, but instead is inferred from King model fits. 
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1.3.3 Time scales 



The two fundamental time scales of a self-gravitating system are the dynamical time scale 
tdyn, and the relaxation time scale t T \. 

The dynamical time scale is the time required for a typical star to cross the system; it is 
also the time scale on which the system (re)establishes dynamical equilibrium. A convenient 
formal definition in terms of conserved quantities is 



t. 



dyn 



GM 5/2 



(10) 



where E = T + U is the total energy of the clu ster. In virial equilibrium, IT + U — and 
this expression assumes the more familiar form ( ISpitzer 19871 ) 
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-1/2 
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2 x 10 4 yr 
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10 6 M m 
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1 pc 
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12) 



The relaxation time, £ r i, is typically much longer than t<jyn- It the time scale on which two- 
body encounters transfer energy between individu al stars and c ause the system to establish 
thermal equilibrium. The local relaxation time is ( ISpitzer 19871 ) 



2\3/2 



tr 



(v 2 ) 



15.4G 2 mplnA' 



(13) 



where m is the local mean stellar mass and p is the local density. The value of the parameter 
A is OAN for the theoretical case where all stars have the same mass and are distributed ho- 



mogeneously with an isotropic velocity distribution ( ISpitzer 19871 ). iGiersz fc Heggie (19941 ) 
find empirically A ~ 0.11 A for systems where all stars have the same mass. For systems 
with a significant range of stellar masses, the effective value of A may be considerably smaller 
than this value. 

For a cluster in virial equilibrium we can replace all quantities by their cluster-wide 
averages, writing (v 2 ) = GM/2r vir and p ~ 3M/87rr^ ir , where we ignore the distinction 
between virial and half-m ass quantitie s, so r nm ~ ?Vir- We thus obtain the "half-mass" 
two-body relaxation time (ISpitzer 19871 ) 



^rh 



0.065(v 2 ) 3 / 2 
G 2 (m)plnA 



0.14 



Cl/2( m )V2 lnA 



N 
7hiA 



t, 



dyn > 



(14) 
(15) 
(16) 



where (m) = N/ M is the global mean stellar mass and N is the total number of stars in the 
cluster. If, for simplicity, we adopt In A = 10 as appropriate for the range in cluster masses 
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of interest in this review, Eq. [TH becomes 



Finally, we note that in real stellar systems the one-parameter simplicity of Eq.HH] is 
broken by the introduction of a third time scale independent of the dynamical properties of 
the cluster — the stellar evolution time scale t$ ~ 10 Myr for YMCs, but a complication here 
is that ts depends on time. This simple fact underlies almost all of the material presented 
in this review. 



2 Properties of young massive star clusters 



2.1 General characteristics 

Traditionally, astronomers have drawn a clear distinction between the relatively young and 
low mass Milky Way open clusters associated with the Galactic disk and the old and massive 
globular clusters that reside mostly in the bulge and halo (see Fig. [Hand Tab. [p. According 
to our definition in ^TJ the Milky Way hosts several clusters that fill the gap between these 
populations, in terms of both mass and density (Fig. |2]), indicating that the formation of 
clusters with masses comparable to old GCs is not restricted to the early universe. This 
becomes even more evident when we loo k at the Magellani c clouds, which host ma ssive 
clusters spanning a broad range of ages ( jHunter et al. 20031 Ide Grijs fc Anders 20061). a s 
well as many YMCs that have received considerable attention since the 1960s ( jHodge 19611 ) . 

The ages of the YMCs in the Magellanic Clouds are comparable to those of many Milky 
Way open clusters (up to a few hundred megayears), but the masses and core densities of 
these clusters exce ed those of open clus ters in the Milky Way, in some cases by several orders 
of magnitude (e.g. lElson fc Fall 1985bl see Tabled]). A prominent example is R136, whose 
core was once believed to be a singl e stellar object at least 2000 times more massive and about 
10 8 times brighter than the Sun. IWeigelt fc Baier (19851 ) unambiguously resolved it into a 
group of stars, and now we know it is a cluster of ~ 10 5 young stars, rather than a singl e ex- 
traordinary object (IMassey fc Hunter 19981 . lAndersen et al. 20091 ICampbeU et al. 2010L and 
left panel of Fig. TJ. YMCs such as R13 6 are responsible for the giant HII regions found 
in other galaxies ( iKennicutt fc Chu 19881 and right panel in Fig. [TJ and appear to be much 
more common phenomena than was previously thought. 

YMCs more massive than R 136 were already known several decades ago. Ground-based 
observations revealed numero us "bluish knots" and "super-star clusters" in the starburst 
galaxies M82 and NGC 1569 (Ivan den Bergh 197ll . lArp fc Sandage 1985), and in the rela - 
tively local ongoing galaxy mergers NGC 7252 and NGC 3597 (jSchweizer 19821 iLutz 19911 ). 
HST has since confirm ed many extragalacti c YMCs, starting with those in the interact- 
ing galaxy NGC 1275 ( IHoltzman et al. 19921 ). These and later well studied examples have 
been cited as promising cand idates for the latest generation of "young globular clusters" 
dWhitmore 20031 lLarsen 2006h . 

We review here the basic characteristics of known YMCs and attempt to assess the 
similarities and differences between the young and old populations (see also §5.4l for a number 
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Figure 7: Left: A region of 50 x 50 pc 2 around the ~ 10 5 M cluster R136 in the 30 Doradus 
region of the LMC, at a distance of ~ 50kpc. Right: Many young star clusters forming 
in M83 at a distance of 3.6 Mpc. Credit: NASA, ESA, the Wide Field Camera 3 Science 
Oversight Committee, and the Hubble Heritage Team (STScI/AURA) and F. Paresce R. 
O'Connell (R136) and R. O'Connell (M83). 
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of well studied cases). In the limited space available in this section we highlight only the 



Whitmore (2003b and Larsen (2006. 


2009b) 


in depth in 


de Grijs & Parmentier (2007 


)• 



We have collected data on clusters for which estimates of the age, (photometric) mass 
and half-light (or effective) radius exist in the literature. If more information on structural 
parameters is available, we determine r vir using the relation presented in Fig. HJ If only 
an estimate of r eS is available, we determine r vir using the Plummer value of r vir /r cff = 
16/(37r) ~ 1.6. We then calculate td yn from r vir and the mass using Eq. [12j We assume here 
that the observed clusters are in virial equilibrium, which results in an overestimate for t& yn 
for unbound (expanding) clusters . The data are presented in Tables El [3] and HI 



As illustrated in Figure 2 of iPfalzner (20091 ). young star clusters appear to show two 
evolutionary sequences: (1) a compact configuration starting at a radius of ~ 0.5 pc, which we 
refer to as (dense) clusters, and (2) a sequence with radii ~ 5 pc, which w e call associati ons. 



(These two groups are referred to as "starburst" and "leaky" clusters by Pfalzner 20091 .) A 
more quantitative distinction may be found in the ratio age/idyn, since for bound objects this 
ratio should be large and fo r unbound obje cts it is expected to remain small. We note that 



the two groups identified by IPfalzner (20091 ) can be divided in this way, with a boundary at 
age/tdyn ~ 3. This division is indicated in the tables. More physical interpretations of these 
sequences are presented in §H1 



2.2 Mass segregation 

Nearby clusters that can be resolved into individual stars are excellent laboratories for studies 
of stellar populations, and their derived properties can be used as input parameters to models 
of cluster evolution. Two well known exa mples of resolved YMCs are the Arches and Quin - 
tuplet clusters near the Galactic center (iFiger et al. 20021 . iFiger. McLean fe Morris 19991 ). 



The Arches has an age of just 1-2 Myr and a central density of ~ 10 5 M M pc 3 ; its stel- 



lar initial mass function (I MF) has been the topic of much debate ( jStolte et al. 20021 12005 



Kim et al. 20061 ). (See also lHarfst. Portegies Zwart fc Stolte 20091 who reconstruct the clus- 
ter's initial conditions by iterated iV-body simulations.) We will refrain from commenting in 
detail on this debate, and instead refer to the recent review bv lBastian. Covev fc Mever (2010h 
for an in-depth discussion of this topic. The current consensus seems to be t hat "evidence fo r 
IMF variations is absent, although this is not evidence for their absence." A Salpeter (1955ll R 
dist ribution cannot be ruled out for the stellar IMF of the Arches cluster (IKim et al. 20061 ). 



Stolte et al. (20021 ) find evidence for mass segregation within the Arches cluster based 
on a steepening o f the stellar mass function (MF) wit h increasing distance from the cluster 
center. However, lEspinoza. Selman &: Melnick (20091 ) suggest that the apparent mass seg- 
regation could simply be an observational bias, and that a Salpeter MF in the core is still 
consistent with their observations. The differential extinction across the cluster and the total 
visual extinction of Ay m 30 magnitudes make this a challenging cluster to study. See also 



d Here we use the term "Salpeter" as a possibly universal power-law mass function (with exponent -2.35) 
for stars more massive than ~ 1 M©. At present, given the range of detectable masses, it is not possible to 
distinguish more finely between the various mass functions commonly adopted by theorists. 
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(2) 


(3) 
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0.06 
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3.50 
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1.20 
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0.55 


6.41 


NGC 3603 


5 
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0.15 


0.70 


2.00 
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0.17 


11.62 
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6 
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2.00 






3.40 


0.93 


4.29 
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6 
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4.50 


4.70 




1.50 




3.70 


2.55 


0.34 


35.22 


RbGC 02 


8 


17.00 


4.60 


4.80 




2.70 




3.40 


4.58 


0.73 


23.18 


RbGC 03 


9 


18.00 


4.50 






5.00 






8.49 


2.07 


8.68 


Trumpler 14 


10 


2.00 


4.00 




0.14 


0.50 


2.00 




0.85 


0.12 


17.15 


Wd 1 


11 


3.50 


4.50 


4.80 


0.40 


1.00 


4.00 


5.80 


1.74 


0.19 


18.27 


Wd 2 


4 


2.00 


4.00 






0.80 






1.36 


0.24 


8.48 


hPer 


4 


12.80 


4.20 
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3.57 


0.80 
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4 
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4.10 






2.50 






4.24 


1.16 
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CYgOB 


4 


2.50 


4.40 
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11.58 


0.17 


I Lac 1 


4 
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3.40 
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35.14 


61.97 


0.23 


Lower Cen-Crux 


4 


11.50 


3.30 
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25.46 


42.89 


0.27 


NGC 2244 


4 
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4.90 


0.41 
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4.40 
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NGC 7380 


4 
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0.29 
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0.72 
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1.70 


3.60 
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8.26 


0.21 
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4 
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0.25 
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3.60 
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54.30 


0.27 
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4 
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0.18 



Table 2 

distinction based o n age/t 



Stolte et al. (2002): 4: 


6: 


Fieer et al. (20061): 


Clark et al. (2009 


); 10: 


12: 


Brandner et al. (2008) 



*™. 1: Fierer. McLean & Morris (1999) 


2: 


Fierer et al. (2002h: 3 


Pfalzner (20091): 5: 


Haravama. Eisenhauer & Martins (2008^ 


7: 


Davies et al. (2008): 8: 


Davies et al. (2007): 9 


Ascenso et al. (2007): 


11: Meneel & Tacconi-Garman (2007) 


13: Hillenbrand & Hartmann (1998). 
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Table 3: Same as Table [21 but now for the Local Group. 1: 


Hunter et al. (1995) 


2: Mackev & Gilmore (2003 


): 3: McLauehlin & van der Marel (20051): 4 


Andersen et al. (2009): 5: 


Vansevicius et al. (2009): 6: 


Barmbv et al. (2009) 


7: Perina et al. (2009): 8: Wvder. Hodsre & Zucker (2000): 9 


Ghandar. Bianchi & Ford (2000h: 10: Mafz-Apellaniz (2001): 11: Sabbi et al. (2008). 
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Ascenso. Alves fc Lago (20091 ) for a general discussion on observational challenges in studies 



on mass segregation in the dense central regions of clusters. 

Currently the most massive young cluster known in our Galaxy is Westerlund 1 (IClark et al. 20051 ). 
Due to its relative proximity (~4kpc), its lower extinction (although still a considerable 
Ay « 10 magnitudes), and slightly lower intrinsic stellar density, it represents a somewhat 
easier target than the Arches and Quintuplet for studies of its stellar content. The MF and 
structural parameters of Westerlund 1 and the somewhat less massive NGC 3603 have been 
investigated by several space-based instruments and from the ground using adaptive optics. 
Evidence for an overa bundance of massive stars within the ha l f-mass radius has bee n found 
in both these clusters (IHarayama. Eisenhauer fc Martins 20081 ; iBrandner et al. 20081 ) . Prob- 
ably the best kn own example of a mass segregated cluster is the (much closer) Orion Nebula 
Cluster (ONC) (jHillenbrand fc Hartmann 19981 . iHuff fc Stahler 20061 ). Radial variations of 
the st ellar MF have also been found for slightly ol der (~ 10 — 25Myr) star clusters in the 
LMC fide Griis et al. 2002al . iGouliermis et al. 20041 ). 

Clusters at much larger distan ces have also been used for studies of mass segregation. 



McCrady. Graham fc Vacca (20051 ) show that cluster F in the starburst galaxy M82 appears 



to be smaller at red wavelengths than at blue wavelengths, which they attribute to mass seg- 
regation. The reasoning is that the most massive stars are red (super)giants and if these are 
more cen trally concentrated th e cluster would appear smaller in the red. A note of caution is 
added by lBastian et al. (20071 ). however, who point out that M82-F is a complicated case be- 
cause of differential extinction across the clus ter. A negative cor re lation between radius and 
wavelength was also found for NGC 1569-B f lLarsen et al. 20081 ). lOaburov fc Gieles (20081 ) 
modeled mass segregated clusters and projected them in different filters, and showed that 
mass segregation can only explain differences in measured r c in different filters at the 5% 
level. This because the massive stars are much more luminous than the average cluster 
member and dominate the light profile at all wavelengths, so if a cluster is mass segregated 
it should appear smaller at all wavelengths. It is thus not clear why NGC 1569-B and M82-F 
should appear smaller at redder wavelengths. 

If stars with masses greater than some mass m are found to be more centrally con- 
centrated than the average stellar mass (m), and if the cluster is much younger than the 
dynamical friction time scale for stars of mass m (~ (m)t r h/m; see Eq. [20] in §|3]), then an 
obvious conclusion is that the observed mass segregation is primordial — that is, the cluster 
formed with stars more massive than M preferentiall y closer to th e center. This is gener - 
ally consistent with simulations of cluster formation (IKlessen 200 ll iBonnell fc Bate 20061 ). 
However, the time scale argument only holds for a spherically symmetric stellar system in 
virial equilibrium. If clusters form through mergers of smaller sub-clumps, the degree of 
mass segregation of the merger product is higher than would be expec t ed from these sim- 



ple dynamical arguments ([McMillan. Vesperini fc Portegies Zwart 20071 ; lAUison et al. 2009 



Moeckel &: Bonnell 20091 ). From studies of star clusters that are still embedde d in their natal 



molecular clouds, it is e vident that stars form in a clumpy hierarchical fashion (ILada fc Lada 20031 . 



Gutermuth et al. 20051 ). and that some merging has to occur during the early (embedded) 



evolution. A YMC with a clear lack of mass segregation would therefore be an very inter- 
esting object to study. 

The degree of mass segregation at a young age has important consequences for the further 
evolution of a cluster, as we discuss in more detail in 331 and 
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2.3 Structural Parameters 



Studies of the structural parameters of YMCs often involve the rich cluster population of 
the LMC. Due to their vicinity, these clusters can be spatially resolved and studied from the 
ground. As mentioned in §1.3.21 unlike the old GCs, whose surface brightness profiles are 
best described by (tidally truncated) King profiles, the surface brightness profiles of YMCs 
are best described by power- law profiles with a core (Eq. EFF87). A typical range for 
the power-law index is 2.2 < 7 < 3.2 (EFF87; iMackey &: Gilmore 20031 ). Similar slopes 



are found for extragalactic clusters, 
dLarsen 20041 ) . 



with tentative evidence for an increase of 7 with age 



Several studies have reported a striking incr e ase of r r with age for L MC clusters (lElson. Freeman fc Laue 



de Grijs et al. 2002bl : IMackey fc Gilmore 20031 1 . iBastian et al. (20081 ) finds a similar increase 



for r c of massive clusters in M51 and in a compilation of literature data. iBrandner et al. (20081 ) 
have shown that r e g of YMCs in the Milky Way increases with age. Fig. [8] shows the in- 
crease of r c and r e fj with age for all the clusters in Tables EHU There may be a rather strong 
selection bias toward dense objects in the studies that consider r c , i.e., there may be young 
clusters with large r c which simply have not been classified as star clusters or do not have a 
r c measurement available. This becomes clear immediately when we look at r e g, where we 
have many more values available for the associations (open symbols). The values of r c are 
not determined for these objects, resulting in a deplet ion of points at the top left of F ig. |HJ 



At an age of 100 Myr there is a large spread in r c . lElson. Freeman fc Lauer (19891 ) show 
that mass loss due to stellar evol ution from stellar pop ulations with different IMF slopes gives 
rise to different expansion rates. IMackey et al. (20081 ) invoke different degrees of (primordial) 
mass segregation and retention fractions of black holes to explain a range of growth rates for 
r c (see §3.4.3l for further discussion). From the lower envelope of points in the right panel of 
Fig. [8] it seems that r e g rises monotonically with age, with no obvious increase in the spread 
with age. 

An increase of a factor of 5-10 in radius has dramatic implications for the evolution 
of these clusters, since it implies that in a very short time the densities drop by two or 
three orders of magnitudes. However, if very young clusters are mass segregated ( §2.2)) . the 
observed r c and r e s (measured from the p rojected luminosi t y profile, see ^1.3.21) cou ld be 
considerably smaller than the true values (IFleck et al. 20061 iGaburov fc Gieles 20081 ). We 
will return to this topic, and the various physical mechanisms that drive cluster expansion 
at these ages, in more detail in §3] and §HJ 

Several studies have discussed the lack of any clear correlation between the size of a 
cluster and its mass or luminosity ( IZepf et al. 19991 lLarsen 20041 . IScheepmaker et al. 20071 ). 
If this is how clusters form, it places important constraints on models of cluster formation, 
since m olecular cloud s and dense cores do follow a (virial) scaling relation between mass and 
radius (ILarson 198ll ). A variable star formation efficiency that i s higher for massive clu mps 
could weaken or erase the molecular cloud mass-radius relation (lAshman &: Zepf 20011 ) . 

Fig. [9] shows the masses and radii of all clusters and associations listed in Tables [2H4J 
The sample is divided into two age bins: less than 10 Myr (left panel) and 10 — 100 Myr 
(right panel). We tentatively overplot lines of constant density, which might indicate a trend 
with p hm Rj 10 3±1 M for the younger cluster sample. This mass-radius correlation seems to 
contradict earlier findings, but it is important to bear in mind that we have limited ourselves 
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Figure 8: Core radius (r c , left) and effective radius (r e e, right) as functions of age for 
the clusters in Tables [2} |3] & HI The filled symbols are clusters and the open symbols are 
associations. 



to a narrow range of (young) ages, and show the data for the clusters (age> 3td yn ) with 
different (and more distinct) symbols, whereas most literature studies do not separate the 
data in this way. The distinction between clusters and assocations is of course a density 
cut (Eq. [I2"|) . corresponding to phm ~ 10M Q pc~ 3 at an age of 10 Myr, causing the upper 
envelope of cluster points (filled) to line up with the constant density lines. However, the 
lower envelope of points in the left panel seems more consistent with constant density than 
with constant radius. For the older clusters (right panel) the lower envelope seems quite 
consistent with constant radius. A correlation between mass and radius, or indeed the lack 
of one, would have important implications for the cluster's long-term survival, as we discuss 
in more detail in §H 



2.4 Global properties 



The YMCs in the Magellanic Clouds and other nearby galaxies, such as M31 ( 


Narbutis et al. 2008 




Barmbv et al. 2009. 


Ma et al. 2009 


Vansevicius et al. 2009. 


Hodee et al. 2009. ( 


Caldwell et al. 2009. 


Peacock et al. 2010h and M33 (e.s. < 


3handar et al. 1999; 


Park. Park & Lee 2009 


Saraiedini & Mancone 20( 


San Roman et al. 2009 


) provide opportunities for population studies. Observational studies 



of these systems suffer less from the problematic distance and extinction effects that chal- 
lenge studies of Milky Way star clusters. They provide a more global view of cluster pop 



ulations, and are better targets for studies of cluster age distributions (lElson fc Fall 1985b 



Hodge 19871 . iGirardi et al. 19951. Ide Griis fc Anders 20061 . iParmentier fc de Grijs 20081 1 and 
luminosity functions (jElson fc Fall 1985al ). 

A comparison between Table [2] and Table [3] reveals the striking absence of Milky Way 
YMCs with ages between 10 and 100 Myr, whereas in the SMC and LMC, all YMCs ex- 
cept R136 have ages spanning this range. This may well be an observational effect: Due 
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Figure 9: Mass-radius relation for all clusters (filled symbols) and associations (open sym- 
bols) from Tables El [3] & HJ using r hm = (4/3)r eff . Lines of constant half-mass density, 
p hm = 3M/(87rr hm 3 ), are overplotted. The values of the densities shown correspond to 
0.07Myr < t dyn < 0.7Myr (Eq. [TO]). The clusters are subdivided into two groups: younger 
than 10 Myr (left) and older than 10 Myr (right). 



to extinction and foreground and background stars in the Milky Way, it is hard to discover 
clusters without nebular emission and the luminous evolved stars found in very young sys- 
tems. There might be a population of slightly older (10-100 Myr) YMCs in the Milky Way 
still awaiting discovery. 



2.4.1 The cluster luminosity function 



YMCs have been observed and identified well beyond the Local Group, providing excit- 
ing new opportunities for studies of star formation and population synthesis. They are 
found in abundance in galaxies with high star formation rates, such as merging an d inter- 
acting galaxies (e.g. Holtzman et al. 1992 . Miller et al. 1997 . Whitmore et al. 19991). How- 
ever, they are also found in quiesc ent spirals (e.g. Larsen fc Richtler 2000l ; Larsen 2004 ; 
Cantiello. Brocato fc Blakeslee 20091 ). and there are many similarities between the young 
cluster populations in these different environments. For example, the luminosity function 
(LF), defined as the number of clusters per unit luminosity (dN/dL) is well described by a 
power-law dN/ dL ~ L~ a with (a ps 2) (e.g. Iwhitmore fc Schweizer 19951 . iMiller et al. 19971. 



Larsen 2002 . Ide Griis et al. 2003T) . and tends to be slightly steeper at the bright end (IWhitmore et al. 1999 
Larsen 2002L iGieles et al. 2006ah . 

An appealing property of an a = 2 power-law is the fact that the luminosi ty of the most 



luminous object, L m3iX , increases linearly with the total number of clusters, N d . IWhitmore (20031 ) 



showed that the value of L mSLX for clusters in different galaxies scales as L mBX oc N cl \ with 
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7] 0.75 (see also lLarsen 20021 ) . The index 77 is a proxy for the shap e of the bright end of 
the LF, since for a pure power-law r\ = l/(a — 1) ( Hunter et al. 20031 ). A value of 77 = 0.75 
corresponds to a = 2.4, supporting the finding that the bright end of the LF is steeper than 
L~ 2 , since the L max method traces the brightest clusters. 

If a universal correlation between L max and N c \ exists, it means that L max is the result 
of the size of the sample. It is therefore determined by statistics, and no special physical 
conditions are needed to form bright er clusters. A similar scaling betwe en L rm and the sta r 



formation rate (SFR) was found by IWeidner. Kroupa &: Larsen f2004f) and lBastian (20081 ). 
This has been interpreted in various ways. iBastian (2008) uses the similarity between the 
L max — N c \ and L max — SFR relations to conclude t hat the cluster formation efficienc y is 
roughly constant (~ 10%) over a large range of SFR. IWeidner. Kroupa fc Larsen (20041 ) in- 
terpret it as an increase of the most massive cluster mass, M max , with SFR and they conclude 
that M max is set by the SFR of the host galaxy, much like the Weidner fc Kroupa f2006j) 
relation between the most massive star and the mass of the parent cluster. lLarsen (2009al ) 
showed that the L max clusters have a large range of ages, with the brightest ones on average 
being younger than the fainter ones. This implies a dependence of M/L on L max , such that 
the mass of the most massive cluster (M max ) increases more slowly with N d (and the SFR) 
than L max , and hence that M max does not follow the size of sample prediction of the pure 
power-law LF. 



2.4.2 The cluster initial mass function 

It is tempting to interpret the LF as the underlying cluster mass function. As discussed 
above, it is not trivial to gain information about the mass function from the LF, since the LF 
consists of clusters with different ages, and clusters fade rapidly during their first ~ 1 Gyr due 
to stellar evolution. Determinations of cluster initial mass functions (CIMFs) are rare, since it 
is hard to acquire the cluster ages needed to select the youngest and to convert luminosity into 
mas s. Several CIMF determinations have also found power-law funct ions with indices close to 
-2 (IZhang fc Fall 1999L iMcCrady fc Graham 2007L Elk et al. 20031 ). an d other studies have 



found evidence for a truncation of this power-law at the high-mass end ( IGieles et al. 2006b 



Bastian 20081 lLarsen 2009al IVansevicius et al. 20091 ) . T he functional form of the initial mass 
function for young star clusters is well represented by a ISchechter (19761 ) distribution 

<j>(M) = = AM' 13 exp(-M/M*). (18) 
dm 

Here /3 ~ 2 and the Schechter mass M* is equivalent to the more familiar for the galaxy lu- 
minosity func tion. For Milky- Way type spiral galaxies M* « 2 x 10 5 M (Gi eles et al. 2006b 



jjl %J \J rj A. X vj -T- ^ 

Larsen 2009al ). For interacting galaxies and luminous infrared galaxies, IBastian (20081 ) ob 
tains M* > 1O 6 M . 

The upper panel of Fig. [TU] compares cluster mass functions for several galaxies with the 
Schechter function, Eq. [TSJ The lower panel compares the corresponding logarithmic slopes 
of the data with that of the Schechter function. Clearly the mass function of the Antennae 
clusters extends to higher masses, and this is not only due to the larger number of clusters, 
since the slope is also flatter at large masses. Thus the value of M* seems to depend on 
the local galactic environment — it is possible to form more massive clusters in the Antennae 
galaxies than in more quiescent environments. 
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Figure 10: Top: Comparison of mass func tions of cluster s younger than ~ 1 Gyr in dif- 
ferent galaxies. The results are t aken from Larsen f2009ah fL MC. cluster rich spirals and 
cluster poor spirals'): iGieles (20091 ) (M51); IZhang fc Fall (19991 ) (the Antennae galaxies) and 
Vansevicius et al. f2009h (two versions of the M31 cluster mass function). The cluster mass 
functions in spirals are compared to a Schechter function (Eq. [T8|) with M* = 2.5 x 10 5 M 
(dashed curves); for the Antennae, M* = 1O 6 M (dotted curve) is used. Bottom: the cor- 
responding logarithmic slopes of the mass functions. The dotted and dashed curves are the 
logarithmic slopes of the functions shown in the top panel. 
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Indirect indications for a truncation of the cluster mass function also come from statistical 
arguments. If we temporarily ignore the (exponential) truncation, i.e. we adopt a power 
law with j3 = 2 without the exponential factor in Eq.[TBJ we can relate the cluster formation 
rate to the masses of the most massive clusters obs erved. For an overall star formation rate 
of 5000 M Myr" 1 kpc" 2 in the solar neighborhood flMiller fc Scalo 1979h. comparable to the 



average values found in external Milky- Way-type spirals (IKennicutt 19981 ). and assuming 
that ~ 10% of this mass ends up in bound star clusters, we find that the total mass formed 
in lOMyr in clusters within 4kpc of the Sun is ~ 2 x 1O 5 M . For our assumed power-law 
mass function, the most massive cluster contains ~ 10% of the total mass, so the mass of 
the most massive cluster is a few xlO 4 M . Within a 4kpc circle we find Westerlund 1, with 
a mass of ~ 6 x 10 4 M Q (see Tab.|2]), in reasonable agreement with expectations. 

Assuming the same star formation rate out to a distance of ~ 8 kpc from the sun, which 
is a conservative assumption since the star formation rate toward the Galactic center is 
probably higher than in the solar neighborhood, the expected most massive cluster becomes 
a factor of 8 2 /4 2 higher, or about 1O 5 M . Over a time span of 1 Gyr, clusters with masses 
of ~ 10 7 M should have formed within that same distance, but if such a cluster existed it 
would most likely already have been discovered, unless it has been disrupted, which seems 
unlikely. This implies that for a quiescent environment like the Milky Way Galaxy, the 
truncation of the cluster mass function must occur at considerably lower mass than in a 
starburst environment. Similar arguments hold for external galaxies, where the entire disk 
can be seen, and it is found that the high-mass end of the cluster mass function falls off more 
steeply than a power-law with exponent —2 (see Fig.lTOl). 

2.4.3 Cluster formation efficiency 

The number of globular clusters in a galaxy is often expressed in terms of the specific 
frequency, the number of GCs per unit luminosity of the host. For young clusters, this is 
probably not a very meaningful quantity, since these clusters form with a power-law mass 
function (or Schechter function, as in Eq.[18]), and the luminosity o f the host galaxy depends 



strongly on the age of the field star population. For this reason, lLarsen fc Richtler (20001 ) 
introduce the specific luminosity, T L = 100 L clusters /L galaxy for samples of cluster populations 
in different galaxies. They show that, in the £/-band, Tl increases strongly with the star 
formation rate per unit area (Esfr, see Fig. [11]). This suggests that in galaxies with a higher 
EgpR a larger fraction of the newly formed stars end up in star clusters. The light in this 
blue filter is dominated by (short-lived) hot stars, making Tl a tracer of current star and 
star cluster formation. A compelling as pect of Fig. [TT] is that both axes are independent of 



distance. More recently, iBastian (20081 ) has used the L max -SFR relation to derive a cluster 
formation efficiency of ~ 8%. Both these estimates suffer from (unknown) extinction effects, 
and the best way to derive the cluster formation efficiency would be to compare the fraction 
of the mass that forms in clusters to Esfr for a sample of galaxies. 

The cluster formation efficiency is important for understanding the extent to which YMCs 
can be used as tracers of star formation. An int eresting example is t he lack of clusters with 



ages between 4 Gyr and 12 Gyr in the LMC (e.g. Ivan den Bergh 1991I ). This "age gap" could 
be the result of cluster disruption processes, or due to a global pause in the SFR of the whole 
LMC. From comparison of the SFR history and the age-metallicity relation of LMC field 
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Figure 11: Specific [/-band luminosity (T L ) for various cluster populations, as a function of 
Sspr of the parent galax y. The dashed line shows a linear relation (T^ oc Esfr)- Data from 



Larsen fc Richtler f2000h . 



stars, the latter scenario seems preferred, i.e. the LMC had a low S FR between 4 and 12 Gy r 



ago, and this is reflected in both the field stars and the star clusters ( [Harris fc Zaritsky 20091 ). 
This suggests that YMCs and their age distributions can be powerful tools in determining 
star formation histories in more distant galaxies where individual stars cannot be resolved. 
We return to the interpretation of the age distribution and its consequences for cluster 
lifetimes in SPU 



3 Dynamical processes in star clusters 

Studies of the evolution of a young star cluster split naturally into three phases: (1) the 
first few megayears, during which stars are still forming and the cluster contains significant 
amounts of ambient gas, (2) a subsequent period when the cluster is largely gas-free, but 
stellar mass loss plays an important role in the overall dynamics, and (3) a late stage, during 
which purely stellar dynamical processes dominate the long-term evolution of the cluster. 
An upper limit on the dividing line between phase 1 a nd phase 2 is th e time of the first 



supernovae in the cluster, some 3 Myr after formation ( lEggleton 20061 ). since these expel 
any remaining gas not already ejected by winds and radiation from OB stars. The dividing 
line between phase 2 and phase 3 may be anywhere between 100 Myr and 1 Gyr, depending 
on the initial mass, radius, and density profile of the cluster and the stellar mass function. 

The evolution of the cluster during the first phase is a complex mix of gas dynamics, stellar 
dyn amics, stellar evo l ution, and radiative transfer, and is currently incompletely understood 
(see lElmegreen 20071 . iPrice fc Bate 20091 ). Unfortunately this leaves uncertain many basic 
(and critical) cluster properties, such as the duration and efficiency of the star-formation 
process, and hence the cluster survival probability and the stellar mass function at the 
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beginning of phase 2 (see §1]). 

Together, phases 2 aqnd 3 span the "iV-body simulation" stage familiar to many theo- 
rists. Phase 3 is the domain of traditional dynamical simulations; phase 2 is the "kitchen 
sink" phase, during which stellar interactions, stellar evolution, and large-scale stellar dy- 
namics all contribute (see §4.11 and Appendix A). As discussed in more detail below, the 
processes driving the dynamical evolution during these phases are mostly well known and 
readily modeled, allowing significant inroads to be made into the task of interpreting clus- 
ter observations. However, since the outcome of phase 1 provides the initial conditions for 
phase 2, the proper starting configuration for these simulations remains largely a matter of 
conjecture. Theoretical studies generally consist of throughput experiments, mapping a set 
of assumed initial conditions into the subsequent observable state of the cluster at a later 
age. 

Setting aside the many uncertainties surrounding the early (phase 1) evolution of the 
cluster, in this section we mainly describe the assumed state of the cluster at the start of 
phase 2 and the physical processes driving its subsequent evolution. For better or worse, N- 
body simulations generally assume quite idealized initial conditions (summarized in Tab.|5]), 
with a spherically symmetric, gas-free cluster in virial equilibrium, with all stars already on 
the zero-age main- sequence. 



3.1 Initial conditions 

In the absence of a self-consistent understanding of cluster evolution during phase 1, assump- 
tions must be made abo ut the following key cluster properties before a phase 2 calculation 
can begin ( IKroupa 20081 ). It must be noted that, in almost all cases, the choices are poorly 
constrained by observations. 



Table 5: Commonly adopted initial conditions for particle-based simulations of YMCs. 



cluster 
property 


parameter 


min 


max 


number of stars 


N 
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mass function 
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power law: ( 


A(m) oc m~ 2 - 35 O.1M 
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Q = - 
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0.5 


density distribution 


Plummer, King 






tidal field 


r t /rj 




0.25 


1 


concentration 


W 




1 


16 


binary fraction 


h 







1 


mass ratio 


ip(q) -- 


= 1 





1 


eccentricity 


H(e) = 


= 2e 





1 


orbital period 


r(p) 


oc l/P 


RLOF 


hard 



The stellar mass function <j)(m) = alN/dm is typically ta ken to be a "standard" distri- 
bution derived from studies of the solar neighborhood (e.g. lSalpeter 19551 . Miller fc Scalo 1979 . 



28 



Kroupa 20011 ) although it has been suggested that the mass functions of some YM Cs 
may be deficient in low-mass stars and/or "top-heavy" (ISmith fc Gallagher 20011 ) . in 
the sense that the slope d log N/d log m of the mass function at the high-mass end is 
shallower (i.e. less negative) than the standard Salpeter value of —2.35. 

In addition to the mass function, in many cases minimum and maximum stellar masses 
are imposed. This is necessary for a pure power law, since the total mass would in 
general otherwise diverge, and is often desirable for other distributions, for which 
convergence is not an issue. Often the minimum mass m m i n is chosen to be a relatively 
high, m m i n ^ 1M , to emulate a more massive cluster simply by ignoring the low- mass 
stars. This may suffice if one is interested in clusters younger than 100 Myr (phase 2), 
but for older clusters the lower-mass stellar population is important, e.g. to reproduce 
the proper relaxation time. The maximum mass m max rarely poses a practical problem, 
although there may be some interes ting correlations between the total cluster mass and 
the mass of the most massive star (IWeidner fc Kroupa 20041 ) 



Mass Segregation. Traditionally, dynamical simulations have begun without initial 
mass segregation — that is, the local stellar mass distribution is assumed not to vary 
systematically with location in the cluster. There is no good reason for this, other 
than sim plicity. Evidence for initial mass segregation can b e found in some young clus- 
ters (e.g. Hil lenbrand fc Hartmann 19981. ISabbi et al. 20081 ). simulations of star forma- 



ti on fe.g. iKlessen 20 01. Bo nnell fc Bate 20061). and dynamical evolu tion during phase 



1 ([McMillan. Vesperini fc Portegies Zwart 20071 : lAllison et al. 20091). Several prescrip- 
tions have been used recently for initial mass segregation (ISubr. Kroupa fc Baumgardt 2008 
Baumgardt. De Marchi &: Kroupa 2008 : Vesperini. McMillan &: Portegies Zwart 2009| ). 
They differ in detail, but lead to simila r conclusions, namely that initial mass segrega - 
tion may be critical to cluster survival (jVesperini. McMillan &: Portegies Zwart 20091 ). 
since mass loss from centrally concentrated massive stars can be much more destructive 
than the same mass loss distributed throughout the body of the cluster (see §1]). 

Virial Ratio. Simulations generally begin with a cluster in virial (dynamical) equi- 
librium, with virial ratio Q = —T/U = 1/2. As with most of the other simplifying 
assumptions described here, the principal reason for this choice is reduction of the di- 
mensionality of the initial parameter space, but there is no compelling physical reason 
for it. The gas expulsion that marks the end of phase 1 is expect ed to leave the cluster 
significantly out of equilibrium, and quite possibly unbound (e.g. iHills 19801 ). The time 
scale for a cluster to return to virial equilibrium may be comparable to the time scale 
on which mass loss due to stellar evolution subsequently modifies the cluster structure 
(see gjZp . 

Spatial Density and Velocity Distributions. The initial density profiles of young star 
clusters are p oorly constraine d. Several standard models a re used to m odel the stellar 
distribution: iPlummer (191ll ) and truncated Maxwellian (jKing 19661 ) models are the 
most common . Other distributions, such as isother mal and homogeneous spheres, are 
also used (e.g. IScally. Clarke &: McCaughrean 20051 ). King models provide good fits to 
many observed GCs, although their relation to YMCs is unclear (see §2]). In the absence 
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of strong observational constraints, the stellar velocity distribution is normally taken 
to be non-rotating and isotropic, with (for a Plummer model) the dispersion following 
the local potential with the assumed virial ratio (see §1.3.2p . 

Tidal Field. Clusters do not exist in isolation, but rather are influenced by the local 
tidal field of their parent galaxy In many cases (see §1.3.2j) . the field is modeled explic- 
itly as an external potential or its quadrupole moment relative to the cluster center; 
however simple "stripping radius" prescriptions are also widely used. In practice, in- 
corporating a simple stripping radius instead of a self-consistent tidal field reduces the 
cluster lifetime by about a factor of two if the stripping radius is taken as rj. This 
effect can be mi tigated by adopting a larger c ut-off radius (for exam ple 2rj). With few 
exceptions (e.g. iBaumgardt fc Making 20031 . 1 Wilkinson et al. 20031 ) the parameters of 
the tidal field are held fixed in time, corresponding (for spherical or axisymmetric po- 
tentials) to an orbit at fixed galactocentric radius. For a given orbit and cl uster mass 
the initial ratio of the cluster limiting radius (e.g. the truncation radius of a lKing 1966 



model, see §1.3. 2ft to the Jacobi radius of the cluster in the local tidal field is a free 
parameter, often taken to be of order unity. For a Plummer sphere, which extends 
to infinity, the tidal radius is often implemented as a simple cutoff at some relatively 
large distance. 

Binary Fraction. Binary stars are critical to cluster evolution during phase 2 and espe- 

cially phase 3 (e.g. lHurley. Aarseth fc Shara 20071 ; iPortegies Zwart. McMillan &: Makino 20071 ) 
It is not so clear how important they are during p hase 1, when major structural changes 



are induced by mass loss and mass segregation (jClarke. Bonnell fc Hillenbrand 2000 



Bate. Bonnell fc Bromm 20031 ). There are few, if any, observational constraints on the 
overall binary fraction in YMCs. Open clusters i n the field typical l y have high binary 



fractions, approaching 100% in so me cases (e.g. ISana et al. 2008t iMason et al. 2009 
Bosch. Terlevich fc Terlevich 20091 ) . On the other hand, most recent studies of bi- 
naries in globular clusters suggest binary fractions of between ~ 6% and ~ 1 5 % 



flBellazzini et al. 2002 . ISollima et al. 2007L iMilone et al. 2008L ISommariva et al. 20091 ). 



Binary Secondary Masses. The mass of the secondary star in a binary is typically se 



lected uniformly between some minimum mass and the mass of the primary (IDuquennoy fc Mayor 199 
With this choice, a binary tends to be more massive than the average cluster star, 
resulting in additional mass segregation of the binary population. This effect can be 
removed if desired by randomly selecting primary stars and splitting them into primary 
and secondary components, in which case adding binaries does not affect the mass func- 
tion of cluster members (single stars and binaries), but it does i ntroduce a de viation 
from the initial stellar mass function among binary components (IKroupa 19951 ). 



Orbital Elements of the Binaries. In general, the choices made for binary orbital 
elements tend to be defensive, given the lack of observational guidance. Apart from 
the introduction of a whole new set of initial parameters, the presence of primordial 
binaries also introduces new mass, length, and time scales to the problem, greatly 
complicating direct comparison between runs having different initial conditions. 
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The initial distribution of binary periods is unknown, but is often assumed to fo l- 



low the observed distribution in the solar neighborhood ( jDuquennoy fc Mayor 199ll ). 
which is approximately a Gaussian in log P with mean log P = 4.8 and dispersion 
CHogP = 2.3, with P in days. Ho wever, distributions flat in log P (lAbt 19831 ) and log a 
(IMcMillan. Hut fc Making 19901 ) have also been adopted. Observations of YMCs are of 
little help, as there are at most a handful of binaries with measured orbital parameters 
in YMCs, and those have very short orbital periods and high-mass components. Th e 
eccentricity distribution is usually taken to be thermal ( jDuquennoy &: Mayor 19911 ). 
Binary orbital orientation and initial phase are chosen randomly. 

Higher order multiples. Primordial multiple stars are rarely included in phase-2/3 dy- 
namical simulations. There are a few examples of calculations with primordial triples 
(Ivan den Berk. Portegies Zwart fc McMillan 20071 ) or hierarchical planetary systems 



(ISpurzem et al. 20091 ). The complications of adding primordial multiples greatly in- 



crease the already significant challenges of including binary dynamics and evolution. 



3.2 Multiple stellar populations 



The discovery of multiple p opulations of main-sequence st ars and giants in an increasing 
number of globular clusters ( jPiotto et al. 20051 . iPiotto 20081 ) has led to the realization that 
star clusters are not idealized entities with single well defined stellar populations. In some 
clusters, the observed stellar populations appear to be separated by less than ~ 10 s years, 
well within our age range for young clusters. The existence of multiple populations indicates 
that a second epoch of star formation must have taken place early in the cluster's lifetime 
(see §5.1.1)) . The differences in light-element abundances suggest that the second-generation 
(SG) stars formed out of gas containing matter processed through high-temperature CNO 
cycle reactions in first-generation (FG) stars. 

The two main candidates currently sug gested as possible sources of enriched gas for SG 

formation are rapidly rotating massive stars (IPrantzos &: Charbonnel 20061. iDecressin et al. 20070 

and massive (4-9 M ) Asymptotic Giant Branch (AGB) stars ( Ventura et al. 200 ll . Karakas fc Lattanzio 20 
In order to form the large fraction of SG stars su ggested by observati ons (50% or more of 
the current mass of multiple-population clusters, ICarretta et al. 20081 ). both scenarios re- 
quire that either the IMF of the FG stars was highly anomalous, with an unusually large 
fraction of massive stars, or the FG population had a normal IMF but was in itially at least 



ten times more massive than is now observed. iBaumgardt fc Kroupa (20071 ) have studied 
the subsequent evolu tion and mixing of the two-component cluster in the first scenario. 
D'Ercole et al. (20081 ) have presented simulations of the second, in which the SG stars form 
deep in the potential well of a FG cluster destabilized by early mass loss. Most of the FG 
cluster dissolves, leaving a mixed FG/SG system after a few gigayears. 

If similar processes are operating today, multiple populations should be expected in at 
least some young star clusters, but it is currently not known whether, or to what extent, this 
phenomenon occurs in observed YMCs. For the unresolved extragalactic clusters, multiple 
populations will be hard to confirm, but for clusters in the local group this should be possi- 
ble. At present, however, only one known cluster, Sandage 96, exhibi ts a young (10-16 Myr) 
population together with a relatively old (32-100 Myr) population (jVinko et al. 20091 ) (see 
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also §5.3. ip . Spreads in the main- sequence turn-off have been reported found for interme- 
diate age clusters in the LMC, which has been int erpreted as an age spread of ~ 300 Myr 
( iMackey fc Broby Nielsen 20071 . iMilone et al. 20091 ) . If true, this has exciting implications 
for theories of the formation of star clusters and the general assumption that clusters are 
simple stellar populations. However, iBastian fc de Mink (20091 ) find that the observed main 
sequence turn-off spread could be explained by stellar evolutionary effects induced by ro- 
tation. Obviously if multiple populations are common in YMCs, they significantly impact 
the assumptions made for simulations of phase 2 and the long-term cluster evolution during 
phase 3. Except where noted, we will not explicitly address the possibility of delayed SG 
star formation in this review. 



3.3 Overview of cluster dynamical evolution 

Most numerical studies start with initial conditions as described in ^3.11 To the extent that 



stellar mass loss can be neglected, we can understand the dynamical evolution of a star 
cluster from the fundamental physics of self-gravitating systems, driven by relaxation. 



3.3.1 Evaporation 

The relaxation time (Eq. fToT) is the time scale on which stars tend to establish a Maxwellian 
velocity distribution. A fraction £ e of the stars in the tail of that distribution have velocities 
larger than t> esc and consequently escape. Assuming that this high-velocity tail is refilled 
every t rh , the dissolution time scale is tdis = £ r h/£e- For isolated clusters, v esc = 2 (v 2 ) 1 ^ 2 . 
For a Maxwellian velocity distribution, a fraction £ e = 0.0074 has v > 2 (v 2 ) 1 ^ 2 , and hence 
tdis = 137 t r h- For tidally limited cluster £ e is hig her since iw is lower. For a typical cluster 

density profile £ e ~ 0.033, implying t dis 30t r h (ISpitzer 19871 ). 

The escape fraction £ e is often taken to be constant (e.g. iGnedin fc Ostriker 19971 ). al- 
thouh it depends on rh m (through (t> 2 ) 1//2 ) and also on the stren gth of the tidal field, or r. 



(through Vpxr )- Effectively, £ e depends on the ratio r^/rj (e.g. ISpitzer fc Chevalier 1973 



Wielen 19881 ). iGieles k Baumgardt f2008h show that £ c oc (r hm /rj) 3 / 2 for r hm /rj > 0.05 
(the so-called tidal regime). From Eq. [IB] we then find, for clusters on circular orbits in the 
tidal regime, that t^is oc N/uj (neglecting the slowly varying Coulomb logarithm). For a 
flat rotation curve, tdis o c Re for a cluster of given mass (e.g. IChernoff fc Weinberg 1990 



Vesperini fc Heggie 19971 ) . This linear dependence of tdis on Rq makes it difficult to explain 
the universality of the globula r cluster mass functio n via dynamical evolution of a power-law 
initial cluster mass function (jVesperini et al. 20031 ). but this will not be discussed further 
here; 

Baumgardt (20011 ) showed that td yn also enters into the escape rate and found, for equal- 



mass stars, tdis oc t r h 3 ^ 4 ^d y n 1/4 - The non-linear scaling of the dissolution time with the 
relaxation time results from the fact that a star with sufficient energy to escape may orbit 
the system many times before "finding" on e o f the Lagrangian points, thro ugh which escape 



actually occurs (IFukushige fc Heggie 20001 ) . iBaumgardt fc Making (20031 ) found that this 



scaling also holds for models of clusters with a stellar mass spectrum, stellar evolution, and 
for different types of orbits in a logarithmic potential. Their result for tdis can be summarized 
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as 



t dis ^2Myr 



N 



3/4 



R 



G 



Vn 



kpc/ V220kms" 



-i 



(1 



(19) 



where e is the eccentricity of the orbit. For an eccentric orbit, e > the distance to the galac- 
tic enter Rq is taken as the apogalcticon, whereas Vq is circular velocity, which is constant in a 
logarithmic potential. If the Coulomb logarithm is take n into account, the scaling is approxi- 
mately t& s oc jV°- 65 in the range of about 10 3 to 10 6 M (ILamers. Gieles fc Portegies Zwart 20051 ) 



3.3.2 Core collapse 

Self-gravitating systems are inherently unstable, and no final equilibrium state exists for a 
star cluster. The evaporation of high- velocity stars and the internal effects of two-body relax- 
ation, w hich transfers energy from the inner to the outer regions of the cluster, result in core 
collapse flAntonov 19621 . iLvnden-Bell fc Wood 19681 . ICohn 1980L iLvnden-Bell fc Eggleton 198(1 
Makino 19961 ). During this phase, the central portions of the cluster accelerate toward infi- 
nite density while the outer regions expand. The process is readily understood by recognizing 
that, according to the virial theorem, a self-gravitating system has negative specific heat- 
reducing its energy causes it to heat up. Hence, as relaxation transports energy from the 
(dynamically) warmer central core to the cooler outer regions, the core contracts and heats 
up as it loses energy. The time scale for the process to go to completion (i.e. a core of zero size 
and formally infinite density) is t cc ~ 15t r h for an initial Plummer sphere of identical masses. 
Starting with a more concen trated King (19661 ) distribution shorten s the time of core collaps e 
considerably ( Quinlan 19961 ) . as does a broad spectrum of masses ( Inagaki fc Saslaw 1985 ). 

In systems with a mass spectrum, two-body interactions accelerate the dynamical evo- 
lution by driving the system toward energy equipartition, in which the velocity dispersions 
of stars of different masses would have (mv 2 ) ~ constant. The result is mass segregation, 
where more massive stars slow down and sink toward the center of the cluster on a time 



scale ( ISpitzer 19691 ) 



(m) 



m 



t r i 



(20) 



Portegies Zwart &: McMillan f2002f ) and lGiirkan. Freitag fc Rasio f2004f ) find that, for a typ- 
ical iKroupa (20011 ) mass function, the time scale for the most massive stars to reach the 
center and form a well defined high-density core is ~ 0.2i r i, where t r \ is the relaxation time 
(see Eq. [T3l) of the region of interest containing a significant number of massive stars — the 
core of a massive cluster, or the half-mass radius of a smaller one (in which case t r \ = t r h, 
see Eq. [T6|) . For dense clusters, t s may be shorter than the tim e scale for stellar evolution, or 
for the first supernovae to occur (IPortegies Zwart et al. 19991 ). 

Thus, a collisional stellar system inevitably evolves toward a state in which the most 
massive objects become concentrated in the high-density central core (see §3.4.2|) . Dynam- 
ical evolution provides a natural and effective mechanism for concentrating astrophysically 
interesting objects in regions of high stellar density. 
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3.4 Internal Heating 



On longer time scales, the evolution of clusters that survive the early phases of mass loss 
is driven by the competition between relaxation and a variety of internal heating mech- 
anisms. High central densities lead to interactions among stars and binaries. Many of 
these interactions can act as energy sources to the cluster on larger scales, satisfying the 

relaxation-driven demands of the halo and temporarily stabilizing the core against collapse 

dGoodman fc Hut 19891 : iGao et al. 19~9ll : iMcMillan. Hut fc Makino 199dll99lflHeggie fc Aarseth 1992 



Fregeau et al. 2003h . On long time scales, these processes lead to a slow (relaxation time) 
overall expansion of the cluster, with r v - u oc t 2//3 , a resu lt that follow s from simple consider- 
ations of the energy flux through the half- mass radius (IHenon 1965I ). 

While these processes are important to the long-ter m dynamical evolution (phase 3), their 

relevance is somewhat different during the first 100 Myr (jPortegies Zwart. McMillan fc Makino 20071 ). 
which is largely dominated by stellar mass loss (phase 2) and the segregation of the most 
massive stars. Often, their major effect is to enhance the rate of collisions and the formation 
of stellar exotica. We now consider in turn the following processes: binary heating ( §3.4. 1 j) . 
stellar collisions ( §3.4.2j) . and black-hole heating ( §3.4.31) . 



3.4.1 Binary Interactions 



Irrespective of the way they form, binaries are often described by dynamicists as either 
"hard" or "soft." A hard binary h as binding energy greater than the mean stellar kinetic 
energy in the cluster ( Heggie 19751 ): \Et\ > \{mv 2 ) ~ |(m)(t> 2 ), where (m) and (v 2 ) 1 ^ 2 are 
the local mean stellar mass and velocity dispersion. A binary with mass = mi + m2 and 
semi-major axis a b has energy = -Gm 1 m 2 /2aft, so hard binaries have a& < ah ar d, where 
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Here we implicitely assumed that m\ = m2 = (m) to derive the right-hand expression. The 
hard-soft distinction is helpful when discussing dynamical interactions between binaries and 
other cluster members. However, we note the definition of hardness depends on local cluster 
properties, so the nomenclature changes with environment, and even within the same cluster 
a binary that is hard in the halo could be soft in the core. 

The d ynamical significance of "hard" binaries (see Eq. l2"Tj) has been understood since 
the 1970s (jHeggie 19751 . iHills 19751 iHut &: Bahcall 19831 ) When a hard binary interacts with 
another cluster star, the resultant binary (which may or may not have the same components 
as the original binary) tends, on average, to be harder than the original binary, making 
binary interactions a net heat source to the cluster. Soft binaries tend to be destroyed by 
encounters. For equal-mass systems, the mean energy liberated in a hard-bina ry encounter i s 
proportional to Ej,: (AEf,) = jEb, where 7 = 0.4 for "res onant" intera ctions (jHeggie 19751 ). 
and 7 ~ 0.2 when wider "flybys" are taken into account (jSpitzer 19871 ). 

The liberated energy goes into the recoil of the binary and single star after the interaction. 
Adopting terminology commonly used in this field, we write the binary energy as E b = —hkT, 
(|mt> 2 ) and h ^> 1, so the total recoil energy, in the center of mass frame of 



where |fcT 



the interaction, is jhkT. In the center of mass frame, a fraction 



of this energy goes to 
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the single star (of mass m) and m ™^ m to the binary. For equal-mass stars, these reduce to | 
for the single star and ~ for the binary. Neglecting the thermal motion of the center of mass 
frame, we identify three regimes: 

1. If ^hkT < \mv es( ? = 2m(v 2 ) = 6kT, i.e. h < 45, neither the binary nor the single 
star acquires enough energy to escape the cluster. Binaries in this stage are kicked out 
of the core, then sink back by dynamical friction, in a process that we call "binary 
convection." 

2. If l^hkT > 6kT but \^hkT < 4m(v 2 ) = 12kT, i.e. 45 < h < 180, the single star 
escapes, but the binary is retained. We refer to such a binary as a "bully." 

3. If h > 36/7 = 180, both the binary and the single star are ejected. Such a binary is a 
"self-ejecter." 

These numbers are only illustrative. For a binary with components more massive than 
average, as is often the case, the threshold for bullying behavior drops, while that for self- 
ejection increases. 

Tab.|6]places these considerations in a more physical context. Note that, since the closest 
approach between particles in a resonant intera ction may be as little as a few percent of 



the binary semi- major axis (IHut fc Inagaki 19851 ). the hardest binaries may well experience 



physical stellar collisions rather than hardening to the point of ejection; and a collision 
tends to soften the surviving binary. Alternatively, before their next interaction, they may 
enter the regime in which internal processes, such as tidal circularization and/or Roche-lobe 
overflow, become important. The future of such a binary may be determined by the internal 
evolution of its component stars, rather than by further encounters. 



Table 6: Terminology (first column) and characterization (second and third columns) for 
the various stages of a binary (see text). The subsequent columns give the orbital separation 
a of a binary (in units of AU) with a total mass m\ + m 2 = mj = 10(m) or mj = 100(ra), 
in a cluster with a mass of M = 10 5 M Q and virial radius r vir = 1 pc and r vir = 10 pc. 
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The binary encounter time scale is t enc = (na(v 2 ) 1 ^ 2 ) -1 , where n is the local stellar 
density and a is the encounter cross section (see Eq. l24"l) . If we arbitrarily compute the binary 
interaction cross section as that for a flyby within 3 binary semi-major axes, consistent with 
the encounters contributing to the ISpitzer (19871 ) value 7 = 0.2, and again assume equal 
masses im b = 2m), we find 



to 



8ht r 



(22) 
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where we have used Eq. [13] and taken In A = 10. Thus the net local heating rate per binary 
during the 100% efficient phase (#1 above), when the recoil energy remains in the cluster 
due to "binary convection" is 



£bin = jhkT t e 



O.lkT/tr 



(23) 



that is, on average, each binary heats the cluster at a roughly constant rate. During the 
"bully" phase, the heating rate drops to just over one-third of this value. The limiting value 
of one-third is not reached since the ejected single stars still heat the cluster indirectly by 
reducing its binding energy by a few kT. For "self-ejecting" binaries, the heating rate drops 
almost to zero, with only indirect heating contributing. 

Binary-binary interactions also heat the cluster, although the extra degrees of freedom 
complicate somewhat the above discussion. If the binaries differ widely in semi-major axes, 
the interaction can be handled in the three-body approximation, with the harder binary 
considered a point mass. If the semi-major axes are more comparable, as a rule of thumb 
the harder binary tends to disrupt the wider one (jBacon. Sigurdsson fc Davies 19961 ). 

Numerical experiments over the past three decades have unambiguously shown how initial 
binaries segregate to the cluster c ore, interact, and suppo rt the core against further collapse 
( IMcMillan. Hut &: Makino 19901 ; iHeggie &: Aarseth 19921 ). The respite is only temporary, 
however. Sufficiently hard binaries are ejected from the cluster by the recoil from their last 
interaction (self-ejection, see Tab.E]), and binaries may be destroyed, either by interactions 
with harder binaries, or when two or more stars collide during the interaction. For large 
initial binary fractions, this binary-supported phase may exceed the age of the universe or 
the lifetime of the cluster against tidal dissolution. However , for low initial bin ary fractions, 
as appears to have been the case for the GCs observed today (IMilone et al. 2008|). the binarie s 
can be depleted before the cluster dissolves, and core collapse resumes ( iFregeau et al. 20031 ). 

Thus binary dynamics drives the evolution of the cluster while, simultaneously, the com- 
bination of cluster dynamics and internal stellar processes determine the internal evolution 
of each binary. This interplay bet ween stellar ev olution and stellar dynamics is someti mes 
referred to as star-cl uster ecology (Heggie 19921) or the binary zoo (IDavies et al. 20061 ) or 
stellar promiscuouity ( IHurley fc Shara 20021 ) . 



3.4.2 Stellar Collisions 

In systems without significant binary fractions — either initially or following the depletion 
of core binaries — core collapse may continue to densities at which actual stellar collisions 
occur. In young clusters, the density increase may be enhanced by rapid segregation of the 
most massive stars in the system to the cluster core. Since the escape velocity from the 
stellar surface greatly exceeds the rms speed of cluster stars (8 < 100 in Eq .[!]), collisions 
lead to mergers of the stars involved, with only small fractional mass loss (IBenz fc Hills 19871 . 
Freitag fc Benz 20011 ). If the merger products did not evolve, the effect of collisi ons would be 



to dissipate kinetic energy, and hence cool the system, accelerating core collapse ( IPortegies Zwart et al. 199i 
However, when accelerated stellar evolution is take n into account, the (time averaged ) en- 
hanced mass loss can result in a net heating effect (jChatterjee. Fregeau fc Rasio 20081 ). 
The cross section for an encounter between two objects of masses mi and m<i and radii 
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T\ and r 2 , respectively, is 



a 



nr 



1 + 



2G(mi + m 2 ) 



(24) 



(jHills fc Day 19761 ). where v is the relative velocity at infinity and r = r\ + r 2 . For r <C 
G(mi + m 2 )/-y 2 , as is usually the case for the systems discussed in this review, the encounter 
is dominated by the second term (gravitational focusing), and Eq. (24] reduces to 



a 



2nr- 



Gm 



(25) 



which is nearly independent of the properties of the other stars. 

Collisions between single stars are unlikely unless one (or both) of the stars is very large 
and/or very massive, or the local density is very high. Consider a large, massive star of 
mass m and radius r moving through a field of smaller stars, so the collision cross section is 
dominated by the properties of the massive star. The rate of increase of the star's mass due 
to collisions is 



dM 
~dT 



p c av ~ 1i\Gmrp c jv 
m \ ( r 







Pc 



/10km s _1 

ViooM y v.iooR y vio 6 M / P c 3 y v 



M Q /Myr. (26) 



Thus a massive star (m/M ~ t/R,© ~ 1) in a dense stellar core (p c ~ 10 6 M Q /pc 3 ) will 
experience numerous collisions during its ~ 3 — 5 Myr lifetime. 

The presence of primoridal binaries can significantly increase the chance of a traffic ac- 
cident. Hard binaries (see §3.4.ip are in the gravitational focusing regime, so the binary 
interaction cross section may be obtained by setting r = a in Eq. [25j Such an encounter 
will likely lead to the hardening of the binary and possibly the ejection of the single star 
and also the binary, as just discussed. However, it may also lead to a hydrodynamical 
encounter, i.e. a physical collision between two of the stars. It is quite l ikely that the 
third star will also be engulfed in the collision product (IFregeau et al. 20041 ) . Since bina- 
ries generally have semi-major axes much greater than the radii of the component stars, 
such binar y-mediated collisions play importan t roles in determining the stellar collision rate 
in YMCs ( jPortegies Zwart fc McMillan 20021 ). leading to significant numbers of mergers in 
lower-density, binary rich environments. Massive binaries in young dense clusters tend to be 
collision targets rather than heat sources ( jGiirkan. Freitag fc Rasio 20041 ) . 

In a sufficie ntly dense system, repeated stellar collisions can lead to a so-called "colli- 
sion runaway" (jPortegies Zwart et al. 19991 ). in which a massive star or collision product 
suffers repeated mergers and grows enormously in mass before exploding as a supernova 
flPortegies Zwart fc McMillan 20021 ; IPorteeies Zwart et al. 2004J : iGlirkan. Freitag fc Rasio 20041 ) 
This has frequently been cited as a possible mechanism for producing intermediate-mass 
black holes (IMBHs) in star clusters. However, while the dynamics is simple, numerous 
uncertainties in the stellar evolution and mass loss of the resultant merger product have 
been pointed out in the recent literature, suggesting that the net growth rate, and hence 
the final mass of the resulting IMBH, may be much lower than su ggested by purely dy- 
namical simulations — perhaps as little as a few hundred solar masses ( jYungelson et al. 20081 . 
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Glebbeek et al. 20091 . IVanbeveren et al. 20091 ) . Alternative formation mechanisms for more 
massive IMBHs involve gas accretion onto a see d ^100 M m black hole during a second round 

of star formation early in the cluster's lifetime (jVesperini et al. 20091 ). or r epeated collisions 

between stellar-mass black holes during the phase-3 evolution of the cluster (jMiller fc Hamilton 20021 ) 



3.4.3 Black Hole Heating 

An IMBH in a star cluster can be an efficient source of energy to the stellar system. Stars 
diffuse by two-body relaxation deeper and d eeper into the IMBH' s potential well, and even- 
tually are tidally disrupted and consumed (IBahcall fc Wolf 19761 ) . The energy lost during 
the process heats the system. The heating rate for an IMBH of mass Mbh in a cluster core 
of density p c and velocity dispersion v c is 



bh 



G 5 (m)p 2 c M BH \nA 



(27) 



Although cores are promising environments for the formation of IMBHs, they may not 
be the best place to look for evidence of massive black holes today. Dynamical heating 
by even a modest IMBH is likely to lead to a cluster containing a fairly extended core 
(iBaumgardt. Making fc Hut 20051 ) . Comparing the outward energy flux from stars relaxing 
inward in the (Bahcall-Wolf ) cusp surrounding the IMBH to the out ward flux implied by two- 
body relaxation at the cluster half-mass radius, iHeggie et al. (20071 ) estimate the equilibrium 
ratio of the half-mass (r h m) to the core (r c ) radius in a cluster of mass M. Calibrating to 
simulations, they conclude that for systems with equal mass, except of course the black holes 



0.23 



M 



M, 



BH 



3/4 



(2* 



Trenti et al. (20071 ) has suggested that the imprint of this process can be seen in his 



"isolated and relaxed" sample of simulated open clusters having relaxation times less than 
1 Gyr, a half-mass to tidal radius ratio rh m /r t < 0.1, and an orbital ellipticity of less than 
0.1. Roughly half of the clusters in this sample have core radii substantially larger than 
would be exp ected on the basis of simple stellar dynamics and binary heating. However, 
Hurley (20071 ) has argued that such anomalously large core to half-mass ratios may also be 
explaine d by the presence of a stellar-mass BH binar ies heating the cores of these clusters 
(see also lMerritt et al. 2004 iMackey et al. 20071 . 120081 ). Many of the results discussed above 
and much of our physical understanding of the dynamical evolution of star clusters have 
been developed and calibrated by means of simulations. 



4 The survival of star clusters 

The realization that the majority of star formation occurs in embedded clusters, whereas only 
a small fraction of stars in the Galactic disk currently reside in clusters (see JQ), indicates that 
most clusters and associations are relatively short lived; they disso lve on time scales compa - 



rable to the median age of open clusters in the solar neighborhood (IKharchenko et al. 20051 ). 
which is about 250 Myr (see t O]) . 
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Historically, studies of the lifetimes of star clusters have focused on open clusters in the 
Milky Way. The s carcity of open clusters with ages ^ 1 Gyr was reported independentl y 
in several studies ( von Hoerner 1958 ; van den Bergh 19571 ; Port. Kerr fc Westerhout 1958 ). 
and has been attributed to their short median lifetimes (about 250 Myr; IWielen 197ll ). 
rather than, say, a variation in the formation history or a detection bias toward young 
objects. These short cluster lifetimes have been ex plained as du e to the destructive effects of 
encounters with giant molecular clouds (GMCs) (jSpitzer 19581 ). A typical Galactic cluster 



with a density of ~ 1 M & pc~ 3 can survive the heating due to passing GMCs for about 
250 Myr. The remarkable agreement between the inferred mean lifetime and the expected 
survival time in the Galactic disk implicated GMCs as responsible for the destruction of open 
clusters (see §4.4p . The argument is further supported by the radial offset of the old (few 
Gyrs) open clusters toward the anticent er of the Galactic disk and away fr om the plane of the 
disk, where the density of GMCs is low ( Ivan den Bergh fc McClure 19801 ). as is illustrated in 
Fig. [H The galactic bulge an d spiral structure also contribute, though in lesser extend, to th e 
destruction of open clusters ([Weinberg 1994 iGieles. Athanassoula fc Portegies Zwart 20071 ) . 

By comparing the age distributions of clusters in the Magellanic clouds with those in the 
Milky Way Galaxy, the former pop ulation is found to be on averag e older and also more 
massive than the local population (jElson &; Fall 1985bl . iHodge 19871 ). The higher average 
cluster mass in the sample of Magellanic cloud clusters is a consequence of the difficulty in 
detecting low mass clusters. The apparent longer lifetimes of the Magellanic clusters could 
imply that more massive clusters tend to live longer, although the longer lifetimes could also 
be explained by the lower density of GMCs, the absence of bulges and spiral structures and 
the overall weaker tidal fields in the Magellanic clouds. However, in §4.4.21 we argue that 
GMCs are unlikely to play an important role in the early evolution of a YMC, due to their 
high initial density. The mechanism leading to the destruction of star clusters is therefore 
of major importance for understanding the evolution of star clusters from youth to old age. 



4.1 Simulating star clusters 

The formation, evolution and death of star clusters is a complex problem combining stellar 
dynamics, gas dynamics, stellar evolution and the evolution of the potential of the parent 
galaxy, all of which contribute to the cluster's appearance over the entire lifetime of the 
cluster (see also §3]). Over the past decade, significant progress has been made in modeling 
many of these processes simultaneously in numerical simulations of clusters during phase 2 
and phase 3. A striking omission is the self-consistent treatment of the interaction between 
stars and gas during phase 1. We focus here on simulations of phase-2 and phase-3 clusters, 
first describing treatments of stellar dynamics (see below), then turning to the inclusion of 
other physical processes into the mix. 

A broad spectrum of numerical methodologies is available for simulating the dynamical 
evolution of young star clusters. In approximate order of increasing algorithmic and physical 
complexity, but not necessarily in increasing numerical complexity, the various methods may 
be summarized as follows. 



Static Models are self-consiste nt potential-de n sity pairs f o r specific choices of phase - 
space distribution functions (IPlummer 19 111 . iKing 19661 . iBinney fc Tremaine 20081 ). 
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They have been instrumental in furthering our understanding of cluster structure, 
and provide a framework for semi-analytical treatments of cluster dynamics. However, 
they do not lend themselves to detailed study of star cluster evolution, and we will not 
dis cuss them fur ther here, instead referring the reader to the discussion in §1.3.2[ or 



to flSpitzer 19871 ). 



"Continuum" Models treat the cluster as a quasi-static continuous fluid whose phase- 
space distribution function evolves under the influence of two-body relaxation and 
other energy sources (such as binary heating) that operate on relaxation time scales 
(see Eq.dT 



• Monte Carlo Models treat some or all components of the cluster as pseudo-particles 
whose statistical properties represent the continuum properties of the system, and 
whose randomly chosen interactions model relaxation and other processes driving the 
long-term evolution. 

• Direct N-body Models follow the individual orbits of all stars in the system, auto- 
matically including dynamical and relaxation processes, and modeling other physical 
processes on a star-by-star basis. 

Much of our current understanding of the evolution of star clusters comes from detailed 
numerical simulations, and the above techniques are used for the vast majority if simulations. 
In order to appreciate some of the details presented in §4.21 to §4.4.21 it may be important to 
have some rudimentary understanding of this range of techniques. However, since this is a 
rather technical discussion, we present it in an appendix 



4.2 Early violent gas expulsion 

Possibly the greatest discrepancy between star cluster simulations and observations lies in 
the first few million years of the evolution (phase 1 in §3D- Real star clusters are formed 
in a complicated interaction between gas and gravity, which is imperfectly understood. 
Once a primordial gas cloud starts to condense into stars dynamical evolution also be- 



gins (IBate. Bonnell &: Bromm 20031 ). At the end of the star formation process, probably 
brought about by the developing winds of the most massive stars or the first supernovae, 
the residual gas is ejected from the protostellar cluster. The gas expulsion phase is expected 
to be short — on the order of several dynamical times (Eq. ll2p — and places the remaining 
stellar population in a super-virial state, making the young cluster vulnerable to dissolution. 
The sharp decrease in the number of young and embedded star clusters at an age of a few 
megayears is thoug ht to be a conseque nce of this early process, and is often referred to as 



"infant mortality" (ILada fc Lada 20031 ). 



4.2.1 Theoretical considerations 

The total mass M = M g + M* of the primordial cluster contains contributions from stars 
M* and gas M g . In virial equilibrium, the rms velocity of the cluster is 

<^ 2 > = ^— (29) 

z,/ vir 
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which in observable quantities becomes 



2 GM 

CTiD = , (30) 

V r efi 

with 77 « 10 (see 3TX2]) . 

The response of the cluster to the loss of the residual gas depends on the gas expulsion 
time scale t exp relative to the dynamical time scale td yn of the cluster (see Eq. fTOl) . For many 
young clusters t cxp <C t<jyn> m which case removing the gas shocks the cluster. In the extreme 
case, the positions and velocities of the stars remains fixed during the gas expulsion phase, 
and the response of the cluster to losing a fraction of gas, f e = M*/M, can be calculated 
under the assumption that the mass loss is instantaneous. The rms velocity of the stars 
immediately before and after gas loss is then given by Eq. [29j the stellar positions are also 
unchanged. As a cons equence, th e cluster expands and re-establishes equilibrium at a new 



virial radius given by (IHills 1980h 



T vi r fe _ (31) 



r vir (t = 0) 2/ e -l 

For high star formation efficiency (f e ^ 0.5) the arguments leading to Eq.[3T]may be reason- 
able, as in this case a relatively small fraction of the stars is lost after the gas is expelled. 
However, the energy argument is too simple to determine the survival probability if f e ^ 0.5. 
Losing more than half of the total mass (f e < 0.5) by explosive gas expulsion is devastating 
for the cluster, leading to its complete dissolution in a few dynamical time scales. A low 
star formation efficiency may explain the majority of disrupted young and embedded clus- 
ters, but even for f e ^ 0.5 some small portion of the cluster can in practice remain bound. 
There are several independent arguments for the survival of embedded clusters even for a 
star formation efficiency as low as f e ^ 0.1. 

The most important argument against the above simple analysis (Eq. l3~l|) is the fact that 
the time scale for gas expulsion, t exp , in practice is short but not instantaneous. This time 
scale should not be compared to the (global) half-mass crossing time, but rather to the local 
dynamical time, which depends strongly on the distance to the cluster center. For example, 
the dynamical time scale in the cluster core has t corc /t dyn = (phm/Pc) 1 ^ 2 and this fraction 
ranges from ;$ 0.01 for a concentrated cluster (Wo ~ 12, or c > 2.7) to ~ 0.2 for a shallow 
potential (c < 1 or Wq ~ 5). For concentrated clusters (c £ 2), the gas expulsion occurs more 
or less instantaneously for stars in the outskirts, but stars in the core respond adiabatically 
to the loss of gas. If t exp > td yn , the cluster is likely to survive; in particular, the core may 
respond with just a sligh t expansion (of at most a factor of 2), even in extreme cases, f e £ 0.1 



flGever fe Burkert 200lh . 



Further complicatio ns arise from the clumping of the ga s and stellar distributions in 



the embedded cluster (iFellhauer. Wilkinson fe Kroupa 20091 ). It is likely that the radial 



dependence of the star formation process causes the central part of the cluster to be de - 



pleted in gas, whereas the outskirts are relatively gas rich ( IBonnell. Bate fe Vine 20031 ). 



This radial variation of the star formation efficienc y, combined with the process of com- 



petitive star formation ( IBonnell. Bate fe Vine 20031 ). may actually render the cluster sub 



virial after the gas is ejected, for example if stars formed in the collapsing cloud are dy- 
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namically cold ( ILada. Margulis fe Dearborn 1984 ). as suggested by the outcome of tu rbu- 
lent self-gravitating simulations of cluster formation ( jOffner. Hansen fe Krumholz 20091 ). In 
that case, Eq. ED depends on the fractional deviation from virial equilibrium g vir , to become 
vir/Vvir(t = 0) = fe/ (2 f e — q^ r ), i.e. th e condition for complete disruption becomes f e = q^ iT /2 



( IProszkow et al. 20091 . iGoodwin 20091 ) . The cluster survival prob ability then depends on th e 
entire star formation process, not just on its overall efficiency (IGoodwin fe Bastian 20061 ). 
A further deviation from the simple formulation comes from the effect of high-velocity es- 
capers, which can carry away a considerab le fraction of the cluster's bi nding energy, leaving 
the remaining stars more strongly bound (IBaumgardt fe Kroupa 20071 ). 

Thus the survival probability of an embedded cluster cannot be determined by a single 
parameter, such as the star formation efficiency, as the complete formation process of stars, 
clumps of stars, and the entire cluster comes into play. The process whereby stars form in 
massive star clusters is still poorly explored terrain within astrophysics, and the formation 
of sub-clumps and clusters is even less well charted. 



4.2.2 Observational constraints 

What sounds convincing from a theoretical standpoint is often very hard to support with 
observations. There are a number of interesting observational indications of infant mortality, 
and of the associated dissolution time scales, but since the embedded phase is short (~ 
1 — 2 Myr), many parameters are poorly constrained by observations, and the interpretations 
of models tend to be sensitive to assumptions made about the initial conditions. 

Many observed young clusters, particularly the extragalactic population in Tab. [4] appear 
to be super- virial, which naively would lead to the early termination of the cluster's existence. 
This is most easily seen in the higher value of their dynamical mass Md yn compared to the 
photometric mass M p h t- The former is derived from measurements of the velocity dispersion 
and the radius of the cluster and the use of Eq. HI The latter is d erived from the tota l 
luminosity, calibrated to single stellar population models (see §|2]). iBastian et al. (20061 ) 



determine M p hot and Md yn from a compilation of 19 clusters and find that both independent 
mass estimates are consistent for the somewhat older ( ^ 50 — 100 Myr) clusters , but for 
young (~ 10 Myr) star clusters they find M dyn > M phot . IGoodwin fe Bastian (20061 ) explain 
this as a signature of the primordial gas expulsion, and o f the process of infan t mortality. 

In Fig. [12] we present an updated version of Fig. 5 of IBastian et al. (20061 ). showing the 
ratio of light to dynamical mass of 24 clusters, taken from Table HJ each of which is about 
lOmegayear old (see $2]). The age range is quite narrow because the red supergiant phase, 
which starts around 10 Myr, makes these clusters extremely bright (especially in the near 
infrared), whereas younger clusters are still heavily obscured. 

Many of the clusters in Fig. [12] and Tab.H] appear to have dynamical masses too high 
for their luminosities, and the i ntuitive explanation is th at these clusters are expanding 
and are possibly even unbound (IGoodwin fe Bastian 20061 ). The time needed for a cluster 
to completely dissolve, or to find a new virial equilibrium after gas expulsion, is about 



20 — 30td yn (see for example Fig. 8 in IBaumgardt fe Kroupa 20071 ). where t dyn is the 



initial dynamical time, when the gas and the stars are still bound. As a result, to 'catch' 



an unbound or expanding cluster at 10 Myr, the initial tdyn should be 



corresponds to an half-mass density of p nm < 300 M® pc 3 (Eq. [12] and phn 



> 0.5 Myr. This 
3M/8vrr nm 3 ). 
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Figure 12: Left: Light to dynamical mass ratio for the clusters in Table H] th at have both 



quantities available. Photometric evolution from the lBruzual fc Chariot (20031 ) single stellar 
population models, using a Chabrier IMF, in the V and K bands is indicated by the full 
and dashed lines, respectively. Right: The measured velocity dispersion squared, presented 
as a function of a prediction for this quantity based on Eq. HJ The solid line is a prediction 
of the effect of binaries on M dyn , with 1-a and 2-cr variations due to stochastic fluctuations 



shown as dashed lines based on the results of iGieles. Sana fc Portegies Zwart (20091 ). 



Clusters with shorter initial tdyn (higher density, see Eq.[[2]) have expanded into the field, 
or found a new equilibrium a few megayears after gas expulsion, and are not observable as 
super-virial clusters at 10 Myri 

The initial density of the clusters in Fig. [T2] is unknown, but is likely to have been higher 
in the past that it is today. In addition we may attemp to estimated their initial densities 
from their current densities. In Fig. [9] we show the radii and masses of the clusters under 
discussion, together with lines of constant p hm . The present day densities of the clusters 
younger than 10 Myr are phm ~ 10 3±1 M & pc -3 . The densities in the embedded phase were 
0(1/ /*) higher — where the reduction in cluster mass accounts for a factor l/f e and the 
consequent adiabatic expansion a factor l// 3 (see Eq. [331) . 

The initial dynamical times were therefore a factor of 1/ / e 2 shorter than tdyn after phase 
1. Based on their physical ages of ~ 10 megayears and the fact that the gas ejection phase 
does not last beyond the moment of the first supernova (within ~ 3 Myr, see §4.ip . these 
clusters h ave evolved for at least 10/ to 100/ / 2 initial dynamical times, and therefore must 



be bound (jMcCrady. Gilbert fc Graham 20031 : iMcCrady fc Graham 20071 ). and the observed 



discrepancy between M p hot and Md yn cannot originate from the overall expansion of the 
cluster. 



Fleck et al. (20061 ) showed that the constant r\ that relates o"id and r eS to M dyn (Eq. I3"U1) 



e The models of Goodwin fc Bastian (2006) start with a density in the embedded phase of 



60 Mq pc (tdyn ~ 1 Myr), which is why their clusters are super-virial for at least 25 Myr. 
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is up to a factor of three higher when the cluster is mass segregated. This results in a value 
of Md yn that is too low compared to the true mass when this is not taken into account. So 
mass segregation is not a plausible explanation for the large downward spread of points in 
Fig. EI 

We are still left with the question why the observed velocity dispersions in some of these 
clusters are higher than would be expected from the virial theorem. Several independent 
and implicit assumptions enter the derivation of M& yn and M p h ot , and each of them could 
be wrong. The stellar mass function, for example, could be bottom-heavy, i.e. steeper than 
Salpeter or with an excess of low-mass stars. Such a mass function would result in a velocity 
dispersion in virial equilibrium higher than that of a cluster with a Salpeter IMF, with little 
effect on M phot . However, to explain the observed discrepancy, the cluster mass function must 
deviate substantially from the canonical mass functions. We do not favor this conjecture, 
since the only star cluster for which the mass fun ction was once anticipated to be deficient 



in low- mass stars, the Arches (IStolte et al. 20 05). turns out to have a rather normal mass 



function at least down to 1 M ( IKim et al. 20061 



Another interesting possibility is provided by the binarity of red supergiants, which dom- 
inate the observed luminosity. A relatively high fraction of hard binaries (see §3.4. ip leads to 
an overestimate of the cluster velocity dispersion due to the contribution from their internal 
orbital motion. This leads to an overestimate of a id, and therefore of the mass of the cluster 



(IKouwenhoven fc de Grijs 20081 ). For typical op en clusters, with am ~ 1km s , this can 



only account for a factor of ~ 2 increase of Md yn (IKouwenhoven fc de Grijs 20081 ). However, 
young star clusters are dominated by ~ 13 — 22 M red supergiants, and a binary fraction 
of ~25% among these stars could explain an apparent dynamical mass of up to an order 



of magnitude more than the photometric mass (IGieles. Sana fc Portegies Zwart 20091 ). As a 
consequence the discrepancy between Md yn and M p hot is larger for clusters with high stellar 
velocity dispersion, or with a small ratios of M p hot Aeff (see Eq. 0j), which is consistent with 
the observations. The effect of binaries on am and the ratio of M^ jn / M v \ lot is presented in 
Fig. [T2J 

Star clusters with phm ^ 100 M pc~ 3 at 10 Myr, such as those listed in Tab. H] and shown 
in Fig. [121 have survived the primordial gas expulsion phase and should be considered bound, 
stable, and likely to survive for a long time (see Eq. [T9|) . They enter the next evolutionary 
phase, described in §4.3[ during which stellar mass loss dominates the evolution of the cluster 
(phase 2). 

4.3 Stellar mass loss 

Clusters that survive phase 1 (the embedded phase, §4.2j) continue to lose mass through 
stellar evolution. During phase 2, the most massive ( £ 50 M Q ) stars leave the main-sequence 
within 4.0 Myr, and lose about 90% of their mass by the time they collapse to a black 
hole following a supernova. A 5M Q sta r loses ^ 80% of its mass by the time its core forms 



a white dwarf at about 100 Myr. For a lKroupa f2001l ) IMF between 0.1 M and 100 M the 



total cluster mass decreases by roughly 10%, 20%, and 30% during the first 10, 100, and 
500 Myr. The impact of this will be discussed below. 
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4.3.1 Theoretical considerations 



The time scale for mass loss depends on the mode in which it is achieved; supernova ex- 
plosions, Wolf-Rayet winds and AGB expulsion result in high mass loss rates, whereas the 
general mass loss for an older stellar population is relatively slow. When the time scales for 
mass loss by stellar evolution is considerably longer than t^ ya , the cluster responds adiabat- 
ically, expanding through a series of virial equilibria. For small f e , Eq. EH reduces to 



and therefore 



5M 



M{t = 0) 
M ' 



(32) 



(33) 



r vi r(t = 0) 

Even losing half the mass by slow stellar evolution (which for the canonical IMF would not 
occur within a Hubble time), the cluster would expand by only a factor of two. In the 
instantaneous approximation ( §4.21 and Eq.ETTl). such mass loss would lead to the dissolution 
of the cluster. 

In reality, the situation is more complicated, in particular because of the connection 
between dynamical evolution and stellar mass loss. For real clusters the expansion due to 
stellar mass loss is considerably more severe than suggested above, and can even result in 
complete dis ruption if the cluster is mass segregated befor e the bulk of the stellar evolution 
takes place ( jVesperini. McMillan fe Portegies Zwart 20091 ). Even an initially unsegregated 
cluster can still undergo mass segregation during the period w hen the residual gas is being 
ejected, and certainly during the early evolution of its stars (lApplegate 19861 ) . which can 
also lead to enhanced expansion at later times. The expansion of a mass segregated cluster, 
however, will not be homologous, as the massive (segregated) core stellar population tends 
to lose relatively more mass than the lower-mass halo stars. The result is a more dramatic 
expansion of the cluster core, with less severe effects farther out. 

This effect is illustrated in Fig. [TBI which presents the results of an isolated N = 128/c 
body sim ulation, run on a GRAPE-4 (IMakino fe Taiji 19981 ) using the star lab software envi- 
ronment (IPortegies Zwart et al. 20011 ) . The figure shows the evolution of the core radius with 

and without stellar evolution. Without stellar evolution the core tends to shrink, and eventu- 

ally reaches core collapse, whereas with stellar mass loss the core expands (IPortegies Zwart. McMillan fe Ma 

The combined effects of mass loss by stellar evolut ion and dynamical evolution in the tidal 



field of the host galaxy have b een extensively studied (IFukushige fe Heggie 19951 . iTakahashi fe Portegies Zw 
Baumgardt fe Making 20031 ). These studies show that when clusters expand to a radius of 
~0.5 rj they lose equilibrium and most of their stars overflow rj (see §1.3.2p in a few crossing 
times. 

We conclude that in phase 2 the overall evolution of the cluster is completely dominated 
by expansion due to stellar mass loss. However, since most of the mass loss comes from 
massive stars in the core, the core expansion is considerably larger than expected for the 
global mass function. Phase 2 lasts until the response of the cluster to stellar mass loss 
deminishes and from that moment the cluster can continue to expand until it completely 
dissolves or until the cluster core starts to contract again due to internal dynamical effects, 
at which point phase 3 begins. 
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Figure 13: Evolution of the core radius during phase 2 of an iV-body simulation (N = 
128fc, rvir = 3.2 pc, King Wo = 12, Mass function is Salpeter between 1 M and 100 M 



Portegies Zwart. McMillan fc Makino 20071 ). After about 3Myr stellar mass loss dominates 



the evolution of the cluster core radius. The dotted curve (top curve) is the result of the 
full simulation, with both stellar evolution and binary dynamics included; the irregular 
dashed curve is calculated with the same initial realization but without stellar evolution. 
The wiggely solid curve (bottom) is calculated without stellar mass loss and without binary 
dynamics by collapsing the binaries in single objects. The smooth dashed line shows the 
expected expansion of the core, assuming adiabatic mass loss for a Salpeter initial mass. 
The smooth solid curve is computed assuming mass loss by stellar evolution from a Salpeter 
mass function with a lower limit of 15 M Q , rather than the 0.1 M used in the simulation. 



4.3.2 Observational constraints 

Rapid dissolution due to stellar evolution mass loss occurs when rh m /^j ~ 0.5. All the 
clusters in Fig. [9] have, within a factor of tw o, rhm/r,] ~ .03, suggesting that they are 



probably all stable against stellar mass loss. iPfalzner (20091 ) recognizes two evolutionary 
sequences in young Galactic star clusters, from which she draws a similar conclusion. The 
first sequence of dense "starburst clusters", containing the Arches cluster, NGC 3603 and 
Trumpler 14, starts at a density of ~ 10 5 M Q pc~ 3 at an age of a few megayears. These 
clusters appear to expand at constant in mass, up to an age of 10-20 Myr, where we find 
the red supergiant clusters RSGC01 and RSGC02. At that age the cluster density has 
dropped to ~ lO 3 M pc~ 3 . The second sequence of "leaky clusters" (which is identical to 
our definition of associations, see §2.ip starts at the same age but with much lower densities 
of ~ lOM pc~ 3 , and expand while M oc l/rh m to densities comparable to the field star 
density. The clusters in Fig. [9] may be compared with the red supergiant clusters in the 
Milky Way, i.e . the end point of the dense cluster sequence. The associations discussed by 



Pfalzner (20091 she refers to them as "leaky" clusters) and listed in Table |2] have dynamical 
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Figure 14: Evolution of observed core radii values of clusters, compared to the results of 
an iV-body simulation including mass loss due to stellar evolution. For the simulation, the 
projected surface brightness profile was constructed and r c was determined using and Eq. [5] 
and Eq. [8] (full line). The dashed line shows the theoretical r c from Eq. [7] discussed in §1.3.21 



times that exceed the cluster age, and are expected to be unstable against mass loss by 
stellar evolution. 

In £J2] we di scussed the remarkable increase in the obsery ed core radii of clusters with 
ages 1-50 Myr ( jMackey fc Gilmore 20031 . iBastian et al. 20081 1 . Earlier studies argued that 
this effect was the result of early gas expulsion (phase 1), but as discussed in $2] the time 
scale for finding a new virial equilibrium after gas expulsion is too short to affect the growth 
of r c over such a long period. Mass loss from the young stellar population (phase 2) can 
contribute to some ex tent to the observe d trend, but its effect is probably too weak to 
explain it completely (jMackey et al. 20071 ). A more likely solution is dynamical heating 
from relatively massi ye objects, such as massive stars and stellar mass black holes, sinking 
to the cluster center ( IMerritt et al. 20041 ). 

An additional effect not yet discussed is the difference between the observed r c , usually 
resulting from a fit to the surface brightness profile (Eq. [5]), and the 3D dynamical core 
radius described in §1.3.21 If a young star cluster is mass segregated, possibly already 
from its formation process, then the light is d ominated by the few massive stars in th e 
core, which can lead to an underestimate of r c ( jFleck et al. 20051 . iGaburov &: Gieles 20081 ) . 



Mackey et al. (20071 . 120081 ) showed that when taking this into account, a remarkable increase 
of r c is "observed," while the "real" core radius changes less. 

We illustrate this in Fig. HH The same data points as in Fig. [8] are shown, together with 
the results of an iV-body simulation (lines), where the theoretical value of r c , as defined in 
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51] and the value of r c resulting from a fit of the EFF87 profile (Eq. [5]) to the projected light 
of the simulate d cluster. The c luster consists of N = 65 536 single stars initially distributed 
according to a lPlummer (19111 ) profile with r vir = 2pc. Before stellar evolution was turned 
on, we evolved the cluster for 100 td yn , which corresponds to ~ 0.1 t ID ., to mimic some 
degree of primordial mass segregation, in which stars with masses > 5 M are more centrally 
concentrated than the less massive stars (Eq. [20]) . Since the massive stars dominate the 
light, the observed r c is almost a factor of three smaller than the 3D version at t — 0. The 
observed r c increase by nearly a factor of ten in a few tens of megayears, while the 3D core 
radius expands only by a factor of three. After ~ 30 Myr the two quantities roughly agree. 
In the simulation r hm increased only by 60% due to the stellar evolution mass loss. 

Primordial mass segregation has a profound effect on the evolution of a star cluster, but 
possibly much more relevant for this review is its consequences for cluster observations. The 
assumption that a young ( ^ 10 Myr) cluster is not mass segregated, when in reality it is, 
can dramatically alter observationally derived quantities, such as the cluster size, velocity 
dispersion, density profile, and central density. 



4.4 External perturbations and evaporation 
4.4.1 Theoretical considerations 



One important external disruptive factor, first considered by lSpitzer (19581 see §4.2p . is en- 
counters between clusters and giant molecular clouds. Since GMCs are typic ally more mas - 
sive than clusters, the cluster is more affected by an encounter than the cloud (jTheuns 199ll ). 
The cluster lifetime due to heating by passing clouds is inversely proportional to the volume 
density of molecular gas, p gas , and proportional to the density of the cluster: 



t 
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lGyr 



0.03 M pc 



-3 



Phi, 



10 M pc" 



(34) 



This result is typical for disruption by external tidal perturbations operating on short time 
scales ( ^ tdyn), also known as tidal shocks a nd can also be applied to passages thr ough the 
disc fe.g. lOstriker. Spitzer fc Chevalier 1972), bulge (e. g. Gnedin fc Ostriker 1997 ) and spi- 
ral arms ( IGieles. Athanassoula fc Portegies Zwart 20071 ). Here 0.03 pc~ 3 is the molecula r 
gas density in the solar neighborhood and the c onstant is tak en from lGieles et al. (2006a ). 
which is an update from the seminal result by ISpitzer (19581 ). The dependence of tdis on 
Pg as indicates that the lifetimes of star clusters scales roughly inversely with the observable 
surface density of molecular gas, E gas . 

This result enables us to make order of magnitude estimates of the lifetimes of clusters 
in other galaxies. In spiral galaxies, GMC encounters are especially frequent during the 
early stages of cluster evolution, since clusters form in the thin gaseous disk where p gas is 
high. Older clusters are typically associated with the thick disk, where p gas is low and GMC 
encounters are less frequ ent. Since y oung (<1 Gyr) clusters in spiral galaxies have only a 
small range in radii (e.g. lLarsen 20041 ). more massive clusters tend to have higher densities, 
making them less vulnerable to encounters with GMCs, which e xplains their longer lifetimes 
compared to their lower-mass counterparts (IGieles et al. 2006a ). It is not clear whether the 
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lack of a mass-radius relation is a universal property imprinted at cluster formation, or the 
result of evolution. 



4.4.2 Observational constraints 



If GMC encounters are a dominant disruption process, then it is important to understand 
the mass-radius correlation of clusters, since t^ c oc phm ( §4.4. ip . If clusters form with a 
constant density, i.e. a mass-radius relation of the form rh m oc M 1 / 3 , then their destruction 
time due to GMC encounters is independent of cluster mass. Additional complications arise 
from the time dependence of a mass-radius relation. Clusters older than those shown in Fig. M 
For young clusters (< 10 Myr) there seems to be some positive correlation between mass and 
radius, roughly consistent with a density of 10 3±1 M Q pc~ 3 . Older clusters (> 10 Myr) do not 
seem to exhibit any correlation between mass and radius, consistent with a near-constant 
radius, i.e. r^ m = constant, which could be the consequence of evolutionary effects, such as 
mass segregation and stellar mass loss. This would lead to a destruction time scale by GMCs 
(^dis 10 ) that is dependent of the cluster mass: t^ c oc M/r\ m . In any case, the youngest 
YMCs listed in Tab. [2] are unlikely to be rapidly destroyed by passing GMCs, since t^is 10 
exceeds a Hubble time due to their high densities (see left panel of Fig. [9]). 

While considering mass loss from star clusters it is convenient to distinguish between two 
fundamental processes: evaporation and tidal stripping. Evaporation is the steady loss of 
stars from the cluster driven by the continuous repopulation by relaxati on of the high- velocity 



tail of the M axwellian velocity distribution (see §3.3. 1[ Eq.[14]and, e.g. lAmbartsumian 1938 



Spitzer 19401) . Th i s process has been the subject of numerous comprehensive numeri cal stud- 
ies flSpitzer 19871 lAarseth 20031 . Irleggie k Hut 20031 iBaumgardt k Makino 2o"o3h . Tidal 
stripping is the prompt removal of stars that find themselves outside the cluster Jacobi 
radius (rj, see ll.3.2p due to internal processes such as stellar mass loss or a change in the 
external tidal field, for example as the cluster approaches pericenter in its orbit around its 
parent galaxy. On a ~100 Myr time scale, and for clusters with masses ^ 10 4 M Q , relaxation 
is unlikely to be important, and so tidal stripping dominates the cluster mass loss. 

Until the 1990s, the main targets of cluster disruption studies were the open clusters 
in the Milky Way and the YMCs in the Magellanic clouds. Since then, HST observations 
have established the propert ies of large p opulations containing mo re massive clusters in 



quiescent spiral gal axies (e.g. lLarsen 20021 ). inte racting galaxies (e .g. IWhitmore et al. 1999 



Bastian et al. 20051 ). and merger products (e.g. iMiller et al. 19971 ). The primary tool used 
in studies of cluster disruption is the cluster age distribution. Different groups give different 
weights to the various factors described in §4.21 in interpreting the results, but empirical 
cluster disruption studies follow one of two basic models, in which the disruption is either 
environment dependent or universal. 

In the environment dependent disruption model, the dissolution time follows a simple 
scaling relation with cluster mass and environment, followin g the age and mass distribu - 
tions of luminosity-limited cluster samples in different galaxies ( IBoutloukos fe Lamers 20031 ). 
Variations in dissolution time scales are explained by difference s in the tidal field s trength 
( iLamers. Gieles k Portegies Zwart 20051 ) and the GMC density ( iGieles et al. 2006a ). 

The universal disruption model assumes that internal processes dominate the cluster 
disruption and that roughly 80-90% of all remaining clusters are destroyed during each 
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Figure 15: Age distribution of SMC clusters based on the catalog of IChiosi et al. (20061 ). 
The sample is split into small and large clusters/associations, with the boundary at a radius 
of 6 pc. The histograms are made using a 0.5 dex bin width with different starting values 
(boxcar averaging). 



decade in age, resulting in a (mass limited) age distribution that declines inversely with time 
(oc t~ x ). The main assumptions are that the majority (80 — 90%) of clusters dissolve within 
a few hundred megayears and that the disruption rate does not depend on mass. This model 
however , is calibrated against the cluster population in the Antennae galax ies, as discussed 
in &j~2l fcaJl. Chandar fc Whitmore 20051 : Iwhitmore. Chandar fc Fall 2007h . and while it is 
quite consistent with these observations, it is not clear how applicable it is to other galactic 
environments. 

The various cluster dissolution models have led to som e controversy, resulting in a number 
of spirited discussions at conferences and in the literature. IChandar. Fall fc Whitmore (20061 ) 
demonstrated that the age distribution of SMC clusters declines oc t~ 85 , which is consis- 
tent with their resuls for the Antennae. iGieles. Lamers fc Portegies Zwart (20071 ) are able 
to reproduce these results only if they impose a n incompleteness in the de tection of clus- 
ters, w hich is consistent with the arguments o f Ide Grijs &: Goodwin (20081 ). As a conse- 
qu ence IGieles. Lamers &: Portegies Zwart f20071) conclude, using the same data sample used 
by IChandar. Fall &: Whitmore (20061 ). that the age distribution of massive ( £ 1O 3 5 M ) 



clusters in the LMC younger than a few hundred megayears is not aff ected by internal pro- 
cesses, contradicting the findings of IChandar. Fall fc Whitmore (20061 ). 

The conclusions of Gieles. Lamers fc Portegies Zwart (20071 ) are supported by Boutloukos fc Lamers 
who show that for a constant formation rate of clusters and without disruption, the age dis- 
tribution of a luminosity-limited cluster sample declines oc t - ^ -1 ), where —j3 is the index 
of the cluster initial mass function. The parameter ( describes how the cluster fades with 
age, due to stellar evolution: F\(t) oc , where F x is the flux at wavelength A of a clus- 
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ter with constant mass. For the U and V-bands £ ~ 1.0 and 0.7, respectively, which for 
(3 = 2 results in t~ K , with 0.7 < k < 1.0, due to fading alone. If the distribution of cluster 
masses is described by a Schechter function (see Eq.[TS] in §2.4.2|) the age distribution of a 
luminosity-limited sample that is not affected by disruption is as steep as 0.9 < k < 1.4. 

The internal disruption model for the age distribution relies on the assumption that 
the cluster formation history is constant within the age range considered (a few hundred 



megayears) . iBastian et al. (20091 ) demonstrate that in the Antennae an increase of the cluster 
formation rate as predicted by galaxy merger models could alleviate the need to invoke longer 
duration cluster disruption to explain the decline in the number of clusters with age. 

The discussion of the distribution of the number of clusters with age is complicated by 
the distinction between associations and dense clusters (see $2] and §4.2.2p . which may be 
hard to distinghuish at large distances. This is illustrated in Fig. CG3, which shows the ag e 



distribution of clusters and associations in the central region of the SMC (IChiosi et al. 20061 ). 
The sample is divided into two subsamples of ~ 200 clusters each, based on size. The age 
distribution of the large (r e g > 6pc) clusters falls off much more rapidly than that of the 
more compact clusters. We note that this is consistent with associations being disrupted 
by GMC encounters, since they have (on average) much lower densities (see right panel of 
Fig. ED and are therefore destroyed much quicker (Eq. I3~4"|) . The median radii of the two 
samples are 4.5 pc and 9 pc, respectively, but at a distance of 20Mpc it would be very 
difficult to tell these two groups apart. We suggest that by relaxing the cluster size in the 
observationally selected sample one includes more short-lived associations, which if resolved 
should probably not be considered genuine star clusters. 

4.5 Early evolution and the luminosity function of old globular 
cluster systems 

The LFs of old globular cluster systems is very different from the (power-law) LFs of YMCs 
discussed in §2.4.11 When approximating the L F of old glob ular clusters in the Milky Way 



by a Gaussian, the peak is at My —7.4 mag ([Harris 20011 ). For a constant mass-to-light 
ratio of M/L ~ 2 the mass function, often referred to as globular cluster mass function 
(GCMF), which in fact is a present-day mass function. The GCMF is, like the LF, also 
peaked or bell-shaped (when using logarithmic mass bins). The characteristic mass, or turn- 
over mass, is Mpo ~2x 10 5 M . It has been argued that old globular clusters fo rmed with 



this -bell-shaped- mass function ( Fall &: Rees 19851 . Parmentier fc Gilmore 2007 ). which is 



completely different than the power-law distribution we see for YMCs today. The discrepancy 
between the bell-shaped curve for old globular clusters and the power-law mass function 
for YMCs then could originate from the long timescale phase-3 dynamical evolution, this 
evolution could severre the relatively low-mass (M Mtq) clusters more whereas it leaves 



the high-mass (M ^ Mro) clusters u naffected (jElmegreen fc Efremov 19971 . IVesperini 19971 . 



Baumeardt 19981 . iFall fc Zhang 2001 



The value of Mro, however, is remarkably constant between galaxies (IJordan et al. 20071 ) 



and even within a si ngle galaxy it hardly varies with the distance to the galaxy center 



(IVesperini et al. 20031 ). This is hard to reconsile with the models that start with a power-law 



mass function and seek to explain the present-day form of the GCMF by phase-3 dynamical 
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evolution (jVesperini et al. 20031 ). 

For this review only those solutions to this con undrum are relevant that play an important 



role in the early evolution (i.e. Phase 1 and 2). iBaumgardt. Kroupa fe Parmentier (20081 ) 
provide such a solution; they argue that relatively low-mass clusters are more easily dissolved 
as a result of phase-1 gas expulsion, leaving the most massive clusters unaffected. They argue 
that the binding energy of embedded clusters scales with the square of their total masses (gas 
and stars), whereas the energy input from stellar winds and supernovae depends linearly on 
the mass of all the stars in the cluster. If the star formation efficiency (f e ) is independent of 
the mass of the primordial molecular cloud, more massive clusters require more time to clear 



all the residual gas from the natal cloud. IBaumgardt. Kroupa fe Parmentier (20081 ) find 
that clusters with M 10 5 M Q loose their residual gas on time-scales much shorter than 
a crossing time, which effectively shocks the cluster and drives it to dissolution (see §4.2j) 
More massive clusters respond adiabatically to the loss of the mass of the cloud ( §4.3. lj) . 
resulting in a higher survival probability. The weak dependence of Mtq on environment al 
factors emerges quite naturally from this model IBaumgardt. Kroupa fe Parmentier (20081 ). 

An alternative explanation for the evolution of the ma ss function of YMCs compared to 
that of the GCMF is presented by lVesperini fe Zepf (20031 ). They argue that phase-2 stellar 
mass loss leads to a near universal dissolution rate of low- mass clu sters. Old globular clusters 
of low mass tend to be less concentrated than high- mass clusters ( McLaughlin 20001 ) . If this 
trend is primordial, stellar mass l oss would be more effective i n destroying low- mass c lusters 
compared to high-mass clusters ( IChernoff fe Weinberg 199CM . IVesperini fe Zepf 20031 ). This 
mechanism is, as the phase-1 outgassing, independent of environmental factors, like the dis- 
tance to the galactic center, and would therefore effectively provide a satisfactory explanation 
of the weak dependence of Mxo to the environment. 

Each of the three solution to this conundrum have their advantages and disadvantages. 
It may be evident that this discussion has not yet resulted in a satisfactory solution to the 
problem, but it promises to continue to be an exciting discussion for the foreseeable future. 



5 Lusus Naturae 

Old globular clusters are of interest because of their old age, their assumed relatively ho- 
mogeneous populations, their relative isolation in their parent galaxies, and because of their 
abundance of unusual objects, such as blue stragglers, x-ray binaries, radio pulsars, etc., of- 
ten referred to collectively as stellar exotica, or lusus natural In the disk of the Milky Way, 
such objects form through internal evolutionary processes in individual stars or close binary 
systems. In star clusters, these processes are augmented by stellar interactions, mediated by 
the high encounter frequency in dense cluster cores. 

Stellar encounters generate new channels for the formation of exotic objects, but can also 
catalyze existing channels. For example, a binary encounter may lead directly to a collision 
and the formation of a blue straggler, or its effect may be indirect, perhaps resulting in 
an exchange that eventually (billions of years later) leads to the formation of a low-mass 
x-ray binary. Repeated encounters can transform bin aries and multiple stellar systems 



multiplying the channels for the production of exotica ( IDavies 19951 . iHurley fe Shara 2002 



f Lusus Naturae is Latin for the freaks of nature, mutants or monsters. 
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Davies et al. 20061 ) . A clear understanding of the formation and ev olution of the se objects 
can provide insight into the past dynamical evolution of the cluster (IDavies 20091 ) . 

Many of the stellar exotica observed in old globular clusters today are the results of 
processes that began when the cluster was young. In some cases, they are the products 
of the interplay between dynamical and evolutionary processes involving stars and binaries 
during the first ~ 100 Myr of the clus ter lifetime. The primordial seeds for man y lusus 
naturae were planted during this period (jPortegies Zwart. McMillan fc Making 20071 ). Later 
stages (during phase 3, see §E} of dynamical evolution, such as core collapse, may produce 
additional generations of exotica. Since the observed YMCs will age to become old globular 
clusters (see they provide a convenient testbed for the study of the progenitors of stellar 
exotica. 

The progenitors of observed lusus naturae in old globular clusters are not necessarily easy 
to identify in the young cluster population, alth ough in some cases the evolutionary link is 
well established (IGlebbeek. Pols fc Hurley 20081 ). There may well be entire populations of 
peculiar objects in young star clusters that do not lead to observable interesting objects at 
later stages, and some objects destined for peculiarity may look perfectly ordinary at early 
times. An example of the latter is the dormant blue straggler population consisting of stars 
that were rejuvenated by mass transfer or collisions while still on the main sequence, and 
now lurk among their fellow main-sequence stars until they remain b ehind after the others 
traverse the Hertzsprung gap (jPortegies Zwart. Hut fc Verbunt 19971 ) . 



5.1 Binary Stars 

Exotic objects in star clusters are closely related to binaries, as they often form via internal 
binary evolution or during dynamical interactions between binaries and other stars. Ex- 
amples are the formation of blue stragglers ( §5.1.ip . colliding wind binaries ( §5.1.2)) . and 
anomalous x-ray pulsars ( §5.2.1)) . all of which require the presence of binary stars in the 
system. In some cases, such as the slow evolution of an accreting X-ray pulsar that leads 
from a low-mas s x-ray binary to a binary millisecond pulsar, the evolutionary track is read- 
ily established (IBhattacharya fc van den Heuvel 19911 ). In others, however, the punctuated 
equilibrium through which these objects evolve makes it virtually impossible to catch the key 
transitions as they occur. Examples are collisions between stars, or the common-envelope 
phase in the Darwin-Riemann instability of a contact binary. 

We will distinguish two fundamental types of binaries in star clusters: (1) "primordial" 
binaries, which formed conte mporaneously with the stars in the cluster as a crucial part 
of the star- formation process (IGoodman fe Hut 19891 ). and (2) "dynamical" binaries, which 
formed later via stellar interactions, often long after the component stars reached the main 
sequence. One can wonder to what extend this limited terminology is still usable for binaries 
that experiences one or more exchange interactions. It is for example, quite possible that two 
stars that were initially single end-up in a binary after two exchanges. As a practical matter 
we would still consider such a binary primordial. The second group of dynamical binaries 
may be further divided into two sub-categories — binaries formed by conservative three- and 
four-body stellar dynamical interactions, and binaries formed by dissipative two-body tidal 
capture. The latter process was introduced by lFabian. Pringle fc Rees (19751 ) to explain the 
relatively high specific frequency of low-mass X-ray binaries in old globular clusters. It has 
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fallen somewhat out of favor since the late 1980s, but it may be entering a revival of sorts 
Ogilvie fc Lin (2004 ). Many of the curious objects discussed in this § are related to binarity, 
either primordial or tidal, although we tentatively will use this distinction in the formation 
process. 



5.1.1 Blue Stragglers 

A blue straggler is a star which exceeds the cluster main-sequence turnoff in both temperature 
and luminosity, but which is not on the horizontal branch. Blue stragglers populate the 
region blueward of the turnoff, as if they lagged behind on the cluster main sequence while 
th e other stars age d. The first (34) blue stragglers were discovered in the globular cluster M3 
by Sandage (19531 ). A t least 8 plaus ible explanations have been proposed for the formation 
of blue stragglers (see Leonard 19891 ). Of these, two are currently in favor: 



direct merger between two stars (IHills fc Day 19761 ). 

mass transfer in a close semi-detached binary star ( McCrea 19641 ). 



The latter scenario is supported by the discovery of two blue stragglers, in the young open 
star clusters N GC 663 and NGC 6649, which have been found to be the donors in Be/X-ray 
binary systems (IMarco. Negueruela fc Motch 20071 ) . The discovery of a blue straggler in the 
old open cluster M67, which appears to be a bout 2.5 times more massive than the turn-off 
mass, favors the former view ([Leonard 19961 ). 

Both favored mechanisms for blue straggler formation appear plausible in YMCs. How- 
ever, no blue stragglers have yet been identified in any observed YMC, although this may 
be explained by the absence of a clear turn-off in the resolved clusters, which makes the 
identification (and the definition) of a blue straggler impractical. There are, however a 
number of "odd" stars in YMCs that might possibly evolve to resemble blue stragglers in 
the future. Objects consistent with this broadened definition inclu de four 03 If/WN6-A 



stars in the star cl uster R136 in the 30 Doradus region of the LMC ([Campbell et al. 1992 
Brandl et al. 1996h . 



Several interesting correlations exist betwee n the numbers of blue stragglers in old glob- 
ular clusters and the numbe rs of red giants (jFerraro. Fusi Pecci fc Bellazzini 19951 ). and 
also with the binary fraction (ISollima et al. 20081 ; Knigge. Leigh fc Sills 20091 ) . In addition, 
Davies. Piotto fc de Angeli (20041 ) find that the number of blue stragglers is be independent 
of integrated absolute magnitude My of the cluster, and use this fact to argue that both 
production mechanisms are relevant. 



5.1.2 Colliding wind binaries 

Binaries containing two massive stars, such as Wolf-Rayet stars, with strong stellar winds 
often exhibit intense radio and/or x-ray emission. Since this process requires copious stellar 
mass loss in a fast wind, these sources do not occur in old globular clusters, but YMCs 
appear to be excellent hosts for such systems. Several young and dense star clusters exhibit 
x-ray and radio emission from colliding w ind bi naries. In some re l atively nearby cases, R136 
(IPortegies Zwart. Poolev fc Lewin 2002 ). Wdl dClark et al. 20051 ICrowther et al. 2006h . and 
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the Arches and Quintuplet clusters (ILang et al. 20051 ). the counterparts of the radio and x- 
ray sources have been identified. 



5.2 Compact objects 



A compact object is generally produced in a supernova explosion (see § 15.3. ip . During such 
event the compact object, either a neutron star or a black hole, is likely to receive a high ve- 
locity, often referred to as a 'kick'. The direction of the kick is ill constrained, but the velocity 
distri bution of this kick is we l l determined by observing the velocities of nearby young neutron 
stars (ILyne fc Lorimer 1994 ; ICordes fc Chernoff 19981 ; lArzoumanian. Chernoff fc Cordes 20021) 



and the galactic scale -height of black hole x-ray transients (IWhite fc van Paradijs 1996 



Gualandris et al. 20051 ) . It seems that neutron stars receive a considerably higher velocity 
kick than black holes. As a result the majority (~ 95 %) of neutron st ars are likely ejected 



from any YMC (IDrukier 19961 ; iPfahl. Rappaport &: Podsiadlowski 20021 ) whereas most black 
holes are retained. However, since neutron stars are considerably more massive than the 
mean mass of a star in a cluster, once formed and retained they tend to stay in the cluster. 
At a later age compact objects can subsequently dominat e the cluster evolution. For black 



holes this happens at an age of about 100 Myr to a Gyr (IKulkarni. Hut fc McMillan 1993 



Mackev et al. 20071). whereas for neutron stars the age of domination starts around a Hubble 
time (IPortegies Zwart. McMillan fc Makino 20071 ) . The evolution of compact objects in star 
clusters deserves a separate review, but here we will present a very brief summary of the 
current observational evidence for compact objects in YMCs. 

A number of YMCs have been observed in the ra dio and x-ray, rey ealing a wealth of 
sources, even richer than found in old globular clusters (IClark et al. 20081 ). Among the x-ray 
sources is a large population of accreting neutron stars, stellar-mass black holes and possibly 
intermediate-mass black holes. Because of crowding in the central regions of these clusters, 
where most of the x-ray sources are found, very few sources have optical counterparts. 

With the adopted age limit of 100 Myr, only relatively few white dwarfs have formed — 
about as many as neutron stars — and cataclysmic variables are not expected. The youth 
of these clusters seem to make it unlikely that any low-mass x-ray binaries or millisecond 
pulsars will be found. 



5.2.1 Magnetars 

Shortly after a supernova (within ~ 10 5 years), a newly formed neutron star may become 
observable as a magnet a r, which can have a magnetic field strength exceeding ~ 10 15 gauss 



( Kouveliotou et al. 1999h . The population of magnetars is subdivided into two classes: soft 
gamma-ray repeaters (SGRs) and anomalous x-ray pulsars (AXPsjE As the products of 
supernovae, one might naively expect these objects to reside mainly in YMCs, and indeed 
about half (3 of 8) of the known SGRs and one-tenth (1 of 10) of the AXPs are known to 
reside in such systems. These are remarkably high fractions, given that only 0.05% of the 
stellar mass of the Galaxy resides in star clusters (see §T]). 



The single cluster AXP is CXOU J164710.2-455216 f lMuno et al. 20061 ) in the Galactic 
young star cluster Westerlund 1. It exhibits a 20-ms burst with energy ~ 10 37 erg (15-150 



=see |http : //www . physics . mcgill . ca/ ^pulsar /magnet ar/main . html 
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keV), and spins down at a rate of PjP ~ — 1CT 4 (jMuno et al. 20071 ). quite typical of a 
magnet ar. 

Several interesting sources are hosted by young star clusters in the LMC. These include 
the microquasar LS I +61° 303 , which may have been ejected from the cluster IC 1805 
([Mirabel. Rodrigues &: Liu 20041 ). The relatively low-density young star c luster SL 463, too 
small for inclusion in this review, seems to be associated with SGR 0526-66 (IKlose et al. 20041 ) 
at a projected distance of ~ 30 pc. Two other SGRs associated with relatively low-density 
young star clusters a re SGR 1806-20, at a projected distance of ~ 0.4 pc fro m the core 
of its parent cluster (IMirabel. Fuchs & Chatv 2000l: ICorbel k Eikenberrv 20o4) . and SGR 
1900+14, at ~ 0.8 pc (IMirabel. Fuchs fc Chatv 2000l : IVrba et al. 2000h . This latter SGR 
has a measured proper motion of 70 mas/yr away from the cluster, suggesting that it indeed 
was ejected from the parent cluster ( DeLuca et al. 2009h. which is relatively old (14 ± 1 Myr) 
compared to the usual SGR-producing stars (jDavies et al. 20091 ). The actual association 
between cluster and SGR is hard to establish in this case since the distances to both objects 
are ill constrained. 



5.2.2 Ultra-luminous X-ray sources 



Old globular clusters are known to host an enormous excess of low-mass x-ray binaries 
compared to the rest of the Gala xy. Much of this excess is attributed to the dy namical 
environment in dense cluster cores (IFabian. Pringle fc Rees 19751 ; iPooley et al. 20031 ) . Young 
star clusters are sites of intense dynamical activity, so it is not surprising that YMCs also 
host many x-ray sources. The majority of x-ray point sources in external galaxies appear to 
be associated with young star clu sters, as is the case for example in the Antennae system 
(NGC4038/39) IZezas et al. (20021 ) . The nature of most of these x-ray sources is unknown, 
and we can only guess at their origin. 

We limit ourselves here to the most striking x-ray sources, the subclass of ultra-luminous 
x-ray sources (ULXs), which are characterized by x-ray luminosity L x ^ 1.3 x 10 39 erg/s, 
the maximum isotr opic luminosity that can be produced by a 10 M Q black hole accret- 
ing pure hydrogen (IKing et al. 20011 ) . For practical reasons we round the threshold up to 
L x ^ 10 40 erg/s, mainly to ensure that such luminosities are unlikely to be produced by 
stellar-mass black holes accreting at the Eddington rate from a main-sequence companion 
star. These ULXs are responsible for the brightest stellar x-ray sources in the sky. 

Several models attempt to explain the high x-ray luminosity of the ULXs, but at present 
there is no consensus in the community on the source of the x-rays. The current leading 
models are: 

• anisotropic (collimated or beamed) emission from an accreting stellar-mass bla ck hole 



(although p orosity, turbulence, and bubbles also provide interesting alternatives. IKing et al. 2001 



King 20021 ) 



accretion from an evolved star onto a stellar-mass black hole, which can, in principle, 
lead to an accreti on rate higher than from a main-sequence star, and therefore a higher 
x-ray luminosity (IMadhusudhan et al. 20061 ). 



accretion from a c ompanion star onto an "inte rmediate-mass black hole" (IMBH), with 



mass £ 



100 M (IPorteries Zwart et al. 2004T ). 
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ULXs tend to be hosted by starburst and spiral galaxies (IMakishima et al. 20001 ) . Some 
of the brightest are associated with YM Cs; a leading examp le is the ULX in the star cluster 
MGG 11 in the starburst galaxy M82 (IKaaret et al. 20011) . The association with YMCs 
argues in favor of an accreting black hole of ~ 1000 M (IPortegies Zwart et al. 2004 ). The 
object in MGG 11 is particularly inte resting, as it shows a strong qua si-periodic oscillation 
in the 50-100 mHz frequency range (jStrohmayer &: Mushotzky 20031 ). providing a strong 
argument against beamed emission, and supporting the hypothesis that the x-ray luminosity 
comes from an accreting black hole of 200-5000 M . Further support is provided by the 
detected periodic variation of ~ 62 days, whi ch can be explained if the black hole is orbited 
by a 22-25 M Q Roche-lobe filling donor star (IPatruno et al. 20061). 

ULXs have been asso ciated with YMCs i n NGC 5204 dLiu. Bregman fc Seitzer 2004T) . the 
starburst galaxies M82 flKaaret et al. 200lh and NGC1313 flGrise et al. 20081) the edge-on 
spiral NGC 4565 JWu et al. 2002]), the intera cting galaxies M51 (lLmet_aL2fi02|) , NGC4038 /39 
(Antennae iFabbiano. Zezas fc Murray 200 l[) and ESQ 350-40 (Cartwheel iGao et al. 20031) . 



and the type 1.5 Seyfert galaxy NGC 1275 ( IGonzalez-Martm. Fabian &: Sanders 20061 ). The 
higher abundance of ULXs in active, starburst, and interacting galaxies may be related to the 
empirical fact that YMCs tend t o form in these environments; 6 0% of ULXs are associated 
with active star-forming regions ( jSwartz. Tennant &: Soria 20091 . although the definition of 
a ULX used here is somewhat faint). 

If the counterpart of a ULX hosts an IMBHs, it is likely to be the the acceptor from a 
windy or Roche-lobe overflowing massive star, as seems to be the case in the 10 39 -10 41 erg/s 
ULX in NGC 5204, where the donor is ide ntified as a B0 lb supergiant with a 10-day 
orbital period ( iLiu. Bregman &: Seitzer 20041) . and in ULX M51 X-7, which has has an even 
shorter orbital period of only 2.1 hr ( ILiu et al. 20021 ). although no stellar companion has been 
identified. The black holes that may be responsible for the ULXs NGC 1313 X-l and X-2 
(0.2-10.0 keV, assu ming a distance of 3.7 Mpc) may have masses in the range 100-1000 M 
( IMiller et al. 20031 ). and in these cases YMCs have been identified as optical counterparts. 



5.3 Explosive events 
5.3.1 Supernovae 

Supernova a re relatively rare events, occurring a bout once every 100 years in a galaxy like the 
Milky Way ( Cappellaro. Evans fc Turatto 19991). Although type I supernovae are u nlikely to 
occur in star clusters younger than 100 Myr ( IPfahl. Scannapieco &: Bildsten 20091). at these 
ages the seeds may be planted for a rich future of type I events (jShara &: Hurley 20021 ). 

Since most massive stars tend to reside in clusters, as discussed in JTJand §0 it is probable 
that the majority of type Ib/c and type II supernovae occur in star clusters. However, 
since most supernovae occur in distant galaxies it is hard to find optical counterparts, and 
very few associations of supernovae wit h young star clu sters have been reported. Based 
on the peculiar metallicity of SN 1987A, lEfremov ( 19911) argue that the star originated in 
the young LMC cluster MKM90. SN 2006gv flFolev et al. 2006h . which is a candidate for 
the fo rmation of an IMBH, may have been th e result of a collision runaway ( §3.4.2j) in a 
YMC ( IPortegies Zwart fc van den Heuvel 2007b . However, perhaps the strongest case is the 
peculiar type Hp supernova SN 2004dj (probably produced by a 12-20 M star) in the spiral 
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galaxy NGC 2403; as it faded, the star cluster Sandage-96 reappeared (jWang et al. 20051 ). 



5.3.2 Gamma-ray bursts 



Any word written about gamma-ray bursts is likely to trigger its own burst of e-mail, but we 
cannot resist the temptation to devote a few lines to this fascinating transient phenomenon. 
Gamma-ray bu rsts come in two types, of short ( £ 2 s) and long duration, respectively 
Meszaros (20021 ) . Several theories exist which point to either old globular or young star 
clusters as possible hosts for ga mma-ray bursts. Colliding com pact objects are often cited 
as sources for the short bursts (INarayan. Piran &: Shemi 19911) . Long bursts are thought 
to be hosted by massive star forming regions ( Paczynski 19981) . Of pa rticular interest is 



the elusive relation between the long bursts and YMCs (jEfremov 20001 ) . The models for 
long-duration gamma-ray bursts should be particu larly applicable to YMCs, as they require 
rapidly rotating high- mass stars (jHeger et al. 20031 ). which could be achieved quite naturally 
by stellar collisions in a YMC (see §3.4.2j) . 



5.4 Summary of exotica 

We have discussed several examples of exotica in YMCs, but other curiosities remain. Many 
of these exotic objects are well studied in old globular clusters, but similar scrutiny is so 
far lacking in their younger siblings. Rather than providing a detailed description of each 
of the oddities found, we simply mention a few developments and recent discoveries of what 
today we call exotic objects, which one day we may consider "normal." The following list 
summarizes a number of peculiar clusters, noteworthy because of a unique source or object. 
The list is far from complete, but at the very least it indicates the diversity of objects found 
in YMCs. 



R136 Contains some 13 colliding wind binaries (jPortegies Zwart. Pooley fc Lewin 20021 ) 
and possibly 3 blue stragglers (even though no clear turn-off can be distinguished). In 
addition, the star cluster shows an OH (1720 MHz) Mase r, which is probably related to 
the surrounding nebula rather than the star cluster itself ( [Roberts fc Yusef-Zadeh 20051 ). 



The Arches cluster contains 10 radio point sources (ILang et al. 20051 ) and several 
colliding wind binaries. 



MGG11 is a YMC in M82 which may contain a ULX (IKaaret et al. 20011 ). alhtough 
the Chandra error box is slightly offset. 



The Quintuplet cluster contains the Pistol star (jFiger et al. 19981 ). a candidate for 
the most massive star in the Galaxy, a t a projected distance of about 1 pc, as well as 
9 radio point sources ( ILang et al. 20051 ). 



Westerlund 1 hosts the anom alous x-ray pulsar CXOU J164710.2-455216 flMuno et al. 20061 ). 
and a wealth of x-ray sources (jClark et al. 20081 ). 



Westerlund 2 hosts the massive Wolf-Raye t binary WR 20a, containing two WN6ha 
stars, at a distance of 1.1 pc from its center (IRauw et al. 20051 ) . 
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NGC 663 and NGC 664 9 contain blue stragglers which f ound to be the donors of 
Be/X-ray binary systems (IMarco. Negueruela fc Motch 20071 ). 



Sandage-96 is a ~ 96, 000 M m star cluster in NGC 2403 in which a type Hp super- 
nova was det ected (Wang et al. 20051 ). This star cluster also exhibits multiple stellar 
populations f Vinko et al. 2009^ 



5.4.1 Planetary nebulae and supernova remnants 

The nuclear evolution of a sufficiently massive star ( £ 8 M ) is associated with a supernova 
explosion (see §5.3. f p . while lower mass stars fizzle into the background after a short, bright 
post-AGB phase. But following these lower mass events a roughly spherical gas shell — a 
planetary nebula — remains visible for a much longer time than either the supernova or the 
post-AGB star, illuminated by the central stellar remnant and shocks as the out-flowing gas 
encounters the interstellar medium. 

Every star experiences either a supernova or a post AGB phase and a cluster of 5 x 10 4 
stars will experience some 350 supernovae, leading to an equal number of remnants, while 
during the first 100 Myr a similar number of planetary nebulae will form. The formation 
rate of planetary nebulae and supernova remnants is thus ~ 7/Myr. With an observable 
lifetime for a nebula of about 10 yr, we naively expe ct to see one nebula for ev ery ~ 14 star 
clusters. Among 80 YMCs in the Milky Way Galaxy lLarsen &: Richtler (20061 ) found 3 with 
a planetary nebula, consistent with our naive estimate. No planetary nebulae or supernova 
remnants have so far been found in any of the YMCs in the Local Group. But since only 16 
of the clusters listed in Tab. are sufficiently old to have formed white dwarfs from single 
stellar evolution, the lack of planetary nebulae may be a statistical fluctuation. 



5.4.2 Brown dwarfs and planets 

We will say little here about brown dwarfs and planets in YMCs since the observations are 
sparse and the general topic deserves its own review. However, a few words are in order. 

Young dense star clusters are generally too distant fo r planets to be detect able using cur- 
rent methods, and no planets have been found to da te (lUdrv fc Santos 20071) . Old globular 
clusters seem to be defici ent in planets (e.g. 47 Tuc, iGilliland et al. 20001 ). possibly because 
of their low metallicities (IWeldrake et al. 20051). How ever, the recent HARPS discovery of a 
planet around a metal poor star (IMayor et al. 20091 ) makes this argument less convincing. 
The planet found or biting the 11 ms pulsar B 1620-26 in the metal poor e nvironment of the 
globular cluster M4 ( Backer 1993 ; Thorsett. Arzoumanian fc Taylor 1993| ). was formed by a 
different mechanism than planets around solar-type stars, and we do not (yet) expect such 
planets in YMCs. Star clusters in high-metallicity environments, such as Westerlund 1 and 
NGC 3603, could host a large population of planets, but none have yet been found. 

We see no reason why YMCs should be deficient in planets. However, the high interaction 
rate in a dense cluster could make a planetary system short-li ved. Disruption of a planetary 
system may leave the planet separated from its parent star ( Spurzem et al. 2009h. Several 
such free-floating objects have been found in the Orion Trapezium cluster ( ILucas fc Roche 20001 ) 
Once dislocated from its parent star, a planet will easily escape the cluster, though. 
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5.4.3 YMCs in a galactic context 



Due to supernovae and winds from massive stars , YMCs are the sources of cluster winds 
(jSifich. Tenorio-Tagle &: Rodriguez- Gonzalez 20041 ). whic h may trigger larger-scale winds and 
chimneys, as in the Perseus arm of the Milky Way Galaxy ( INormandeau. Taylor &: Dewdney 19961 ). 
Although little studied from the point of view of cluster evolution, this is an important topic 
that provides a possible and rather natural link between the evolution of a YMC and that 
of its parent galaxy. 



6 Concluding Remarks: Young Globular Clusters? 

The discovery of large numbers of young massive star clusters, particularly in other galaxies, 
over the last decade has led to the realization that such clusters are responsible for a signifi- 
cant fraction of all current star formation in the local universe. The study of these systems, 
and especially their lifetimes, against various stellar evolutionary and dynamical processes, 
is therefore of critical importance to several branches of stellar and galactic astrophysics. 

Star clusters appear to form with a cluster mass function described by a power-law with 
index —2. This mass function seems to be the same for both open and globular clusters, and 
does not depend significantly on the local galactic environment or on the specific character- 
istics of the giant molecular clouds from which the clusters formed. An important parameter 
for a young bound cluster appears to be its age relative to its current dynamical time scale 
tdyn- For unbound stellar agglomerates or associations, t& yri exceeds the system age, indicat- 
ing that the cluster is either extremely young or dissolving into the tidal field of its parent 
galaxy. For the typical bound star clusters discussed in this review, t& jn is smaller than the 
current age; for an association it is the other way around. 

YMCs evolve and eventually dissolve due to the combined effects of a number of physical 
processes, the most important (for clusters that survive the early expulsion of their natal 
gas) being mass loss due to stellar evolution. The most massive clusters, such as those found 
in the Antennae system, have expected lifetimes comparable to the age of the universe, and 
we could well imagine that the Antennae will someday be a medium-sized elliptical galaxy 
with an extended population of intermediate-age clusters having overall properties quite 
comparable to the old globular clusters seen in other ellipticals. 

The seeds of the exotica observed in many present-day globular clusters were initiated 
during the infancy of those systems, in the strongly coupled mix of stellar dynamics and 
stellar evolution that characterized their early evolution (in particular dring phase 2). YMCs 
destined to survive to ages comparable to the globular clusters appear to contain much richer 
populations of stellar exotica (per unit mass) than are found in the field, and may provide 
important testbeds for this unique period in cluster evolution. Our limiting cluster age of 
100 Myr is chosen in part to include the period when this ecological interplay is strongest. 

From an observational point of view, little is known of the formation and evolution of 
stellar exotica in YMCs, mainly because such clusters are relatively rare. The Milky Way 
contains only half a dozen, and the closest lies ~ 4kpc away, too distant for detailed study 
of the individual stars in its central region. However, a number of studies have drawn 
connections between YMCs and exotic objects, such as unusual supernovae, magenetars, 
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x-ray binaries, and ultraluminous x-ray sources, possibly placing YMCs in the same league 
as the old GCs, which are rich in such objects. 

Finally, we return to the term "young globular cluster," which we introduced in £Q] but 
have deliberately avoided throughout this review. Is this an appropriate term for a YMC? 
A fairly common definition of a globular cluster, similar to that found in many textbooks, 
appears in the Oxford English Dictionary (2009): 11 A roughly spherical cluster of stars, 
typically seen in galactic halos, containing large numbers of old, metal-poor stars." This 
definition contains several qualifiers that would seem to exclude YMCs, but in many ways 
the definition itself seems to be a relic of a bygone era. 

Our own Galaxy contains a significant population of bulge GCs, and there may well 
be a disk component among GCs in other galaxies, such as the LMC. In any case, it is 
quite conceivable that YMCs will in many cases come to populate the halos of their parent 
galaxies. Since there appears to be little dependence of YMC properties on current galactic 
location, it seems unreasonable to exclude a YMC from the GC definition based solely on 
this property. The "low metallicity" qualifier in the definition also seems a side effect rather 
than a requirement. If we required simply that globular clusters be old — or rather long-lived, 
as in our YMC definition — the metallicity becomes redundant, merely reflecting the epoch 
at which they formed. Furthermore, the LMC GCs display a broad range of ages, so any 
cutoff on age or metallicity would be arbitrary. Lastly, any massive cluster more than a few 
tens of dynamical times old will necessarily have a smooth, roughly spherical appearance, 
regardless of metallicity or location. 

Applying all of these arguments, we can strip away all of the extraneous qualifying terms 
in the above GC definition until all that remains is the word "massive." In this view, 
globular clusters are simply "old massive clusters," the logical descendants of YMCs in the 
early universe. 

Acknowledgments 

We are particularly grateful to Douglas Heggie, S0ren Larsen, Nate Bastian and Fred Ra- 
sio for their careful reading of an earlier version of this manuscript and their useful (and 
often quite critical) comments. In addition we received very helpful comments from Holger 
Baumgardt, Melvyn Davies, Richard de Grijs, Diederik Kruijssen and Hans Zinnecker. The 
results of the iV-body simulation shown in Fig. [T41 was done on one of the GRAPE-6 BLX64 
boards of the European Southern Observatory in Garching. This work was made possible 
with the financial support of DFG cluster of excellence Origin and Structure of the Universe 
(www.universe-cluster.de), NASA (grant NNX07AG95G), the National Science Foundation 
(NSF, via grant AST-0708299), Nederlandse Onderzoekschool voor Astronomie (NOVA), het 
Leids-Kerkhoven Bosscha Fonds (LKBF), the Kavli Institute for Theoretical Physics (KITP, 
through their hosptality and NSF grant PHY05-51164), and in particular of the Neder- 
landse Organizatie voor Wetenschappelijk Onderzoek (NWO, via grants #643.200.503 and 
#639.073.803) and the Nederlandse Computer Faciliteiten (NCF, via project #SH-095-08). 



61 



References 

Aarseth SJ. 2003. Gravitational N-Body Simulations. Gravitational N-Body Simulations, 
by Sverre J. Aarseth, pp. 430. ISBN 0521432723. Cambridge, UK: Cambridge University 
Press, November 2003. 

Abel T, Bryan GL, Norman ML. 2002. Science 295:93-98 

Abt HA. 1983. ARA&A 21:343-372 

Ahmad A, Cohen L. 1973a. J. of Comp. Phys. 12:289 

Ahmad A, Cohen L. 1973b. ApJ 179:885-896 

Allison RJ, Goodwin SP, Parker RJ, de Grijs R, Portegies Zwart SF, Kouwenhoven MBN. 
2009. ApJL 700:L99-L103 

Ambartsumian VA. 1938. Sci. Mem. Leningrade State Univ. #22, ser. Math. Sci. (astron- 
omy) 4:19 

Anders P, de Grijs R, Fritze-v. Alvensleben U, Bissantz N. 2004. MNRAS 347:17-28 

Andersen M, Zinnecker H, Moneti A, McCaughrean MJ, Brandl B, et al. 2009. ApJ 
707:1347-1360 

Antonov VA. 1962. Solution of the problem of stability of stellar system Emden's density law 
and the spherical distribution of velocities 

Applegate JH. 1986. ApJ 301:132-144 

Arp H, Sandage A. 1985. A J 90:1163-1171 

Arzoumanian Z, Chernoff DF, Cordes JM. 2002. ApJ 568:289-301 
Ascenso J, Alves J, Lago MTVT. 2009. A&A 495:147-155 
Ascenso J, Alves J, Vicente S, Lago MTVT. 2007. A&A 476:199-215 
Ashman KM, Zepf SE. 2001. A J 122:1888-1895 

Backer DC. 1993. In Planets Around Pulsars, ed. J. A. Phillips, S. E. Thorsett, & S. R. Kulka- 
rni, vol. 36 of Astronomical Society of the Pacific Conference Series 

Bacon D, Sigurdsson S, Davies MB. 1996. MNRAS 281:830-846 

Bahcall JN, Wolf RA. 1976. ApJ 209:214-232 

Barmby P, Perina S, Bellazzini M, Cohen JG, Hodge PW, et al. 2009. A J 138:1667-1680 
Barnes J, Hut P. 1986. Nat 324:446-449 
Bastian N. 2008. MNRAS 390:759-768 

62 



Bastian N, Covey KR, Meyer MR. 2010. ARA&A , in press karXiv.lOOlWJMl 
Bastian N, de Mink SE. 2009. MNRAS 398:L11-L15 

Bastian N, Gieles M, Goodwin SP, Trancho G, Smith LJ, et al. 2008. MNRAS 389:223-230 

Bastian N, Gieles M, Lamers HJGLM, Scheepmaker RA, de Grijs R. 2005. A&A 431:905- 
924 

Bastian N, Konstantopoulos I, Smith LJ, Trancho G, Westmoquette MS, Gallagher JS. 2007. 
MNRAS 379:1333-1342 

Bastian N, Saglia RP, Goudfrooij P, Kissler-Patig M, Maraston C, et al. 2006. A&A 448:881- 
891 

Bastian N, Trancho G, Konstantopoulos IS, Miller BW. 2009. ApJ 701:607-619 

Bate MR, Bonnell IA, Bromm V. 2003. MNRAS 339:577-599 

Battinelli P, Capuzzo-Dolcetta R. 1991. MNRAS 249:76-83 

Baumgardt H. 1998. A&A 330:480-491 

Baumgardt H. 2001. MNRAS 325:1323-1331 

Baumgardt H, De Marchi G, Kroupa P. 2008. ApJ 685:247-253 

Baumgardt H, Kroupa P. 2007. MNRAS 380:1589-1598 

Baumgardt H, Kroupa P, Parmentier G. 2008. MNRAS 384:1231-1241 

Baumgardt H, Makino J. 2003. MNRAS 340:227-246 

Baumgardt H, Makino J, Hut P. 2005. ApJ 620:238-243 

Bell EF, Zucker DB, Belokurov V, Sharma S, Johnston KV, et al. 2008. ApJ 680:295-311 

Bellazzini M, Fusi Pecci F, Messineo M, Monaco L, Rood RT. 2002. A J 123:1509-1527 

Belleman RG, Bedorf J, Portegies Zwart SF. 2008. New Astronomy 13:103-112 

Benz W, Hills JG. 1987. ApJ 323:614-628 

Bettwieser E, Sugimoto D. 1984. MNRAS 208:493-509 

Bhattacharya D, van den Heuvel EPJ. 1991. PhysRep 203:1-124 

Bik A, Lamers HJGLM, Bastian N, Panagia N, Romaniello M. 2003. A&A 397:473-486 

Binney J, Tremaine S. 2008. Galactic Dynamics: Second Edition. Princeton University Press 

Bisseling RH. 2004. Parallel Scientific Computation: A Structured Approach using BSP and 
MPI. Oxford University Press, pp. 324. ISBN 0-19-852939-2, March 2004. 

63 



Bonnell IA, Bate MR. 2006. MNRAS 370:488-494 

Bonnell IA, Bate MR, Vine SG. 2003. MNRAS 343:413-418 

Bosch G, Terlevich E, Terlevich R. 2009. AJ 137:3437-3441 

Boutloukos SG, Lamers HJGLM. 2003. MNRAS 338:717-732 

Brandl B, Sams BJ, Bertoldi F, Eckart A, Genzel R, et al. 1996. ApJ 466:254 

Brandner W, Clark JS, Stolte A, Waters R, Negueruela I, Goodwin SP. 2008. A&A 478:137- 
149 

Bruzual G, Chariot S. 2003. MNRAS 344:1000-1028 

Caldwell N, Harding P, Morrison H, Rose JA, Schiavon R, Kriessler J. 2009. A J 137:94-110 

Campbell B, Hunter DA, Holtzman JA, Lauer TR, Shayer EJ, et al. 1992. AJ 104:1721- 
1742 

Campbell MA, Evans CJ, Mackey AD, Gieles M, Alves J, et al. 2010. MNRAS , in press 
karXiv:1002.0288\) 

Cantiello M, Brocato E, Blakeslee JP. 2009. A&A 503:87-101 

Cappellaro E, Evans R, Turatto M. 1999. A&A 351:459-466 

Carretta E, Bragaglia A, Gratton RG, Lucatello S. 2008. ArXiv e-prints 

Casertano S, Hut P. 1985. ApJ 298:80-94 

Chabrier G. 2003. PASP 115:763-795 

Chandar R, Bianchi L, Ford HC. 2000. AJ 120:3088-3097 

Chandar R, Bianchi L, Ford HC, Salasnich B. 1999. PASP 111:794-800 

Chandar R, Fall SM, Whitmore BC. 2006. ApJL 650:L111-L114 

Chatterjee S, Fregeau JM, Rasio FA. 2008. In IAU Symposium, ed. E. Vesperini, M. Giersz, 
& A. Sills, vol. 246 of IAU Symposium 

Chernoff DF, Weinberg MD. 1990. ApJ 351:121-156 

Chiosi E, Vallenari A, Held EV, Rizzi L, Moretti A. 2006. A&A 452:179-193 

Clark JS, Muno MP, Negueruela I, Dougherty SM, Crowther PA, et al. 2008. A&A 477:147- 
163 

Clark JS, Negueruela I, Crowther PA, Goodwin SP. 2005. A&A 434:949-969 

Clark JS, Negueruela I, Davies B, Larionov VM, Ritchie BW, et al. 2009. A&A 498:109-114 

64 



Clarke CJ, Bonnell IA, Hillenbrand LA. 2000. Protostars and Planets IV :151 

Colin H. 1979. ApJ 234:1036-1053 

Colin H. 1980. ApJ 242:765-771 

Corbel S, Eikenberry SS. 2004. A&A 419:191-201 

Cordes JM, Chernoff DF. 1998. ApJ 505:315-338 

Crowther PA, Hadfield LJ, Clark JS, Negueruela I, Vacca WD. 2006. MNRAS 372:1407- 
1424 

Davies B. 2009. ArXiv e-prints 

Davies B, Figer DF, Kudritzki R, Trombley C, Kouveliotou C, Wachter S. 2009. ArXiv 
e-prints 

Davies B, Figer DF, Kudritzki RP, MacKenty J, Najarro F, Herrero A. 2007. ApJ 671:781- 
801 

Davies B, Figer DF, Law CJ, Kudritzki R, Najarro F, et al. 2008. ApJ 676:1016-1028 
Davies MB. 1995. MNRAS 276:887-905 

Davies MB, Amaro-Seoane P, Bassa C, Dale J, de Angeli F, et al. 2006. New Astronomy 
12:201-214 

Davies MB, Piotto G, de Angeli F. 2004. MNRAS 349:129-134 
de Grijs R, Anders P. 2006. MNRAS 366:295-307 

de Grijs R, Fritze-v. Alvensleben U, Anders P, Gallagher JS, Bastian N, et al. 2003. MNRAS 
342:259-273 

de Grijs R, Gilmore GF, Johnson RA, Mackey AD. 2002a. MNRAS 331:245-258 

de Grijs R, Gilmore GF, Mackey AD, Wilkinson MI, Beaulieu SF, et al. 2002b. MNRAS 
337:597-608 

de Grijs R, Goodwin SP. 2008. MNRAS 383:1000-1006 

de Grijs R, Parmentier G. 2007. Chinese Journal of Astronomy and Astrophysics 7:155-186 
de Wit WJ, Testi L, Palla F, Zinnecker H. 2005. A&A 437:247-255 

Decressin T, Meynet G, Charbonnel C, Prantzos N, Ekstrom S. 2007. A&A 464:1029-1044 
Dehnen W. 2000. ApJL 536:L39-L42 

Deiters S, Spurzem R. 2001. Astronomical and Astrophysical Transactions 20:47-50 



65 



DeLuca A, Caraveo PA, Esposito P, Hurley K. 2009. ApJ 692:158-161 

D'Ercole A, Vesperini E, D'Antona F, McMillan SLW, Recchi S. 2008. MNRAS 391:825-843 

Dias WS, Alessi BS, Moitinho A, Lepine JRD. 2002. A&A 389:871-873 

Djorgovski SG, Meylan G, eds. 1993. Structure and dynamics of globular clusters, vol. 50 of 
Astronomical Society of the Pacific Conference Series 

Dorband EN, Hemsendorf M, Merritt D. 2003. Journal of Computational Physics 185:484- 
511 

Drukier GA. 1996. MNRAS 280:498-514 

Drukier GA, Fahlman GG, Richer HB. 1992. ApJ 386:106-119 
Duquennoy A, Mayor M. 1991. A&A 248:485-524 

Ebisuzaki T, Makino J, Fukushige T, Taiji M, Sugimoto D, et al. 1993. Publ. Astr. Soc. 
Japan 45:269-278 

Efremov YN. 1991. Pis ma Astronomicheskii Zhurnal 17:404-409 

Efremov YN. 2000. Astronomy Letters 26:558-564 

Eggleton P. 2006. Evolutionary Processes in Binary and Multiple Stars 

Eggleton PP, Fitchett MJ, Tout CA. 1989. ApJ 347:998-1011 

Einsel C, Spurzem R. 1999. MNRAS 302:81-95 

Elmegreen BG. 2007. ApJ 668:1064-1082 

Elmegreen BG, Efremov YN. 1997. ApJ 480:235 

Elson RAW, Fall SM. 1985a. PASP 97:692-696 

Elson RAW, Fall SM. 1985b. ApJ 299:211-218 

Elson RAW, Fall SM, Freeman KG 1987. ApJ 323:54-78 

Elson RAW, Freeman KC, Lauer TR. 1989. ApJL 347:L69-L71 

Espinoza P, Selman FJ, Melnick J. 2009. ArXiv e-prints 

Fabbiano G, Zezas A, Murray SS. 2001. ApJ 554:1035-1043 

Fabian AC, Pringle JE, Rees MJ. 1975. MNRAS 172: 15P 

Fall SM, Chandar R, Whitmore BC. 2005. ApJL 631:L133-L136 

Fall SM, Rees MJ. 1985. ApJ 298:18-26 



66 



Fall SM, Zhang Q. 2001. ApJ 561:751-765 

Fellhauer M, Wilkinson MI, Kroupa P. 2009. MNRAS 397:954-962 
Ferraro FR, Fusi Pecci F, Bellazzini M. 1995. A&A 294:80-88 

Figer DF, MacKenty JW, Robberto M, Smith K, Najarro F, et al. 2006. ApJ 643:1166-1179 
Figer DF, McLean IS, Morris M. 1999. ApJ 514:202-220 

Figer DF, Najarro F, Gilmore D, Morris M, Kim SS, et al. 2002. ApJ 581:258-275 

Figer DF, Najarro F, Morris M, McLean IS, Geballe TR, et al. 1998. ApJ 506:384-404 

Fleck J, Boily C, Lancon A, Heggie D, Deiters S. 2005. In SF2A-2005: Semame de 
VAstrophysique Francaise, ed. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani 

Fleck J J, Boily CM, Lancon A, Deiters S. 2006. MNRAS 369:1392-1406 

Foley RJ, Li W, Moore M, Wong DS, Pooley D, Filippenko AV. 2006. Central Bureau 
Electronic Telegrams 695:1 

Forte JC, Vega EI, Faifer F. 2009. MNRAS 397:1003-1020 

Foster I, Kesselman C, eds. 2004. The Grid: Blueprint for a New Computing Infrastructure. 
Published by Morgan Kaufmann 

Fregeau JM, Cheung P, Portegies Zwart SF, Rasio FA. 2004. MNRAS 352:1-19 

Fregeau JM, Giirkan MA, Joshi KJ, Rasio FA. 2003. ApJ 593:772-787 

Fregeau JM, Rasio FA. 2007. ApJ 658:1047-1061 

Freitag M, Benz W. 2001. A&A 375:711-738 

Freitag M, Benz W. 2005. MNRAS 358:1133-1158 

Freitag M, Rasio FA, Baumgardt H. 2006. MNRAS 368:121-140 

Fujii M, Iwasawa M, Funato Y, Makino J. 2007. Publ. Astr. Soc. Japan 59:1095- 

Fukushige T, Heggie DC. 1995. MNRAS 276:206-218 

Fukushige T, Heggie DC. 2000. MNRAS 318:753-761 

Gaburov E, Gieles M. 2008. MNRAS 391:190-196 

Gaburov E, Harfst S, Portegies Zwart SF. 2009. New Astronomy 14:630-637 
Gaburov E, Lombardi JC, Portegies Zwart S. 2008. MNRAS 383:L5-L9 
Gao B, Goodman J, Cohn H, Murphy B. 1991. ApJ 370:567-582 

67 



Gao Y, Wang QD, Appleton PN, Lucas RA. 2003. ApJL 596:L171-L174 
Geyer MP, Burkert A. 2001. MNRAS 323:988-994 
Gieles M. 2009. MNRAS 394:2113 

Gieles M, Athanassoula E, Portegies Zwart SF. 2007. MNRAS 376:809-819 

Gieles M, Baumgardt H. 2008. MNRAS 389:L28-L32 

Gieles M, Lamers HJGLM, Portegies Zwart SF. 2007. ApJ 668:268-274 

Gieles M, Larsen SS, Bastian N, Stein IT. 2006a. A&A 450:129-145 

Gieles M, Larsen SS, Scheepmaker RA, Bastian N, Haas MR, Lamers HJGLM. 2006b. A&A 
446:L9-L12 

Gieles M, Portegies Zwart SF, Baumgardt H, Athanassoula E, Lamers HJGLM, et al. 2006c. 
MNRAS 371:793-804 

Gieles M, Sana H, Portegies Zwart SF. 2009. MNRAS, in press (jarXiv:0911.155T\ ) 

Giersz M. 1998. MNRAS 298:1239-1248 

Giersz M. 2001. MNRAS 324:218-230 

Giersz M. 2006. MNRAS 371:484-494 

Giersz M, Heggie DC. 1994. MNRAS 268:257 

Giersz M, Heggie DC. 2009. MNRAS 395:1173-1183 

Giersz M, Spurzem R. 2003. MNRAS 343:781-795 

Gies DR. 1987. ApJS 64:545-563 

Gilliland RL, Brown TM, Guhathakurta P, Sarajedini A, Milone EF, et al. 2000. ApJL 
545:L47-L51 

Girardi L, Chiosi C, Bertelli G, Bressan A. 1995. A&A 298:87 

Glebbeek E, Gaburov E, de Mink SE, Pols OR, Portegies Zwart SF. 2009. A&A 497:255-264 
Glebbeek E, Pols OR, Hurley JR. 2008. A&A 488:1007-1015 
Gnedin OY, Ostriker JP. 1997. ApJ 474:223 

Gonzalez-Martin O, Fabian AC, Sanders JS. 2006. MNRAS 367:1132-1138 
Goodman J, Hut P. 1989. Nat 339:40-42 
Goodwin SP. 2009. Ap&SS :108 



6s 



Goodwin SP, Bastian N. 2006. MNRAS 373:752-758 

Gouliermis D, Keller SC, Kontizas M, Kontizas E, Bellas- Velidis I. 2004. A&A 416:137-155 
Grise F, Pakull MW, Soria R, Motch C, Smith IA, et al. 2008. A&A 486:151-163 
Gualandris A, Colpi M, Portegies Zwart S, Possenti A. 2005. ApJ 618:845-851 
Giirkan MA, Freitag M, Rasio FA. 2004. ApJ 604:632-652 

Gutermuth RA, Megeath ST, Pipher JL, Williams JP, Allen LE, et al. 2005. ApJ 632:397- 
420 

Gvaramadze VV, Bomans DJ. 2008. A&A 490:1071-1077 
Hamada T, Iitaka T. 2007. ArXiv Astrophysics e-prints 
Harayama Y, Eisenhauer F, Martins F. 2008. ApJ 675:1319-1342 

Harfst S, Gualandris A, Merritt D, Spurzem R, Portegies Zwart S, Berczik P. 2007. New 
Astronomy 12:357-377 

Harfst S, Portegies Zwart S, Stolte A. 2009. ArXiv e-prints 

Harris J, Zaritsky D. 2009. A J 138:1243-1260 

Harris WE. 1996. AJ 112:1487 

Harris WE. 2001. In Saas-Fee Advanced Course 28: Star Clusters, ed. L. Labhardt & 
B. Binggeli 

Heger A, Fryer CL, Woosley SE, Langer N, Hartmann DH. 2003. ApJ 591:288-300 

Heggie D, Hut P. 2003. The Gravitational Million-Body Problem: A Multidisciplinary Ap- 
proach to Star Cluster Dynamics. The Gravitational Million-Body Problem: A Multidisci- 
plinary Approach to Star Cluster Dynamics, by Douglas Heggie and Piet Hut. Cambridge 
University Press, 2003, 372 pp. 

Heggie DC. 1975. MNRAS 173:729-787 

Heggie DC. 1992. Nat 359:772-773 

Heggie DC. 2001. In The Restless Universe, ed. B. A. Steves & A. J. Maciejewski 
Heggie DC, Aarseth SJ. 1992. MNRAS 257:513-536 
Heggie DC, Giersz M. 2008. MNRAS 389:1858-1870 

Heggie DC, Hut P, Mineshige S, Makino J, Baumgardt H. 2007. Publ. Astr. Soc. Japan 
59:L11-L14 

Heggie DC, Inagaki S, McMillan SLW. 1994. MNRAS 271:706 

69 



Henon M. 1965. Annales d'Astrophysique 28:62 

Henon M. 1973. In Saas-Fee Advanced Course 3: Dynamical Structure and Evolution of 
Stellar Systems, ed. G. Contopoulos, M. Henon, & D. Lynden-Bell 

Hillenbrand LA, Hartmann LW. 1998. ApJ 492:540 

Hills JG. 1975. AJ 80:809-825 

Hills JG. 1980. ApJ 235:986-991 

Hills JG, Day CA. 1976. Astrophys. Lett. 17:87 

Ho LC, Filippenko AV. 1996. ApJL 466:L83 

Hodge P. 1987. PASP 99:724-729 

Hodge PW. 1961. ApJ 133:413 

Hodge PW, Krienke OK, Bellazzini M, Perina S, Barmby P, et al. 2009. AJ 138:770-779 

Hoekstra A. G, Portegies Zwart S, Bubak M, Sloot P. 2008. Towards Distributed Petascale 
Computing. Petascale Computing: Algorithms and Applications, by David A. Bader (Ed.). 
Chapman & Hall/CRC computational science series 565pp. (ISBN: 9781584889090, ISBN 
10: 1584889098) 

Holtzman JA, Faber SM, Shaya EJ, Lauer TR, Groth J, et al. 1992. AJ 103:691-702 
Huff EM, Stahler SW. 2006. ApJ 644:355-363 

Hunter DA, Elmegreen BG, Dupuy TJ, Mortonson M. 2003. AJ 126:1836-1848 
Hunter DA, O'Connell RW, Gallagher JS, Smecker-Hane TA. 2000. AJ 120:2383-2401 
Hunter DA, Shaya EJ, Holtzman JA, Light RM, O'Neil Jr. EJ, Lynds R. 1995. ApJ 448:179 
Hurley JR. 2007. MNRAS 379:93-99 

Hurley JR, Aarseth SJ, Shara MM. 2007. ApJ 665:707-718 
Hurley JR, Pols OR, Tout CA. 2000. MNRAS 315:543-569 
Hurley JR, Shara MM. 2002. ApJ 570:184-189 

Hurley JR, Tout CA, Aarseth SJ, Pols OR. 2001. MNRAS 323:630-650 
Hurley JR, Tout CA, Pols OR. 2002. MNRAS 329:897-928 
Hut P, Bahcall JN. 1983. ApJ 268:319-341 
Hut P, Djorgovski S. 1992. Nat 359:806-808 
Hut P, Inagaki S. 1985. ApJ 298:502-520 

70 



Hut P, Makino J, McMillan S. 1988. Nat 336:31-35 

Hut P, Shara MM, Aarseth SJ, Klessen RS, Lombardi Jr. JC, et al. 2003. New Astronomy 
8:337-370 

Inagaki S, Saslaw WC. 1985. ApJ 292:339-347 

Innanen KA, Harris WE, Webbink RF. 1983. AJ 88:338-360 

Jordan A, McLaughlin DE, Cote P, Ferrarese L, Peng EW, et al. 2007. ApJS 171:101-145 
Joshi KJ, Rasio FA, Portegies Zwart S. 2000. ApJ 540:969-982 

Kaaret P, Prestwich AH, Zezas A, Murray SS, Kim DW, et al. 2001. MNRAS 321:L29-L32 

Karakas A, Lattanzio JC. 2007. Publications of the Astronomical Society of Australia 24:103- 
117 

Kennicutt Jr. RC. 1998. ARA&A 36:189-232 
Kennicutt Jr. RC, Chu YH. 1988. AJ 95:720-730 

Kharchenko NV, Piskunov AE, Roser S, Schilbach E, Scholz RD. 2005. A&A 440:403-408 

Kim E, Einsel C, Lee HM, Spurzem R, Lee MG. 2002. MNRAS 334:310-322 

Kim E, Lee HM, Spurzem R. 2004. MNRAS 351:220-236 

Kim SS, Figer DF, Kudritzki RP, Najarro F. 2006. ApJL 653:L113-L116 

King AR. 2002. MNRAS 335:L13-L16 

King AR, Davies MB, Ward MJ, Fabbiano G, Elvis M. 2001. ApJL 552:L109-L112 
King IR. 1962. AJ 67:471 
King IR. 1966. AJ 71:64-75 
Klessen RS. 2001. ApJ 556:837-846 

Klose S, Henden AA, Geppert U, Greiner J, Guetter HH, et al. 2004. ApJL 609:L13-L16 
Knigge C, Leigh N, Sills A. 2009. Nat 457:288-290 

Kouveliotou C, Strohmayer T, Hurley K, van Paradijs J, Finger MH, et al. 1999. ApJL 
510:L115-L118 

Kouwenhoven MBN, de Grijs R. 2008. A&A 480:103-114 

Kroupa P. 1995. MNRAS 277:1491 

Kroupa P. 2001. MNRAS 322:231-246 

71 



Kroupa P. 2008. In Lecture Notes in Physics, Berlin Springer Verlag, ed. S. J. Aarseth, 
C. A. Tout, & R. A. Mardling, vol. 760 of Lecture Notes in Physics, Berlin Springer 
Verlag 

Kulkarni SR, Hut P, McMillan S. 1993. Nat 364:421-423 

Kupi G, Amaro-Seoane P, Spurzem R. 2006. MNRAS 37LL45-L49 

Lada CJ, Lada EA. 2003. ARA&A 41:57-115 

Lada CJ, Margulis M, Dearborn D. 1984. ApJ 285:141-152 

Lamers HJGLM, Gieles M, Bastian N, Baumgardt H, Kharchenko NV, Portegies Zwart S. 
2005. A&A 441:117-129 

Lamers HJGLM, Gieles M, Portegies Zwart SF. 2005. A&A 429:173-179 

Lang CC, Johnson KE, Goss WM, Rodriguez LF. 2005. AJ 130:2185-2196 

Larsen SS. 2002. AJ 124:1393-1409 

Larsen SS. 2004. A&A 416:537-553 

Larsen SS. 2006. In Planets to Cosmology: Essential Science in the Final Years of the Hubble 
Space Telescope, eds. M Livio, S Casertano 

Larsen SS. 2009a. A&A 494:539-551 

Larsen SS. 2009b. Phil. Trans. R. Soc. London A , in press ( arXiv: 0911.0 796) 



Larsen SS, Origlia L, Brodie J, Gallagher JS. 2008. MNRAS 383:263-276 

Larsen SS, Richtler T. 2000. A&A 354:836-846 

Larsen SS, Richtler T. 2004. A&A 427:495-504 

Larsen SS, Richtler T. 2006. A&A 459:103-111 

Larson RB. 1981. MNRAS 194:809-826 

Leonard PJT. 1989. AJ 98:217-226 

Leonard PJT. 1996. ApJ 470:521 

Liu JF, Bregman JN, Irwin J, Seitzer P. 2002. ApJL 581T93-L96 
Liu JF, Bregman JN, Seitzer P. 2004. ApJ 602:249-256 

Lombardi JC, Thrall AP, Deneva JS, Fleming SW, Grabowski PE. 2003. MNRAS 345:762- 
780 

Lucas PW, Roche PF. 2000. MNRAS 314:858-864 

72 



Lutz D. 1991. A&A 245:31-40 

Lynden-Bell D, Eggleton PP. 1980. MNRAS 191:483-498 
Lynden-Bell D, Wood R. 1968. MNRAS 138:495 
Lyne AG, Lorimer DR. 1994. Nat 369:127-129 
Lynga G. 1982. A&A 109:213-222 

Ma J, Fan Z, de Grijs R, Wu Z, Zhou X, et al. 2009. AJ 137:4884-4896 
Mackey AD, Broby Nielsen P. 2007. MNRAS 379:151-158 
Mackey AD, Gilmore GF. 2003. MNRAS 338:85-119 

Mackey AD, Wilkinson MI, Davies MB, Gilmore GF. 2007. MNRAS 379:L40-L44 
Mackey AD, Wilkinson MI, Davies MB, Gilmore GF. 2008. MNRAS 386:65-95 
Madhusudhan N, Justham S, Nelson L, Paxton B, Pfahl E, et al. 2006. ApJ 640:918-922 
Maiz-Apellaniz J. 2001. ApJ 563:151-162 

Maiz-Apellaniz J. 2002. In Extragalactic Star Clusters, eds. DP Geisler, EK Grebel, D Min- 
niti, vol. 207 of IAU Symposium 

Makino J. 1996. ApJ 471:796 

Makino J. 2002. New Astronomy 7:373-384 

Makino J. 2005. ArXiv Astrophysics e-prints 

Makino J, Aarseth SJ. 1992. Publ. Astr. Soc. Japan 44:141-151 

Makino J, Hut P, Kaplan M, Saygm H. 2006. New Astronomy 12:124-133 

Makino J, Taiji M. 1998. Scientific simulations with special-purpose computers : The GRAPE 
systems. Scientific simulations with special-purpose computers : The GRAPE systems /by 
Junichiro Makino & Makoto Taiji. Chichester ; Toronto : John Wiley & Sons, cl998. 

Makishima K, Kubota A, Mizuno T, Ohnishi T, Tashiro M, et al. 2000. ApJ 535:632-643 

Marco A, Negueruela I, Motch C. 2007. In Massive Stars in Interactive Binaries, eds. N St.- 
Louis, AFJ Moffat, vol. 367 of Astronomical Society of the Pacific Conference Series 

Mason BD, Hartkopf WI, Gies DR, Henry TJ, Helsel JW. 2009. AJ 137:3358-3377 

Massey P, Hunter DA. 1998. ApJ 493:180 

Mayor M, Udry S, Lovis C, Pepe F, Queloz D, et al. 2009. A&A 493:639-644 
McCrady N, Gilbert AM, Graham JR. 2003. ApJ 596:240-252 

73 



McCrady N, Graham JR. 2007. ApJ 663:844-856 
McCrady N, Graham JR, Vacca WD. 2005. ApJ 621:278-284 
McCrea WH. 1964. MNRAS 128:147 
McLaughlin DE. 2000. ApJ 539:618-640 

McLaughlin DE, Anderson J, Meylan G, Gebhardt K, Pryor C, et al. 2006. ApJS 166:249- 
297 

McLaughlin DE, van der Marel RP. 2005. ApJS 161:304-360 
McMillan S, Hut P, Makino J. 1990. ApJ 362:522-537 
McMillan S, Hut P, Makino J. 1991. ApJ 372:111-124 

McMillan SLW. 1986. In The Use of Supercomputers in Stellar Dynamics, ed. P. Hut & 
S. L. W. McMillan, vol. 267 of Lecture Notes in Physics, Berlin Springer Verlag 

McMillan SLW, Aarseth SJ. 1993. ApJ 414:200-212 

McMillan SLW, Vesperini E, Portegies Zwart SF. 2007. ApJL 655:L45-L49 

Mengel S, Lehnert MD, Thatte N, Genzel R. 2002. A&A 383:137-152 

Mengel S, Lehnert MD, Thatte NA, Vacca WD, Whitmore B, Chandar R. 2008. A&A 
489:1091-1105 

Mengel S, Tacconi-Garman LE. 2007. A&A 466:151-155 

Merritt D, Piatek S, Portegies Zwart S, Hemsendorf M. 2004. ApJL 608:L25-L28 
Meszaros P. 2002. ARA&A 40:137-169 

Meurer GR, Heckman TM, Leitherer C, Kinney A, Robert C, Garnett DR. 1995. AJ 
110:2665 

Miller BW, Whitmore BC, Schweizer F, Fall SM. 1997. A J 114:2381 
Miller GE, Scalo JM. 1979. ApJS 41:513-547 

Miller JM, Fabbiano G, Miller MC, Fabian AC. 2003. ApJL 585:L37-L40 

Miller MC, Hamilton DP. 2002. MNRAS 330:232-240 

Milone AP, Bedin LR, Piotto G, Anderson J. 2009. A&A 497:755-771 

Milone AP, Piotto G, Bedin LR, Sarajedini A. 2008. Memorie della Societa Astronomica 
Italiana 79:623 



74 



Mirabel IF, Fuchs Y, Chaty S. 2000. In Gamma-ray Bursts, 5th Huntsville Symposium, ed. 
R. M. Kippen, R. S. Mallozzi, & G. J. Fishman, vol. 526 of American Institute of Physics 
Conference Series 

Mirabel IF, Rodrigues I, Liu QZ. 2004. A&A 422:L29-L32 

Moeckel N, Bonnell IA. 2009. MNRAS 400:657-664 

Moll SL, Mengel S, de Grijs R, Smith LJ, Crowther PA. 2007. MNRAS 382:1877-1888 

Moore B, Quinn T, Governato F, Stadel J, Lake G. 1999. MNRAS 310:1147-1152 

Muno MP, Clark JS, Crowther PA, Dougherty SM, de Grijs R, et al. 2006. ApJL 636:L41- 
L44 

Muno MP, Gaensler BM, Clark JS, de Grijs R, Pooley D, et al. 2007. MNRAS 378:L44-L48 
Narayan R, Piran T, Shemi A. 1991. ApJL 379:L17-L20 

Narbutis D, Vansevicius V, Kodaira K, Bridzius A, Stonkute R. 2008. ApJS 177:174-180 

Nitadori K, Makino J. 2008. New Astronomy 13:498-507 

Nitadori K, Makino J, Hut P. 2006. New Astronomy 12:169-181 

Normandeau M, Taylor AR, Dewdney PE. 1996. Nat 380:687-689 

Offner SSR, Hansen CE, Krumholz MR. 2009. ApJL 704:L124-L128 

Ogilvie GI, Lin DNC. 2004. ApJ 610:477-509 

Oort JH, Kerr FJ, Westerhout G. 1958. MNRAS 118:379 

Ostlin G, Cumming RJ, Bergvall N. 2007. A&A 461:471-483 

Ostriker JP, Spitzer LJ, Chevalier RA. 1972. ApJL 176:L51 

Paczynski B. 1998. ApJL 494:L45 

Park W, Park HS, Lee MG. 2009. ApJ 700:103-113 

Parker RJ, Goodwin SP. 2007. MNRAS 380:1271-1275 

Parmentier G, de Grijs R. 2008. MNRAS 383:1103-1120 

Parmentier G, Gilmore G. 2007. MNRAS 377:352-372 

Patruno A, Portegies Zwart S, Dewi J, Hopman C. 2006. MNRAS 370:L6-L9 
Peacock MB, Maccarone TJ, Knigge C, Kundu A, Waters CZ, et al. 2010. MNRAS :28 
Perina S, Barmby P, Beasley MA, Bellazzini M, Brodie JP, et al. 2009. A&A 494:933-948 



75 



Pfahl E, Rappaport S, Podsiadlowski P. 2002. ApJ 573:283-305 
Pfahl E, Scannapieco E, Bildsten L. 2009. ArXiv e-prints 
Pfalzner S. 2009. A&A 498:L37-L40 

Piotto G. 2008. Memorie della Societa Astronomica Italiana 79:334 
Piotto G, Villanova S, Bedin LR, Gratton R, Cassisi S, et al. 2005. ApJ 621:777-784 
Piskunov AE, Kharchenko NV, Roser S, Schilbach E, Scholz RD. 2006. A&A 445:545-565 
Piskunov AE, Schilbach E, Kharchenko NV, Roser S, Scholz RD. 2008. A&A 477:165-172 
Plummer HC. 1911. MNRAS 71:460-470 

Pooley D, Lewin WHG, Anderson SF, Baumgardt H, Filippenko AV, et al. 2003. ApJL 
591:L131-L134 

Portegies Zwart SF, Baumgardt H, Hut P, Makino J, McMillan SLW. 2004. Nat 428:724- 
726 

Portegies Zwart SF, Belleman RG, Geldof PM. 2007. New Astronomy 12:641-650 

Portegies Zwart SF, Hut P, Verbunt F. 1997. A&A 328:130-142 

Portegies Zwart SF, Makino J, McMillan SLW, Hut P. 1999. A&A 348:117-126 

Portegies Zwart SF, McMillan S, Groen D, Gualandris A, Sipior M, Vermin W. 2008. New 
Astronomy 13:285-295 

Portegies Zwart SF, McMillan S, Harfst S, Groen D, Fujii M, et al. 2009. New Astronomy 
14:369-378 

Portegies Zwart SF, McMillan SLW. 2002. ApJ 576:899-907 

Portegies Zwart SF, McMillan SLW, Hut P, Makino J. 2001. MNRAS 321:199-226 

Portegies Zwart SF, McMillan SLW, Makino J. 2007. MNRAS 374:95-106 

Portegies Zwart SF, Pooley D, Lewin WHG. 2002. ApJ 574:762-770 

Portegies Zwart SF, Takahashi K. 1999. Celestial Mechanics and Dynamical Astronomy 
73:179-186 

Portegies Zwart SF, van den Heuvel EPJ. 2007. Nat 450:388-389 
Portegies Zwart SF, Verbunt F. 1996. A&A 309:179-196 
Prantzos N, Charbonnel C. 2006. A&A 458:135-149 
Price DJ, Bate MR. 2009. MNRAS 398:33-46 

76 



Proszkow EM, Adams FC, Hartmann LW, Tobin J J. 2009. ApJ 697:1020-1031 
Quinlan GD. 1996. New Astronomy 1:255-270 

Rauw G, Crowther PA, De Becker M, Gosset E, Naze Y, et al. 2005. A&A 432:985-998 
Read JI, Wilkinson MI, Evans NW, Gilmore G, Kleyna JT. 2006. MNRAS 366:429-437 
Roberts DA, Yusef-Zadeh F. 2005. AJ 129:805-808 

Sabbi E, Sirianni M, Nota A, Tosi M, Gallagher J, et al. 2008. AJ 135:173-181 
Safronov VS. 1969. Evoliutsiia doplanetnogo oblaka. 
Salpeter EE. 1955. ApJ 121:161 

San Roman I, Sarajedini A, Garnett DR, Holtzman JA. 2009. ApJ 699:839-849 

Sana H, Gosset E, Naze Y, Rauw G, Linder N. 2008. MNRAS 386:447-460 

Sandage AR. 1953. AJ 58:61-75 

Sarajedini A, Mancone CL. 2007. AJ 134:447-456 

Scally A, Clarke C, McCaughrean MJ. 2005. MNRAS 358:742-754 

Schechter P. 1976. ApJ 203:297-306 

Scheepmaker RA, Haas MR, Gieles M, Bastian N, Larsen SS, Lamers HJGLM. 2007. A&A 
469:925-940 

Schilbach E, Roser S. 2008. A&A 489:105-114 

Schweizer F. 1982. ApJ 252:455-460 

Shapiro SL. 1985. In Dynamics of Star Clusters, ed. J. Goodman & P. Hut, vol. 113 of IAU 

Symposium 

Shara MM, Hurley JR. 2002. ApJ 571:830-842 

Silich S, Tenorio-Tagle G, Rodriguez-Gonzalez A. 2004. ApJ 610:226-232 

Sills A, Deiters S, Eggleton P, Freitag M, Giersz M, et al. 2003. New Astronomy 8:605-628 

Smith LJ, Gallagher JS. 2001. MNRAS 326:1027-1040 

Smith LJ, Westmoquette MS, Gallagher JS, O'Connell RW, Rosario DJ, de Grijs R. 2006. 
MNRAS 370:513-527 

Sollima A, Beccari G, Ferraro FR, Fusi Pecci F, Sarajedini A. 2007. MNRAS 380:781-791 
Sollima A, Lanzoni B, Beccari G, Ferraro FR, Fusi Pecci F. 2008. A&A 481:701-704 

77 



Sommariva V, Piotto G, Rejkuba M, Bedin LR, Heggie DC, et al. 2009. A&A 493:947-958 

Spitzer L. 1987. Dynamical evolution of globular clusters. Princeton, NJ, Princeton Univer- 
sity Press, 1987, 191 p. 

Spitzer Jr. L. 1975. In Dynamics of the Solar Systems, ed. A. Hayli, vol. 69 of IAU Symposium 

Spitzer LJ. 1940. MNRAS 100:396 

Spitzer LJ. 1958. ApJ 127:17 

Spitzer LJ. 1969. ApJL 158:L139 

Spitzer LJ, Chevalier RA. 1973. ApJ 183:565-582 

Spitzer LJ, Hart MH. 1971. ApJ 164:399 

Spurzem R. 1999. Journal of Computational and Applied Mathematics 109:407-432 
Spurzem R, Giersz M, Heggie DC, Lin DNC. 2009. ApJ 697:458-482 
Stodolkiewicz JS. 1982. Acta Astronomica 32:63-91 
Stodolkiewicz JS. 1986. Acta Astronomica 36:19-41 

Stolte A, Brandner W, Grebel EK, Lenzen R, Lagrange AM. 2005. ApJL 628:L113-L117 

Stolte A, Grebel EK, Brandner W, Figer DF. 2002. A&A 394:459-478 

Strohmayer TE, Mushotzky RF. 2003. ApJL 586:L61-L64 

Swartz DA, Tennant AF, Soria R. 2009. ApJ 703:159-168 

Takahashi K. 1996. Publ. Astr. Soc. Japan 48:691-700 

Takahashi K. 1997. Publ. Astr. Soc. Japan 49:547-560 

Takahashi K, Portegies Zwart SF. 1998. ApJL 503:L49 

Takahashi K, Portegies Zwart SF. 2000. ApJ 535:759-775 

Theuns T. 1991. Memorie della Societa Astronomica Italiana 62:909-914 

Thorsett SE, Arzoumanian Z, Taylor JH. 1993. ApJL 412T33-L36 

Trenti M, Ardi E, Mineshige S, Hut P. 2007. MNRAS 374:857-866 

Udry S, Santos NC. 2007. ARA&A 45:397-439 

Subr L, Kroupa P, Baumgardt H. 2008. MNRAS 385:1673-1680 

Vallee JP. 2008. AJ 135:1301-1310 



78 



van den Bergh S. 1957. ApJ 125:445 

van den Bergh S. 1971. A&A 12:474 

van den Bergh S. 1991. ApJ 369:1-12 

van den Bergh S, McClure RD. 1980. A&A 88:360-362 

van den Berk J, Portegies Zwart SF, McMillan SLW. 2007. MNRAS 379:111-122 

Vanbeveren D, Belkus H, van Bever J, Mennekens N. 2009. Ap&SS :114 

Vansevicius V, Kodaira K, Narbutis D, Stonkute R, Bridzius A, et al. 2009. ApJ 703:1872- 
1883 

Ventura P, D'Antona F, Mazzitelli I, Gratton R. 2001. ApJL 550:L65-L69 
Vesperini E. 1997. MNRAS 287:915-928 

Vesperini E, D'Ercole A, D'Antona F, McMillan SLW, Recchi S. 2009. In Bulletin of the 
American Astronomical Society, vol. 41 of Bulletin of the American Astronomical Society 

Vesperini E, Heggie DC. 1997. MNRAS 289:898-920 

Vesperini E, McMillan SLW, Portegies Zwart S. 2009. ApJ 698:615-622 

Vesperini E, Zepf SE. 2003. ApJL 587:L97-L100 

Vesperini E, Zepf SE, Kundu A, Ashman KM. 2003. ApJ 593:760-771 

Vinko J, Sarneczky K, Balog Z, Immler S, Sugerman BEK, et al. 2009. ApJ 695:619-635 

von Hoerner S. 1957. ApJ 125:451 

von Hoerner S. 1958. Zeitschrift fur Astrophysik 44:221-242 
von Hoerner S. 1963. Zeitschrift fur Astrophysik 57:47-82 

Vrba FJ, Henden AA, Luginbuhl CB, Guetter HH, Hartmann DH, Klose S. 2000. ApJL 
533:L17-L20 

Wang X, Yang Y, Zhang T, Ma J, Zhou X, et al. 2005. ApJL 626:L89-L92 

Weidner C, Kroupa P. 2004. MNRAS 348:187-191 

Weidner C, Kroupa P. 2006. MNRAS 365:1333-1347 

Weidner C, Kroupa P, Larsen SS. 2004. MNRAS 350:1503-1510 

Weigelt G, Baier G. 1985. A&A 150:L18-L20 

Weinberg MD. 1994. AJ 108:1403-1413 

79 



Weldrake DTF, Sackett PD, Bridges TJ, Freeman KC. 2005. ApJ 620:1043-1051 
White NE, van Paradijs J. 1996. ApJL 473:L25 

Whitmore BC. 2003. In A Decade of Hubble Space Telescope Science, eds. M Livio, K Noll, 
M Stiavelli 

Whitmore BC, Chandar R, Fall SM. 2007. A J 133:1067-1084 
Whitmore BC, Schweizer F. 1995. AJ 109:960-980 

Whitmore BC, Zhang Q, Leitherer C, Fall SM, Schweizer F, Miller BW. 1999. A J 118:1551- 
1576 

Wielen R. 1971. A&A 13:309-322 

Wielen R. 1988. In IAU Symp. 126: The Harlow-Shapley Symposium on Globular Cluster 
Systems in Galaxies, eds. JE Grindlay, AGD Philip 

Wilkinson MI, Hurley JR, Mackey AD, Gilmore GF, Tout CA. 2003. MNRAS 343:1025- 
1037 

Wu H, Xue SJ, Xia XY, Deng ZG, Mao S. 2002. ApJ 576:738-744 
Wyder TK, Hodge PW, Zucker DB. 2000. PASP 112:1162-1176 

Yungelson LR, van den Heuvel EPJ, Vink JS, Portegies Zwart SF, de Koter A. 2008. A&A 
477:223-237 

Zepf SE, Ashman KM, English J, Freeman KC, Sharpies RM. 1999. AJ 118:752-764 
Zezas A, Fabbiano G, Rots AH, Murray SS. 2002. ApJ 577:710-725 
Zhang Q, Fall SM. 1999. ApJL 527:L81-L84 

A Appendix: Dynamical algorithms 

Even though the fundamental physics is not hard to understand, simulating star clusters 
is not a trivial matter. Significant complications arise due to the long-range nature of 
the gravitational force, which means that every star in the cluster is effectively in constant 
communication with every other. The result is a potentially enormous number of interactions 
that must be calculated as we follow the time evolution of the system, leading to high 
computational cost. Further complications arise from the enormous range in spatial and 
temporal scales inherent in a star cluster. Computers, by the way they are constructed, have 
difficulty in resolving such wide ranges, and many of the software problems in simulations of 
self-gravitating systems arise from this basic limitation. The combination of many physical 
processes occurring on many scales, with high raw processing requirements, makes numerical 
gravitational dynamics among the most demanding and challenging areas of computational 
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science. In this Appendix we discuss some of the issues involved in the numerical modeling 
of massive star clusters. 

The broad spectrum of numerical methodologies currently available is summarized in 
§4.11 Here we present some of the details of these simulation techniques, starting with 
"Continuum" Models, followed by Monte Carlo Models and Direct N-body Models. We 
end this section with brief discussions of computer hardware and multi-core codes, and 
in particular Kitchen Sink Models and Future Prospects. 



A.l Continuum methods 



The two leading classes of continuum models are gas-sp here (jLynden-Bell fc Eggleton 1 980. 



Bettwieser fc Sugimoto 19841 iDeiters fc Spurzem 20011) and Fokker-Planck (ICohn 1 979; Sh apiro 1985 . 
Chernoff fc Weinberg 1990l : lDrukier. Fahlman fc Richer 19921 iTakahashi 19961 . ll997l : lTakahashi fc Portegi 



1C 



methods. They have mainly been applied to spherically symmetric s ystems, although axisym- 



metric extensions to rotating systems have a lso been implemented (jEinsel fc Spurzem 1999 



Kim et al. 20021 ; iKim. Lee Spurzem 20041 ). and some limited ex periments with rudimen- 
tary binary treatments have also been performed (IGao et al. 19911 ) . 

Both approaches start with the collisional Boltzmann equation as the basic description 
for a stellar system, then simplify it by averaging the distribution function /(x, v) in different 
ways. Gas-sphere methods proceed in a manner closely analogous to the derivation of the 
equations of fluid motion, taking velocity averages to construct the moments of the distri- 
bution: p = f d 3 v /(x, v), u = f d 3 v v/(x, v), a 2 = ~ f d 3 v v 2 /(x, v), etc. Application of a 
closure condition leads to a set of equations identical to those of a classical conducting fluid, 
in which the conductivity depends inversely on the local relaxation time. Fokker-Planck 
methods transform the Boltzmann equation by orbit-averaging all quantities and recasting 
the equation as a diffusion equation in E — J space, where E is stellar energy and J is angu- 
lar momentum. Since both E and J are conserved orbital quantities in a static, spherically 
symmetric system, two-body relaxation enters into the problem via the diffusion coefficients. 

These methods have been of enormous value in developing and refining theoretical insights 
i nto the fundamental physical processes driving the dynamical evolution of stellar systems 
(IBettwieser fc Sugimoto 19841 ). However, as the degree of realism demanded of the simu- 
lation increases — adding a mass spectrum, stellar evolution, binari es, etc. — the algorithms 



rapidly become cumbersome, inefficient, and of questionable validity (IPortegies Zwart fc Takahashi 19991 ) . 
As a result, they are generally not applied to the young stellar systems of interest here. The 
major approaches currently used for simulating young massive clusters are particle-based 
Monte Carlo or direct iV-body codes. 



A. 2 Monte Carlo methods 

Depending on one's point of view, Monte Carlo methods can be regarded as particle al- 
gorithms for solving the partial differential equations arising from the continuum models, 
or approximate schemes for determining the long-term average gravitational interactions 
of a large collection of particles. The early techniques developed in the 1970s and 1980s 
flSpitzer k Hart 197ll . iHenon 19731 . ISpitzer 19751. IStodolkiewicz 1981 fl98~6h fall into the for- 
mer category, but recent studies, in particular ( Giersz 19981 : Joshi. Rasio &: Portegies Zwart 200ol : 
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Freitae & Benz 2001: 
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treatment of the "background" stellar population with a Monte Carlo realization of the or- 
bits and interactions of binaries and other objects of interest. These approaches have allowed 
the first simulations of an entire globular cluster, from a very early (although gas depleted) 
phase to complete dissolution. 

Monte Carlo methods are designed for efficient computation of relaxation effects in colli- 
sional stellar systems, a task which they accomplish by reducing stellar orbits to their orbital 
elements — energy and angular momentum — effectively orbit averaging the motion of each 
star. Relaxation is modeled by randomly selecting pairs of stars and applying interactions 
between them in such a way that, on average, the correct rate is obtained. This may be imple- 
mented in a number of ways, but interactions are generally realized on time scales comparable 
to the orbit-averaged relaxation time. As a result, Mo nte Carlo schemes can be orders of mag - 



nitude faster than direct iV-body codes. For example, iJoshi. Rasio fc Portegies Zwart (20001 ) 
report a CPU time scaling for their Monte Carlo scheme of 0(N 1A ) for core-collapse prob- 
lems, compared to iV 3 for iV-body methods, as discussed below. To achieve these speeds, 
however, the geometry of the system must be simple enough that the orbital integrals can be 
computed from a star's instantaneous energy and angular momentum. In practice, this limits 
the approach to spherically symmetric systems in virial equilibrium, and global dynamical 
processes occurring on relaxation (or longer) time scales. 



A. 3 iV-body methods 

iV-body codes incorporate detailed descriptions of stellar dynamics at all levels, using direct 



integration of the in dividual (Newtonian) stellar equations of motion for all stars (lAarseth 2003 



Heggie &: Hut 20031 ). Their major attraction is that they are assumption- free, in the sense 
that all stellar interactions are automatically included to all orders, without the need for any 
simplifying approximations or the inclusion of additional reaction rates to model particular 
physical processes of interest. Thus, problems inherent to Monte Carlo methods (see §A.2j) . 
related to departures from virial equilibrium, spherical symmetry, statistical fluctuations, the 
form of (and indeed the existence of) phase space distribution functions, and the possibility 
of interactions not explicitly coded in advance, simply do not arise, and therefore do not 
require fine-tuning as in the Monte Carlo models. 

The price of all these advantages is computational expense. Each of the N particles must 
interact with every other particle a few hundred times over the course of every orbit, each 
interaction requires O(N) force calculations, and a typical (relaxation time) run spans O(N) 
orbits (see Eq. [TBI) . The resulting 0(N 3 ) scaling of the total CPU time means that, even 
with the best time-step algorithms, integrating even a fair ly small system of, say, N ~ 10 5 
stars requires sustained teraflops speeds for several months (IHut. Makino fc McMillan 19881 ). 
Radically improved performance can be achieved by writing better software, or by building 
faster computers (or both). In fact, the remarkable speed-up of iV-body codes over the last 
four decades has mainly been due to advances in hardware, and in a lesser extend due to 
software. 

Substantial performance improvements were realized by adopting better (individual) time 
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stepping schemes (as opposed to earlier shared time step schemes), in which particles ad- 
vance using steps appropriate to their individual orbits, r ather than a single step for all. 
Further gains were made by utilizing neighbor schemes ( Ahmad fc Cohen 1973bh . which 
divide the force on every particle into irregular (rapidly varying) and regular (slowly vary- 
ing) parts, due (loosely speaking) to nearby and more distant bodies. By recomputing the 
regular force at every particle step, but extrapolating the more expensive O(N) regular 
force for most time steps, and recomputing it only on longer time scales, significant im- 
pr qvements in efficiency have been realize d. A multi-level generalization of this approach 
by iDorband. Hemsendorf fc Merritt (20031 ) is incorporated into the collisonal NBODY6++ 
(ISpurzem 19991 ) . 

Another important alg orithmic improvemen t was introduced in the mid-1980s with the 
development of tree codes (jBarnes fc Hut 19861 ). which reduce the force calculation complex- 
ity from O(N) to O(logA^). Despite their algorith mic efficiency, tree codes h ave not been 
widely used in modeling collisional systems (but see iMcMillan fc Aarseth 19931 ). This seems 
principally to be because of lingering technical concerns about their long-term accuracy in 
systems dominated by relaxation processes and their performance in clusters with large dy- 
namic ranges in densities and time scales, even though these objections may no longer be 
well founded flMoore et al. 19991 . iDehnen 200oh . Very promising direct-treecode methods 
have recently been developed to model the dynamical interaction between a cl uster and the 
surrounding galactic population (IFujii et al. 20071 . iPortegies Zwart et al. 2009l ) . 



A. 4 Parallelization 



Individual time step schemes are general ly hard to optimize on parallel mach ines. For those 
architectures, block time step schemes (IMcMillan 19861 . iMakino et al. 20061 ) offer substan- 
tially better performance. By rounding each star's "natural" step down to the nearest neg- 
ative integer power of two, such a scheme effectively discretizes the time variable, allowing 
the possibility that large blocks of stars will be "next" on the time step list, and so can be 
efficiently integrated in parallel. 

The two most important parallel integration techniques are the ring and copy algo- 
rithms. Both have advantages and disadvantages, but the execution times for each, on 
computers with p processors, scale as N/p, while the communication times scale as p 
(IHarfst et al. 20071 ). Both algorithms are implemented in a range of A^-body codes, in- 
cluding N BODY6++ (IDorband. Hemsendorf fc Merritt 20031 ) and the kira integrator in 
Starlab ( Portegies Zwart et al. 20081 ). The two-dimensional lattice parallelization for di- 
rect A^-body kernels has comparable CPU time scaling, but the communication has a weaker 
scaling (oc l/^/p), enabli ng the code to maintain satisf actory performance even on computers 
with p ^ 10 3 processors ( Makino 2002 . Bisseling 20041 ). So far, however, this scheme has not 
been implemented in a production A^-body code. 

An interesting further step is to use a widely distributed grid of computers ( Foster fc Kesselman 20041 ). 
In this extreme form of parallel computing the computational bottleneck often shifts from 
the 0(N 2 ) force calculation (see §A.3j) to communication (latency and bandwidth). However, 
even in the worst-case scenario the communication costs scale oc N, so, for a sufficiently larg e 
number of stars even intercontinental grid computing can be practical ( iHoekstra et al. 20081 ). 
In addition, if (as seems likely — see §A.6I) future simulation environments will combine a 
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range of codes in addition to pure stellar dynamics to address the evolution of YMCs in 
detail, grid computing may provide the solution to the problem of limited supply of local 
computer resources. This is particularly relevant if the desired algorithms for solving stellar 
dynamics, stellar evolution, hydrodynamics, etc, require a diversity in computer architectures 
that may not be locally available. 



A. 5 Hardware acceleration 



A quantum leap in gravitational iV-body simulation speed came from the introduction of 
special purpose computers. All iV-body codes, including neighbor schemes and treecodes, 
suffer from the cost of computing inter-particle forces at every step along the orbit. A tech- 
nological solution in widespread use is the "GRAPE" (short for "G RAvity PipE") series o f 
machines developed by Sugimoto and co-workers at Tokyo University (jEbisuzaki et al. 19931 ). 
Abandoning algorithmic sophistication in favor of simplicity and raw computing power, 
these m achines achieved high pe rformance by mating a fourth-order Hermite integration 
scheme (IMakino fc Aarseth 19921 ) with special-purpose hardware in the form of highly par- 
allel, pipelined "Newtonian force accelerators" implementing the computation of all inter- 
particle forces entirely in hardware. Operationally, the GRAPE hardware is simple to pro- 
gram, as it merely replaces the function that computes the (regular) force on a particle 
by a call to the hardware interface libraries; the remainder of the user's A-body code is 
unchanged. 

The effect of GRAPE on simulations of stellar systems has been nothing short of revolu- 
tionary. Today, GRAPE-enabled code lies at the heart of almost all detailed A-body simula- 
tions of star clusters and dense stellar systems. GRAPE-like optimizations of the innermost 
force calculation operations using the IA-32 Streamin g SIMP Extensions 2 (SSE2), su itable 
for use on any Intel/AMD processor, are described by lNitadori. Making &: Hut (20061 ). 

Recently, Graphics Processing Units (GPUs) have ac hieved speeds and price /performance 

levels previously attainable only by GRAPE systems fseelPortegies Zwart. Belleman fc Geldof 20071: 

Hamada fc Iitaka 20071 ; iBelleman. Bedorf &: Portegies Zwart 20081 ; iGaburov. Harfst fc Portegies Zwart 20( 
for recent GPU implementations of the GRA PE interface). In addition, the programming 
model for GPUs (as well as the GRAPE-DR, IMakino 20051 ). means that many other kinds 
of algorithms can (in principle) be accelerated, although, in practice, it currently seems that 
CPU-intensive operations such as direct A-body force summation show substantially better 
acceleration than, say, treecodes running on the same hardware. It appears that commodity 
components may be poised to outpace special-purpose computers in this specialized area of 
computational science, just as they have already done in general-purpose computing. 



A. 6 The kitchen sink 



Consistent with our growing understanding of the role of stellar and binary interactions in 
collisional stellar systems, the leading programs in this field are "kitchen sink" packages 
that combine treatments of dynamics, stellar and binary evolution, and stellar hydrody- 
namics within a single simulation. Of t hese, the most widely used are the iV-body codes 



NBODY flHurlev et al. 2001 



Aarseth 20031 ). kira which is part of the starlab package (e.g . 



Portegies Zwart et al. 20011 ) . and the Monte Carlo codes developed by Giersz (IGiersz 19981 . 
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Heggie fc Giersz 2008llGiersz fc Heggie 20091 ). Freitag ( IFregeau et al. 20031 ; Freitag. Rasio fc Baumgardt 2i 
and Fregeau (IFregeau fc Rasio 20071 ). 

Despite the differences in their handling of the large-scale dynamics, as just outlined, 
these codes all employ conceptually similar approaches to stellar and binary evolution and 
collisions. All use approximate descriptions of stellar evolution, generally derived from look- 
up tables based on the deta iled evolutionary models of lEggleton. Fitchett fc Tout (19891 ) and 
Hurley. Pols &: Tout (20001). They also rely on semi-analytic or heuristic rule-based treat- 
ments of binary evolution ( IPortegies Zwart Verbunt 1996c iHurley. Tout fc Pols 20021 ) , con- 
ceptually similar from code to code, but significantly different in detail. 

In most cases, collisions are implemented in the simple "sticky-sphere" approximation, 
where stars are taken to collide (and merge) if they approach within the sum of their effec- 
tive radii. The effective radii may be calibrated using hydrodynamical simulations, and mass 
loss may be included in some approximate way. Freitag's Monte Carlo code, geared mainly 
to studies of galactic nuclei, uses a more sophisticated approach, interpolating encounter 
outco mes from a pre-comp uted grid of smoothed particles hydrodynamics (SPH) simula- 
tions (IFreitag &: Benz 20051 ). An interesting alternative, thoug h currently only opera tional 
in AMUSE (see is the "Make Me A Star" pac kage fMMAS:lLombardi et al. 2003^ and 

its extension "Make Me a Massive Star" (MMAMS ^Gaburov. Lombardi fc Portegies Zwart 2008l fl. 
MMA(M)S constructs a merged stellar model by sorting the fluid elements of the original 
stars by entropy or density, then recomputing their equilibrium configuration, using mass 
loss and shock heating data derived from SPH calculations. 

Small-scale dynamics of multiple stellar encounters, such as binary and higher-order en- 
counters, are often handled by look-up from pre-computed cross sections or — more commonly — 
by direct integration, either in isolation or as part of a larger iV-body calculation. Codes 
employing direct integr ation may also include post-Newtonian terms in the interactions be- 
tween compact objects (IKupi. Amaro-Seoane fc Spurzem 20061 ). 



A. 7 Future prospects 

The very comprehensiveness of kitchen-sink codes gives them the great advantage of appli- 
cability to complex stellar systems, but also the significant disadvantage of inflexibility. By 
selecting one of these codes, one is implicitly choosing a particular hard-coded combination 
of dynamical integrator, stellar and binary evolution schemes, collision prescription, and 
treatment of multiple dynamics. The structure of these codes is such that implementing a 
different algorithm within the larger framework is difficult at best, and practically impossible. 

However, studies of dense stellar systems force interactions between programs that were 
never intended to interact with other programs, and by extension require new communication 
channels between the programmers responsible for them. Closely rel ated to this effort is the 
"MUSE" (MUltiscale, MUltiphysics Software Environment) proiectPl Jportegies Zwart et al. 2009h . 
and its successor AMUSE (Astrophysical Multipurpose Software Environment), two am- 
bitious open-source efforts in code integration. (A) MUSE aims at the self-consistent in- 
tegration of dynamics, collisions, stellar evolution, and other relevant physical processes, 



h See http : //webpub . allegheny . edu/employee/ j / j alombar/mmas/ 
'See http : / / castle . strw . leidenuniv . nl/ 



J http : //muse . li 
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thereby realizing o ne vision of the MODEST0 community (IHut et al. 20031 . ISills et al. 20031 



Davies et al. 20061 ). The long-term goal is a comprehensive environment for modeling dense 



stellar systems, including multiphysics/legacy codes and flexible interfaces to integrate ex- 
isting software (written in many languages) within this unifying environment. 

Aside from future developments in modular simulation environments, the recent appear- 
ance of programmable high-performance hardware has spurred the develo pment of new algo- 

rithms for implementing the force-evaluation operations in N-body codes. ( INitadori fc Making 20081 ) 
have developed extensions of the standard fourth-order Hermite scheme to higher (sixth and 
eighth) orders, and Gaburov and Nitadori (2010, in pr eparation) have incorporated these 
methods into the well-known [Ahmad fc Cohen (1973al ) neighbor scheme. Other intrigu- 
ing future prospects come from hybridization of direct methods with hierarchical tree and 
particle-mesh algorithms, an approach currently under development by Nitadori and Makino 
(2009, private communication), but not yet operational. 



k MODEST stands for MOdeling DEnse STellar systems, and can be found at 
tp : / / www . manybody . org/mode st . 
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This figure "a.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/1002.1961vl 



